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Abstract 

In this thesis, I investigate aspects of local Hamiltonians in quantum computing. First, 
I focus on the Adiabatic Quantum Computing model, based on evolution with a time- 
dependent Hamiltonian. I show that to succeed using AQC, the Hamiltonian involved must 
have local structure, which leads to a result about eigenvalue gaps from information the- 
ory. I also improve results about simulating quantum circuits with AQC. Second, I look 
at classically simulating time evolution with local Hamiltonians and finding their ground 
state properties. I give a numerical method for finding the ground state of translationally 
invariant Hamiltonians on an infinite tree. This method is based on imaginary time evo- 
lution within the Matrix Product State ansatz, and uses a new method for bringing the 
state back to the ansatz after each imaginary time step. I then use it to investigate the 
phase transition in the transverse field Ising model on the Bethe lattice. Third, I focus on 
locally constrained quantum problems Local Hamiltonian and Quantum Satisfiability and 
prove several new results about their complexity. Finally, I define a Hamiltonian Quantum 
Cellular Automaton, a continuous-time model of computation which doesn't require control 
during the computation process, only preparation of product initial states. I construct two 
of these, showing that time evolution with a simple, local, translationally invariant and 
time-independent Hamiltonian can be used to simulate quantum circuits. 
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Chapter 

Introduction 



Today's computers, working in binary, churn terabytes of classical information for us every 
day, solving computational problems. Nevertheless, for some of these problems, the required 
resources for even the best algorithms available today scale unfavorably with the problem 
size. Such computational difficulty can be turned into a resource of its own, motivating 
cryptographic schemes. Today's favorite, RSA, is based on the hardness of factoring large 
numbers. It can withstand attacks by even bigger and faster computers, as making the 
input size of the problem larger increases the required resources exponentially. Is there 
then a chance of ever solving hard problems? For many optimization problems like 3- 
SAT, the expected answer is no. There are problems thought to be easier than 3-SAT, 
such as graph isomorphism, the problem to determine whether two graphs are isomorphic 
or not. Many scientists have devoted many hours to obtain an efficient algorithm for it 
and did not succeed, but the jury is still out on whether one exists. Another problem 
thought hard classically is factoring. In 1994, Shor found an algorithm which runs in time 
polynomial, rather than exponential, in the number of digits of the number to be factored. 
However, Shor's algorithm requires a quantum computer, the target of considerable effort by 
many research groups, with nuclear spins, optical lattices, quantum dots, superconducting 
circuits or trapped ions as their building blocks. In fact, a small, single-purpose, NMR-based 
quantum computer was built by Isaac Chuang's group, ran Shor's algorithm, and factored 
15. Needing to go much farther than this proof-of-principle quantum computer, today's 
effort focuses on developing and building scalable and fault-tolerant quantum computers. 

The world is quantum-mechanical. It is interesting to investigate what the laws of 
physics would allow us to compute, if we could control a quantum system and use it as 
a computer. We would like to know the power and the limits of computers that would 
use quantum instead of classical information. From understanding the details of chemical 
reactions to finding ground state properties of quantum spin glass models, there are many 
problems whose nature is quantum mechanical. It is possible that a quantum computer 
could help solve these. Not only that, there are many classical problems where algorithmic 
speedups for quantum computers over their classical counterparts have been found, and 
there is hope that many are yet to be discovered. 

Quantum mechanical systems behave according to the Schrodinger equation, governed 
by the Hamiltonian operator of the system, which describes its interactions. In our experi- 
ence, the interactions between particles involve only a few particles at a time. A Hamiltonian 
describing such a system is called local. More specifically, if the particles are also spatially 
close, I will call the Hamiltonian spatially local, and emphasize this distinction throughout 
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the text. In this thesis, I ask — and answer — questions about the power of local Hamiltoni- 
ans, whether how to simulate them, what could we do with them if we had them, or how 
hard it is to find their properties in the first place. 

0.1 Outline and Summary of Results 

After reviewing classical and quantum computing and introducing local Hamiltonians, il- 
lustrated by Feynman's Hamiltonian Computer and Kitaev's Local Hamiltonian problem 
in Chapter [H I present the following new results in the next four chapters: 

Chapter [2} Adiabatic Quantum Computation investigates the AQC model of quan- 
tum computation, which is based on time-evolving a system with a time-dependent Hamil- 
tonian from a specific initial ground state. I define a general Hamiltonian Computer model, 
show that it is universal for quantum computing, and show that it also implies the univer- 
sality of AQC for quantum computation. I improve the estimates on the runtime of AQC 
simulations of the quantum circuit model using the results I obtain later in Chapter [H 
Then, based on [3l], I prove two theorems about what not to do when designing quantum 
adiabatic algorithms, in essence saying that the Hamiltonians one uses must be local, if one 
hopes to succeed. 

Chapter [3t Matrix Product States on Infinite Trees is based on ^69j and contains 
a new numerical method for investigating the ground state properties of translationally in- 
variant Hamiltonians on an infinite tree (Bethe lattice). This method is based on imaginary 
time evolution within the Matrix Product State ansatz. I use it to analyze two specific 
interaction models, first of which is the Ising model in transverse field on the Bethe lattice. 

Chapter |4l Quantum Satisfiability analyzes the complexity of the quantum analogue 
of classical Satisfiability. I give several novel clock constructions useful for proofs of QMA 
completeness of Hamiltonian problems and for proofs of universality of Hamiltonian com- 
puters. Based on [70], I show that the 3-local Hamiltonian problem stays QMA-complete 
even without high-norm penalty terms required in previous constructions and that Quan- 
tum (3,2,2)-SAT is QMAi complete. I then present two previously unpublished results. 
First, Quantum 2-SAT on a line with 11-dimensional particles is QMAi complete. Second, 
the Quantum 3-SAT Hamiltonian is universal for quantum computing, when used in my 
Hamiltonian Computer model. 

Chapter [5l Hamiltonian Quantum Cellular Automata is based on [7l]. I define 
HQCA, a model of quantum computation with a local, time-independent, translationaly 
invariant and problem-independent Hamiltonian. I present two such automata in ID ge- 
ometry, one built from particles with dimension d = 10, and the other from particles with 
dimension d = 20. These results are then related back to the universality of adiabatic 
quantum computation. 

Finally, I summarize the thesis in Chapter [6l and discuss further research directions. Ad- 
ditional material can be found at the end of the thesis. Appendix lAl investigates Kitaev's 
propagation Hamiltonian, Appendix [B] works out the details for a Trotter-like formula for 
classically simulating an adiabatic algorithm. Appendix [C] is an analysis of the continuous 
time quantum walk on a line and cycle. Appendix [P] concerns the diffusion of free fermions 
in ID, and Appendix lEl looks at transitions in a two-level system. 
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Chapter 1 



Preliminaries 



This preliminary Chapter is aimed at a reader with httle background in classical computing 
and complexity, and only basic knowledge of quantum mechanics. Kitaev's Classical and 
Quantum Computation [56] and Nielsen & Chuang's Quantum Information and Computa- 
tion [73] are great starting points to investigate the topics mentioned herein in much more 
detail. The advanced reader should start straight with Section 11.31 on Local Hamiltoni- 
ans, where I review the basic building blocks for many results in this thesis: Feynman's 
Hamiltonian computer and Kitaev's Local Hamiltonian problem. 

This Chapter is organized as follows. I start with the classical model of computation in 
Section 11.1. 11 review classical complexity in Section 11.1.21 and then look at classical locally 
constrained problems, A:- Satisfiability and Max-/c-Satisfiability, in Section [1.1.31 

In Section [L2l I introduce quantum computation, discussing the quantum circuit model 
and universality in Section II. 2. H quantum complexity in Section 11.2.21 the strengths 
and limits of quantum computers in ll.2.3[ 

Finally, in Section [1.31 1 turn my attention to the main topic of this thesis, local Hamil- 
tonians in quantum computing, introducing Feynman's Hamiltonian computer model in 
Section 11.3.11 and Kitaev's QMA-complete Local Hamiltonian problem in Section 11.3.21 

1.1 Classical Computation 
1.1.1 Introduction to Computation 

I have a computer that I am typing this on. What can I do (compute) with it? I can 
surely add and multiply numbers really fast, but could I also use it to factor large numbers? 
It can effortlessly beat me in chess, but I am up to the challenge in gc0. Could I have 
programmed it to write this thesis? How big would the program have to be? How long 
would it take to execute? The question of computability and effectiveness of computation 
needs to be looked at in a particular model of computing. In 1937, Alan Turing introduced 
his remarkable Turing machine (TM) model [99 \ 1100] . Consider a semi-infinite tape with 
symbols Si G {0, 1, ^} on it (see Figure [Ld]) . a read/write head with a finite memory (internal 
state space) positioned on the left end of the tape, and a finite set of instructions. The head 
is initialized in the state qo and positioned on the left of the tape, reading the symbol si. 
The computation is then performed in sequential steps according to the instruction set. 

^though not against the very recent Monte Carlo programs like MoGo 
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Figure 1-1: The Turing Machine. 



These instructions have the form 

[current state], [symbol being read] — > [next state], [symbol to write], [where to move]. 

Depending on the current state of the machine and the symbol that is under the read/write 
head, the machine can change its state, write a symbol to the tape and move one step to 
the left (denoted by <), one step to the right (l>) or stay at its current position ( • ). It is 
straightforward to run the machine according to the instruction set. The execution stops 
only when the internal state of the head changes to a special state halt. The output of the 
TM is then what is written on the tape (or just the very last symbol written, right under 
the head). 

Here is a simple example of a Turing Machine, with the head starting in the state qo on 
the left of the tape: 



^0, 


- 


qo , 


0, 


>, 


(1.1) 


qi, 


- 


qi , 


0, 


>, 


(1.2) 


qo, 


1 - 


qi , 


1, 


>, 


(1.3) 


qi, 


1 - 


qo , 


1, 


>, 


(1.4) 


qo, 




HALT, 


0, 




(1.5) 


qi. 




HALT, 


1, 




(1.6) 



The first two rules (|l.ip - (|1.2p say that if the head reads a zero, it should keep its state and 
move to the right. Rules (jl.3p - (|1.4p apply when the head reads a 1. The internal state of 
the machine then flips (go ^ qi), and the head again moves to the right. Finally, when the 
head reads ^ , it assumes it has reached the end of the string, writes down or 1 according 
to its internal state, and halts. This TM computes the parity of the string such as 

10010100110 ^ ^ ^ ^ ^ ^ . . . (1.7) 

n 

initially written on the tape. The input string size n is the number of symbols on the tape 
after which the rest of the tape is blank. Observe that this TM takes at most n steps before 
it halts on any input string of size n. 

When a TM halts, it produces an output. In the above example, the output was the 
parity of the input string. However, a general TM does not have to halt on each possible 
input, leaving its output undefined. A function on a set S whose domain is not the whole 
S is called a partial function. A TM is thus computing some partial function (pM : S ^ S 
on S, the set of strings over the alphabet {0, 1, ^}. Taking it from the other end, does a 
TM computing a partial function F exist for any F : 5 — > 5? There is an uncountable 
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Figure 1-2: A binary addition circuit for two two-bit numbers, Z2ZiZq = xiXq © UiUq- 



number of such functions, while the set of possible Turing Machines is countable. This 
means not all functions can be computable, motivating the following definition: A partial 
function F : S ^ S \s computable, if there exists a Turing machine M such that (pM = F- 
An example of an uncomputable function is the halting problem: Does a given TM halt on 
all inputs? [91] A hypothetical halting problem-solving TM would then have to take the 
description of the TM whose halting properties it wants to determine as input. 

On the other hand, for functions that are computable, it is interesting to ask how hard 
it is to compute them. For this, let me first ask how long does a specific Turing Machine 
run before it stops. This question must be related to the input string size n. The running 
time T{n) of a TM is then defined as the largest number of steps it takes for the TM to stop 
on any input string of size n. There can be several different Turing Machines computing 
the same function, corresponding to different approaches to computing F. Some of them 
can be more efficient than others, with a better running time T{n). An efficient TM is 
one whose running time T(n) is small. The most common notion of an efficient TM is one 
whose running time scales only like a polynomial in the input size n (i.e. like T{n) oc 
and not like T(n) oc 2"), and I adapt it in this thesis. On the other hand, functions F can 
be classified into complexity classes depending on what would the running time be for the 
best possible Turing Machine. I return to this topic in Section [1.1. 2i 

The notion of an algorithm, a prescription of what to do to compute something, is 
dependent on the computational model one uses. Turing himself defined what an algorithm 
is by his Turing thesis: 

"Any algorithm can be realized by a Turing machine." 

Suprisingly, the alternative models of computation proposed in the early days of computer 
science turned out to be equivalent to the TM model [9ll [56]. Today's computers do not 
have long tapes with read/write heads inside them. The conventional circuit model of digital 
computers (see Figure [1T2]) . equivalent to the TM model, consists of wires carrying bits of 
information and logical gates like 



X 


NOT X 





1 


1 






(1.8) 
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and two-bit gates 
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(1.9) 



that are applied to them. What is the smahest set of gates one needs to simulate a Turing 
Machine (compute a computable function F)l It turns out that the NAND gate is universal 
in the sense that any circuit can be transformed into another built only from wires and 
NAND gates. 

The single-bit gate NOT is reversible, as its input can be recovered from its output, while 
the gates AND, OR, XOR and NAND are all non-reversible. It turns out that reversible 
circuits are as powerful as their irreversible counterparts. The reversible 3-bit Toffoli gate, 
which flips the third bit if the first two bits are both 1, is universal |97j . This is important, 
because quantum computation is a generalization of classical reversible computation (see 
Section OH). 

In the 1930's, analog computers like the differential analyzer for solving differential 
equations had their programs built in, as their components had to be physically rearranged 
for solving different problems. Similarly, for the TM mentioned so far, a different TM was 
required for computing each partial function F. However, Turing made a big breakthrough, 
paving the way for digital computation, by proving that it is possible to have a single 
universal TM, which can produce the output of any other TM. This is done by having its 
'program' (the description of which particular computation to perform) given as part of 
the input. The universal TM then processes the 'data' part of the tape according to the 
'program'. It is interesting to look for a simplest possible universal TM, and people have 
dedicated quite some effort to it [58l I112j . In Chapter [H I will return to the idea of a 
classical universal reversible TM, and bring it to the world of quantum computing. I then 
investigate simple universal translationally and time-independent Hamiltonians which could 
be used to simulate any quantum computation in a ID system. 

1.1.2 Classical Complexity 

Classical complexity examines how hard it is to compute a function F on strings S, or how 
hard it is to answer yes/no questions (decision problems) about strings S. A simple yes/no 
question about a string is: "Is the input string S a palindrom^" 

For comparing the hardness of problems, it is important to have a notion of reducibility. 
A reduction is a way of converting one problem to another in such a way that the solution 
to the second problem can be used to solve the first problem [91] . That multiplication can 
be reduced to addition, using addition as a subroutine polynomially many times, is a simple 
example of a polynomial reduction. I will frequently use reductions in this thesis, building 
on previous results about hardness of certain problems. 

The complexity of computing a function or a decision problem can be classified by 

• reducing the problem into another that is already classified, 

palindrome is a string that reads the same backwards as it does normaUy, such as "hvedondecaf- 
facednodevil" . 
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• showing that an already classified problem reduces to the one at hand, 

• showing a lower bound l{n) < T{n), 

• finding an algorithm for it that runs in some T[n). 

To prove a lower bound l{n) on the complexity of computing a function F, one must 
show that no algorithm computing F can possibly have running time less than l{n). On the 
other hand, finding efficient algorithms is hard. We can solve many problems efficiently, but 
there are plenty others for which today's best known algorithm scales exponentially with 
the problem size (see Section [TTL3]) . However, we do not know what tomorrow's algorithm 
designers will come up with. The best situation comes when a problem has a lower bound 
on its complexity, and an optimal algorithm that matches it. A trivial example is the 
problem of computing the parity of a string. The Turing Machine described by (jl.ip is 
in fact optimal, matching the lower bound of n steps for a string of length n. A more 
interesting example is the Ordered search problem, where one needs to determine whether 
a specific number appears in an ordered list of n numbers. The list of numbers can be 
written on the one-dimensional tape of a Turing Machine which in the worst case would 
have to go through it all. On the other hand, I can view this problem differently. Consider 
the list of numbers to be hidden in a black box which the TM can query by writing a 
number m on a separate 'query' tape. Upon a query, the black box subroutine writes the 
m-th number from the list on the tape. The TM can then compare this number to the 
one it searches for, do some computation, and run another query. The complexity of the 
algorithm is then determined by the number of queries to the black box the TM needs to 
ask before it determines whether a number is in the black box list or not. The complexity of 
such black-box problems is called query complexity. The best classical algorithm for ordered 
search matches the corresponding lower bound, requiring to query log2 n numbers in the 
list. In Section 11.2.31 I look at Unstructured search, another problem for which we have 
both a lower bound on the query complexity and an optimal algorithm, in the classical and 
in the quantum case as well. 

I now introduce a few important classical complexity classes, all related to decision 
problems, which are simply yes/no questions. 

Definition 1 (Complexity class P). The class of decision problems solvable in polynomial 
time by a Turing machine. 

This is the basic class of efficiently solvable problems. It contains easy problems like 
determining whether two numbers are coprime, but also the problem of determining whether 
a number is prime [3j. On the other hand, it is not known whether it contains factoring, 
the problem of determining whether a number p has a divisor d in some range a < d < h. 

71 

Today's best known classical algorithm for factoring has runtime 2 3", where n is the number 
of bits of p. On the other hand, in 1994, Shor p)OJ discovered a breakthrough quantum 
algorithm with running time T{n) oc n^, whose use of the Fourier transform resulted in 
many of today's quantum algorithms. I refer the reader to Section 11.2.31 for more details 
about efficient quantum algorithms. 

All the problems in P are contained in another class, NP, described by Scott Aaronson 
in the Complexity ZOO [1] as 'the class of of dashed hopes and idle dreams'. Determining 
whether NP is in fact bigger than P or equal to it is a long standing open question in 
Computer Science with a million dollar prize from the Clay foundation awaiting the solver. 
The definition of NP is 
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Definition 2 (Complexity class NP (Non-deterministic, Polynomial time)). NP is the class 
of decision problems such that, if the answer is 'yes', then there is a proof of this fact, of 
length polynomial in the size of the input, that can be verified by a deterministic polynomial- 
time algorithm. On the other hand, if the answer is 'no ', then the algorithm must declare 
invalid any purported proof that the answer is 'yes'. 

It is possible to expand the original Turing Machine model by allowing the TM to use 
an additional random generator, making it a probabilistic Turing Machine. This can be 
done by giving the TM access to a second tape with a random string written on it. The 
question whether this model is equivalent to the original, deterministic TM model is open, 
but believed to be true. The probabilistic Turing Machine model gives rise to the class 
BPP: 

Definition 3 (Complexity class BPP (Bounded-error, Probabilistic, Polynomial time)). 
The class of decision problems solvable in polynomial time by a randomized Turing machine. 

Great progress has been made in randomized algorithms recently, and it is a challenge for 
researchers to derandomize them. I give an example of a randomized algorithm in Section 
11.1.31 (the WALKS AT algorithm for solving 3-SAT). Moreover, quantum computing, the 
topic of this thesis, is inherently randomized. 

Analogously to generalizing P to BPP, the class NP motivates the class MA. Consider an 
interactive protocol between a powerful wizard Merlin and a logic, rational verifier Arthur, 
whose role is to verify purported proofs prepared by Merlin. The difference from NP is that 
Arthur can use a randomized algorithm, rather than verify the proof deterministicalljd: 

Definition 4 (Complexity class MA (Merlin-Arthur)). The class of decision problems solv- 
able by a Merlin-Arthur protocol, which goes as follows. Merlin, who has unbounded com- 
putational resources, sends Arthur a polynomial-size purported proof that the answer to the 
problem is yes. Arthur must verify the proof in BPP (i.e. probabilistic polynomial-time), so 
that if the answer to the decision problem is 

1. 'yes', then there exists a proof which Arthur accepts with probability at least 2/3. 

2. 'no', then Arthur accepts any "proof" with probability at most 1/3. 

An alternative definition requires that if the answer is 'y^s', then there exists a proof 
such that Arthur accepts with certainty. However, the definitions with one-sided and two- 
sided error are equivalent [91]. Moreover, the probabilities 2/3 and 1/3 in the definition of 
MA can be equivalently substituted by 1 — e and e with arbitrarily small e. 

I introduce the quantum analogues of these classes in Section 11.2.21 

1.1.3 Classical locally constrained problems 

Whether deciding on the best class schedule, finding a good route to spend a day touring a 
city or pairing people up at a dance class, our lives are full of local optimization problems. 
Some of them are simple, such as picking a single box of cereal with the lowest price/taste 
ratio from the shelf in the supermarket, some of them harder, like deciding which lenses 
to pack into my camera bag for the next trip. The latter is a question concerning a cost 
function involving only a few elements at a time (how useful will a particular lens in my bag 

^If all random algorithms can be made deterministic (BPP=P), MA equals NP [1]. 



24 



X 




y 




z 




w 










X 




... 






x' 














> 







Figure 1-3: The implication graph for a 2-SAT instance (x V y') A {x' V y') A (y V tf) A (y V 
z') A (z V w'). There is a directed loop containing both y and y' , so the instance is not 
satisfiable. 

be, considering that I already bagged some others) and a global constraint (the size of the 
bag and my willingness to carry it all). In fact, this is an example of the Knapsack problem. 
Another locally constrained problem. Satisfiability, and its quantum analogue, are the focus 
of a large part of this thesis. 

A boolean clause is an expression of the form si V S2 V • • • V s^, where Si are boolean 
variables or their negations. This clause is true for all bit strings siS2 ■ ■ ■ Sk except for 
{00 ... 0}. Consequently, the clause si V Sg, where is the negation of bit S2, is true for all 
assignments siS2 except for {01}. Satisfiability is then defined as 

Definition 5 (A;- SATISFIABILITY, A;-SAT). Consider a collection C = {ci,C2, . . . ,Cm} of 
boolean clauses for a finite binary string S = S1S2 ■ ■ ■ Sn, such that each Ci involves k of the 
string variables. Does a string S satisfying all the clauses in C exist? 

This is an instance of 2-SAT on four bits: 

Cl = X V y', 

c-2 = X Vy , 

C3 = y Vw, (1.10) 

04 = y y z\ 

c^ = z y w' . 

Here y' labels the negation of variable y, and the first clause, ci, is false if 2; = and y = 1. 
The complexity of A:— SAT varies with k. For 2-SAT, we have an algorithm with complexity 
scaling as O(n^). Draw a directed graph of implications between the bits and their negations 
as in Figure ll-3t where for example clause ca implies the connections z =^ yandy' =^ z' . 
If a loop containing both a variable and its negation exists (it is easy to check this), the 
instance of 2-SAT is unsatisfiable. 2-SAT is thus in P. However, no known polynomial 
algorithm for 3-SAT exists. The Cook-Levin theorem says that 3-SAT is NP-complete. A 
problem H is complete for a complexity class, when 

• the problem H is contained in the class, and 

• every other problem in the class can be reduced to H (solved by using the solution to 
n subroutine polynomially many times). 

It is straightforward to see that 3-SAT is in NP. The basic idea for the second direction of 
the NP-completeness proof is to encode the verification circuit in the definition of NP using 
a 3-SAT instance. It is instructive to see how this type of proof works in the classical case. 
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Figure 1-4: Illustration for the Cook-Levin theorem. Transforming a verifier circuit into a 
3-SAT instance. 



before I introduce its quantum analogue in Section I1.3.2I An instance rr of a problem 11 in 
NP is a yes/no question, denoted as x G {Lyes'J Lno), where Lyes is the set of instances with 
the answer 'yes'. Each problem instance has a verifier circuit V associated with it. When 
V checks a purported proof that x £ Lyes, it accepts only a valid proof, never outputting 
1 if the "proof" is wrong. This verifying procedure can be reduced to finding a satisfying 
assignment for a certain 3-SAT instance. Without loss of generality, circuit V consists 
of wires and gates involving at most 3 bits (two inputs, one output). Assign a boolean 
variable to each wire of the circuit as in Figure 11-41 The 3-SAT clauses then check the 
proper evaluation of the circuit. For example, the boolean clauses for y,z,w in Figure [T^ 
(where w is the output of the AND gate on y and z) read 

Cl = y \J z \Jw', (1.11) 

C2 = y \/ z' \/ w', (1.12) 

eg = y' V z V w', (1.13) 

C4 = y' V / V w, (1.14) 

ruling out the assignments yzw G {001,011, 101, 110}. The last ingredient in the proof is a 
clause verifying that the circuit outputs 1: 

Cout = w. (1.15) 

If a proof of X G Lyes which the circuit V accepts exists, an assignments of the bits cor- 
responding to the computation satisfying the 3-SAT instance exists as well. On the other 
hand, if there is no "proof" that V accepts, there is no satisfying assignment to all the 
clauses in the 3-SAT instance. Therefore, if one could solve 3-SAT, one could solve any 
problem in NP. This is also true for /c-SAT with k > 3. 

To illustrate the idea of randomized algorithms mentioned in Section [LL2l let me present 
two randomized algorithms for the SAT problem. The procedures are given in Figure 11-51 
The greedy algorithm GSAT [86] starts with a random bit string and proceeds to flip bits 
which give the greatest decrease in the number of broken clauses. It works great for 2-SAT, 
but not for 3-SAT, as it gets stuck in local minima. The second algorithm, WALKSAT 
[85], goes around this problem by randomly choosing between performing a step of GSAT 
or flipping a random bit in some unsatisfied clause. This randomized algorithm performs 
surprisingly well on many real-life instances of 3-SAT. 
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GSAT 



WALKSAT 



for i = 1 to maxjtries 

S = & randomly generated string 

for j = 1 to max flips 

if S satisfies C, return S 
else flip any variable in S that 

results in the greatest decrease in 
the number of unsatisfied clauses 

end 

end 

return "No satisfying assignment found." 



With probability p: 

follow the standard GSAT 
scheme, making the best 
possible greedy move. 

With probability 1 - p: 

instead of making a GSAT step, 
randomly pick a variable 
occuring in some unsatisfied 
clause and flip it. 



Figure 1-5: The GSAT and WALKSAT algorithms. 

^ IB, IC /\ Al, CI f 




Figure 1-6: A reduction of Graph 3-colorability to 2- Satisfiability for trit-bit pairs, also 
called (3,2)-SAT. The trits with possible values {A,B,C} are denoted by triangles, while 
the bits are denoted by circles. Each clause is represented by the pair of symbols which it 
rules out. For Graph 3-colorability, the excluded assignments are AA, BB and CC, ruling 
out the same color for two neighboring vertices. 

SAT for non-bit strings 

SAT can be generalized to non-bit strings. The simplest variant of this is (3,2)-SAT, with 
each clause involving a trit (a three- valued letter) and a bit. Although it doesn't seem much 
more complicated than regular 2-SAT, already this variant is NP-complete. It belongs to 
NP, because it can be reduced to 3-SAT by encoding the trit into two bits. On the other 
hand, it is also NP-hard, as the NP-complete Graph 3-colorability problem can be reduced 
to (3,2)-SAT as depicted in Figure 11-61 

MAX-fc-SAT 

The last variant of Satisfiability I introduce is MAX-A;-SAT. In regular /c-SAT, the question 
is whether an assignment satisfying all clauses exists or not. Sometimes, showing that no 
such assignment exists can be easy. However, the MAX-/c-SAT problem poses a harder 
question: does an assignment breaking fewer than a clauses exist? The clauses in MAX- 
A;-SAT define a cost function, and solving MAX-/c-SAT for different a would determine 
its minimum value. Although regular 2-SAT is in P, MAX- 2-SAT is NP-complete, so the 
existence of a polynomial algorithm for it is highly unlikely. 
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Classical 


bit strings 


general strings 


A:-SAT 


k = 2 
k>3 


in P 

NP-complete 


(3, 2)-SAT : NP-complete 


MAX-fc-SAT 


k>2 


NP-complete 





Table 1.1: Known complexity for classical satisfiability problems. (3, 2)-SAT in the general 
strings column is the problem where each clause involves one trit (letter with three values) 
and one bit. 



In Table [LT| I summarize the known results about classical satisfiability problems. The 
goal of Chapter H] is to classify the complexity of the quantum analogues of these. I refer 
the reader interested in NP-complete problems and other variants of SAT to the book of 
Carey and Johnson |42j . 



1.2 Quantum Computing 

Up until now, I have considered the model of computing governed by the laws of classical 
physics. One needs to ask what is the computing power of systems governed by quantum 
mechanics. Alternatively, how hard is it to simulate quantum-mechanical systems with 
classical computers? Already Feynman stated these questions in [38J. I assume the reader 
is familiar with quantum mechanics, so the following exposition will be brief, introducing 
only those topics which I later use. 

Classical data can be easily measured, copied and erased. If someone gave me a binary 
string like 1001001, I could easily read it (measure it) and make many copies for myself. 
However, the world is quantum-mechanical. The state space of the simplest quantum system 
- a spin-i particle, a qubit, is much bigger than the two possible values {0, 1} of a classical 
bit. The first postulate of quantum mechanics is 

• Associated to any isolated physical system is a complex vector space with inner prod- 
uct (i.e. a Hilbert space) known as the state space of the system. The system is 
completely described by its state vector, which is a unit vector in the system's state 
space. 

The state space of a qubit is a 2-dimensional complex vector space 7i = C'^. The state of a 
qubit is given by 

= co|0) + ci|l), (1.16) 

where |0) and |1) are basis vectors in the usual Dirac notation and co,ci are complex 
numbers for which |coP + |cip = 1. 

How can a quantum computer be built from a system of qubits? In a classical computer, 
one can apply irreversible gates like NAND or reversible operations like the Toffoli gate. 
Besides measurements which I discuss below, the available operations for a quantum sys- 
tem are unitary transformations. Quantum mechanical systems obey the time-dependent 
Schrodinger equation. The next postulate of quantum mechanics is 
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• The time evolution of the state of a quantum system is described by the time- 
dependent Schrddinger equation 

ih^^im) = H{t)\m), (1-17) 

where H{t) is a Hermitian operator cahed the Hamiltonian. Note that the value of 
the constant h can be absorbed into the definition of H. 

Because the Hamiltonian is Hermitian, the time evolved state \ip{t)) and the initial state 
\ip{0)) are related by a unitary transformation Utfi'. 

\m) = ut,o\m)- (1-18) 

When the Hamiltonian H is time-independent, this unitary transformation is 

Ut,o = e-*^*. (1.19) 

Unitarity of U means that uJq = U^q, and implies that time evolution according to the 
Schrodinger equation is reversible. The transformations one can perform on a quantum 
computer are thus unitary. 

The final ingredient for quantum computation are measurements, used to read out the 
result. They could be used throughout the computation, but in the next Section I argue 
why it is sufficient to use measurements only at its end. Quantum measurements work in 
the following fashion: 

• Consider a quantum observable M, a Hermitian operator acting on the state space of 
the system being measured, expressed in terms of its eigenvectors as 

M = Y,XMf){'Pf\- (1.20) 

i 

If the state of the quantum system is {ip) immediately before the measurement then 
the probability to obtain the result A, is given by 

p{i) = \{cPf'\^)\\ (1.21) 
and the state of the system after the measurement is the corresponding eigenstate 




Imagine now that I want to determine what state the qubit ()1.16p is in. As an example, 
let me choose to measure = |0) (0| — |1) (1|, the z-component of the spin. I will get the 
result +1 with probability |coP, and the result —1 with probability |cip. It takes many 
copies of 1^) and many measurements to determine cq and ci in general. Even if I made 
different measurements, distinguishing nonorthogonal states perfectly is impossible |73j . 

1.2.1 The Quantum Circuit model 

Quantum circuits are a direct generalization of classical reversible circuits, using unitary 
transformations instead of classical reversible gates. A few basic single-qubit unitary gates, 
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Figure 1-7: The components of a quantum circuit. 



expressed as 2 x 2 matrices acting on column vectors in are the Pauli matrices 
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and the Hadamard gate H, the phase gate S and the f gate T: 
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A general single qubit unitary gate can be expressed using the Pauli matrices as 

U 



(1.23) 



(1.24) 



where i;^ is a real number h = {nx,ny,nz) is a unit vector and h ■ a = UxOx + f^y^^y + f^x<^z 
At least one two-qubit gate is required for quantum computation. The following two 
gates apply or to the target qubit, if the control qubit is in the state |1): 



CNOT 



10 
10 
1 
10 



C-Z 



10 
10 
10 
-1 



(1.25) 



The elements of a quantum circuit are usually depicted as in Figure 11-71 At the end of 
the computation, the output of a quantum circuit is read out by performing a few final 
measurements. Although I could perform measurements throughout the computation and 
apply gates conditioned on the outcomes, it is possible to defer measurements until the very 
end of the computation by using controlled gates, as shown in Figure [T^ 



Universality 

The NAND gate is universal for classical circuits and the Toffoli gate is universal for classical 
reversible circuits as mentioned in Section 11.1. 11 On the other hand, can every unitary 
transformation on n qubits be decomposed into many applications of a small finite set of 
gates? Indeed, universal gate sets for quantum computing exist in the sense that they can 
approximate any unitary operation arbitrarily well. A universal gate set must be dense in 
the group SU{n), and the approximations must converge fast. I point the interested reader 
to Appendix 3 of [73\ for a thorough discussion of the Solovay-Kitaev theorem concerning 
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Figure 1-8: A quantum circuit with measurements performed during the computation and 
its version with measurements deferred until the end of the computation. 

the convergence rate. Some well known gate sets universal for quantum computing are: 

• single qubit gates, CNOT [73], 

• single qubit gates, C-Z [73] . 

• H, S, T, CNOT [73]. 

Another universal gate set uses the quantum Toffoli gate: 

• H, Toff oh [88] 

It is intriguing, because it uses only real matrices. It is possible to do quantum mechanics 
with real numbers, as superpositions and minus signs, not imaginary numbers are essential. 
Consider the state (I1.16P of a qubit, and add to the system another qubit labeled r. This 
qubit can be used to work around the requirement of imaginary amplitudes. 

1^) = Re(co)|0)|0), + Im(co)|0)|l), (1.26) 
+ Re(ci)|l) |0), + lm(ci)|l)|l), (1.27) 

is a state of the larger system with only real amplitudes, and has all the information con- 
tained in the original state {ip) (I1.16p . Arbitrary unitary transformations on can be 
then mapped onto real unitary transformations on the larger system [88 1. Finally, a single 
two-qubit gate 

• W, 

the controlled ^-rotation about the y-axis is universal for quantum computation [88]. In 
Chapter |5l I use this single universal gate in my Hamiltonian Quantum Cellular Automaton 
constructions. 

1.2.2 Quantum Complexity 

In recent years, quantum complexity classes have been defined and studied in an attempt 
to understand the capacity and limitations of quantum computers and quantum algorithms 
and their relation to classical complexity classes. The basic quantum complexity class is 
BQP, the analogue of BPP. 

Definition 6 (Complexity class BQP (Bounded error. Quantum, Polynomial time). Con- 
sider a decision problem 11 = {Lyes U Lno); where Ly^s is the set of all instances x of the 
problem 11 with the answer 'yes '. The problem 11 is in BQP, if there exists a uniform family 
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of circuit^^ U, with the following properties. For a problem instance x with the answer 'yes' 
(x G Lyes), the probability of the circuit U outputting 'yes' is greater than |. On the other 
hand, when x G Lno, the probability of the circuit outputting 'yes' is less than |. 

In short, BQP is the class of problems one can solve efficiently using a quantum com- 
puter. The complexity class QMA, also known as BQNP, studied and defined in ^57j and 
[56] . is the quantum analogue of the classical complexity class MA, the probabilistic setting 
of NP. 

Definition 7 (Complexity class QMA). A decision problem U = {Lyes U Lno) of size n is 
in the class QMA if there exists a polynomial time quantum verifier circuit V such that for 
every instance x of the problem 

1. \/x £ Lyes '■ there exists a witness state \'^{x)) such that the computation V \x) ® 

|0 . . . 0)anciUa Ijidds the answcr 1 with probability at least p ; 

2. Va; G Lno ■ for any witness state \(p{x)), the computation V \x) 10 . . . 0)anciiia 
yields the answer 1 with probability at most p — e, 

where p > and e = / poly{n)) . 

The first known QMA-complete problem, local Hamiltonian, is a quantum analogue of 
classical MAX-k-SAT. Building on ideas that go back to Feynman ^38j, Kitaev [56] has 
shown that 5-local Hamiltonian is QMA complete and I present this result in Section I1.3.2I 
A special case of Local Hamiltonian, Quantum A;-SAT, is the topic of Chapter 31 There 
I show that several of its variants are complete for QMAi, the class QMA with one-sided 
error: 

Definition 8 (Complexity class QMAi). The class QMAi is the class QMA with single 
sided error, i.e. with p = 1 in the above definition. 

In the classical case, MA with one-sided error is equivalent to regular MA. It is not 
known whether QMA=QMAi. 

1.2.3 The Strengths and Limits of Quantum Computing 

There are quite a few problems on which quantum computers would algorithmically out- 
perform classical ones. In quantum computing, we can use resources unknown in classical 
computing: superpositions, interference and entanglement. Much has been discovered since 
Shor's factoring algorithm in 1994 [90]. Stephen Jordan has a thorough and up-to-date 
summary of the algorithmic successes of quantum computing in his thesis [51]. However, 
the speedups are rarely exponential. Not only that, solving problems such as Parit'^ cannot 
be sped up on a quantum computer at all [35]. I refer the reader interested in the limitations 
of quantum computers to [HI [21 [TSl [17] . 

Algorithm designers bound the complexity of a problem from above by finding better 
algorithms, while lower bound enthusiasts limit it from below. The quantum variant of 
the black-box problem Unstructured search is an example of a problem for which these 
two approaches meet. The lower bound on its query complexity is known, and an optimal 
quantum algorithm, Grover's search, exists. 

*Here uniform family of circuits means that there is one circuit for each problem size (and not for each 
problem instance), with a number of gates polynomial in the problem size. 
'^Parity: determine the sum modulo 2 of n bits. 
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Definition 9 (Unstructured Search). Consider a quantum black-box (oracle) 



0\s) = I-2\w){w 



(1.28) 



which flips the phase of a single, unknown to you, canonical basis quantum state \w) , and 
acts trivially on any other canonical basis state. Find the state \w) . 

Looking at a single state of n bits unsuccessfully does not give me any information 
about where to look next, making this problem unstructured. Classically, the required 
number of queries scales like the number of possible n-bit strings, 2", and is matched by the 
simplest randomized algorithm one can think of — check the strings randomly until you find 
the marked one. The query complexity of a quantum algorithm for unstructured search is 
the number of necessary calls to O in the worst case. In 1996, Grover discovered his famous 
algorithm requiring only \/2" queries, getting a square root speedup over the best possible 
classical algorithm. Also around that time, Bennett, Bernstein, Brassard and Vazirani found 
a matching lower bound l{n) oc \/2" [17]. Therefore, Grover's algorithm is the best possible 
quantum algorithm for unstructured search. Unstructured search can be also recast in a 
continuous setting, where quantum computation is done by continuous-time Hamiltonian 
evolution instead of application of discrete unitary gates. Instead of a quantum black box, 
in [37j Farhi and Gutmann use an oracle Hamiltonian 



where E determines an energy scale. This oracle Hamiltonian is always on in the system. 
What one can do to find \w) in this Hamiltonian setting is to prepare the system in some 
initial state, and let the system evolve according to a Hamiltonian 



where Hoit) is an additional arbitrary time-dependent Hamiltonian. If measuring the 
system after some time produces \ w), the algorithm is successful. Farhi and Gutmann show 
that the required running time for any such Hamiltonian algorithm (any Hoit) one uses) 
necessarily scales like jE. In Chapter [21 I use this lower bound to show what not to do 
when designing quantum adiabatic algorithms. 

1.3 Local Hamiltonians 

in Quantum Computation 

Finding the properties of relatively simple quantum systems is the topic of much condensed 
matter physics. Even though all the interactions are local, involving only a few neigh- 
boring particles, the resulting possibilities are amazingly rich. In this thesis I focus on 
two properties of systems governed by Local Hamiltonians. First, a (quantum) computer 
can be constructed from a system governed by a simple Hamiltonian. Second, the ground 
states of local Hamiltonians can be hard to find, as they can encode solutions to interesting 
optimization problems. 

Hamiltonian time evolution according to the Schrodinger equation is the underlying 
principle of all quantum computing, producing the unitary operations required the quan- 
tum circuit model. Controlling the Hamiltonian of a quantum system means changing the 




(1.29) 



H{t) = HD{t) + H, 



(1.30) 
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interaction of the particles with each other and their surroundings. For example, in addi- 
tion to the natural interaction of the particles in a quantum system, we can apply magnetic 
field to spins on a lattice, laser pulses to ions in an ion trap or electric fields to electrons 
in quantum dots. External control can allow us to perform the desired unitary transforma- 
tions in the system, making it into a computer. Feynman [38] thought of how to make a 
Hamiltonian Computer for simulating other quantum systems long before people invented 
quantum circuits. I introduce his idea in Section II. 3. H and build on it throughout this 
thesis, especially in Chapters H] and O 

Besides investigating the time evolution of a quantum system, analyzing the spectrum 
of its eigenvalues is another source of interesting problems. It is long known that finding 
ground states of certain Hamiltonians is hard. If we could find ground states of a system 
governed by the Ising model 3D, we could solve NP-complete problems [13] . Finding minima 
of locally constrained optimization problems can be recast into finding the ground state 
energy of a local Hamiltonian. I am interested in showing how hard it is to find it for 
some Hamiltonians. In one case, in Chapter [3l I give a numerical method for finding the 
ground state energy of the translationally invariant Ising model in transverse field on a 
Bethe lattice, based on the Matrix Product State ansatz. On the other hand, building on 
Kitaev's Local Hamiltonian problem I introduce in Section 11.3.21 I classify the complexity 
of finding the ground state of several Hamiltonians in Chapter [H 



1.3.1 Feynman's Hamiltonian Quantum Computer 

Already in 1985, Feynman [38] proposed the following interpretation/implementation of a 
computer using a quantum-mechanical systems. What he originally had in mind was to 
simulate a sequence of reversible classical gates in a quantum mechanical system. However, 
his construction applies equally well to implementing a quantum circuit. Take a sequence 
of L unitary transformations Ut on n work qubits 

U = UlUl-i...U2Ui. (1.31) 

Consider now a system with two registers, 

H = Hwork ® 'Hclock- (1-32) 

The first one holds the n work qubits, which I label qi, . . . ,qn throughout the thesis, while 
the second register holds a pointer particle hopping on a line 0, 1, . . . , L, serving as a clock 
for the computation. I can realize the clock register using L + 1 qubits labeled cq, . . . ,cl, 
and initialize it in the state 

|0), = |l,,0,,Oe,...0,J. (1.33) 
The state of the clock register corresponding to "time" t is 

|t),= |0e„...0,,_,le,0e,+,...0,,>. (1.34) 

A simple time-independent Hamiltonian can facilitate the evaluation of U on the work 
register of this system. As the pointer particle hops on the line 0, 1, . . . , L, performing a 
quantum walk, the corresponding unitary operations are applied (or uncomputed) on the 
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Figure 1-9: Feynman's computer with a hopping pointer particle (a single spin up). When 
the active site moves from C2 to C3, gate C/3 gets applied to work qubits gs and q^. When 
the active spin jumps back from C2 to ci, gate U2 is uncomputed on qubits (72 and 93. 



work qubits as depicted in Figure [T^ The Hamiltonian with the desired dynamics is 

L 

Hprimn = '^Ut(^ ^(t)^(t-l) + ® ^{t)^{t-l)^ (1.35) 
t=l 

where cx^^ = |1) (0|^^ and u^^^ = |0) are the raising and lowering operators for the clock 
register spin q and Ut acts on the corresponding work qubits. When the system is initialized 
in the state 



l^o) = iV'o) ® |0), = IVo) |10 . . . 0), , (1.36) 

where IV'o) is some initial state of the work qubits, time evolution with Hpnmn according 
to the Schrodinger equation brings the state into a superposition of states 

l^t) = |V't)®|i)c, (1-37) 
where iV't) is the state of the work qubits after the first t gates of the circuit: 

\il^t) = UtUt-i...Ui\iPo)- (1.38) 

After some time r, when I measure the clock register of |^') and obtain L, the state of the 
system after the measurement becomes ® l-^)c- '^^^ work register of this state contains 
the desired output of the circuit U . 

I can obtain with high probability using a slightly modified system. First, I pad 
the sequence U (jl.3ip with 2L extra identity operations 

[/'=l3Lll3L-l...IlL+lC/, (1.39) 

and then expand the clock register to 3L + 1 qubits accordingly. All the states \iljt) (|1.38p 
for t > L now are the same, equal to the desired output state |V'l)- Starting from (ll.36p . I 
let the system evolve for a time r chosen uniformly at random between and 0(L log L). 
Because the quantum walk on a line is rapidly mixing, the probability to measure a state 
with the clock register \t > L)^ at time r is close to |, as shown in AppendixO This yields 
the desired state IV'l) of the work qubits. If a measurement of the clock register results in 
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\ipt) |t < L)^, repeat the experiment. 

The clock register representation of the states \t)^ given in (jl.34p is only one of the 
possibilities. I describe a different one in the next Section, and several others in Chapter SI 
Considering this, Feynman's Hamiltonian can be rewritten in a general form, encompassing 
different implementations of \t)^, as 

T 

t=l 

H'p = Ut0Xt,t-i + {Ut®Xt,t-i)K (1.41) 

where Ut acts on the corresponding work register qubits, and the operator Xt^t-i acts on 
the clock register as 

Xt^t^i = \t){t-l\^, (1.42) 

increasing the time from t — 1 to t. In Feynman's original formulation, this corresponds to 
moving the position of the single up spin to the right as |. . . 0100 ■ ■ — > |. • • 0010 ■ ■ ■)^, 
and the operator Xt^t-i can be written as Xt^t-i = |01) (lOl^^ .^ where q_i and q are 
two of the clock register qubits. Analogously, the Hermitian conjugate of Xt^t-i moves the 
time register backwards as 

(Xi,i_i)t = \t-l){tl. (1.43) 

This is my preferred formulation of Feynman's Hamiltonian throughout this thesis. It shows 
up in Chapter [2] where I use it to prove results about the universality of adiabatic quantum 
computing, while in Chapter [J] it serves a basic building block for my results about the 
complexity of Quantum Satisfiability. Finally, in Chapter [5] it inspires my Hamiltonian 
Quantum Cellular Automaton model. 

1.3.2 Kitaev's Local Hamiltonian (LH) Problem 

In this section I introduce Kitaev's Local Hamiltonian and sketch the proof that it is QMA 
complete. In |56j Kitaev defined the following promise problem: 

Definition 10 (Local Hamiltonian). Take a k-local Hamiltonian H = Yl]=i-^j o,cting on 
n qubits, composed of r = poly{n) terms, each of those acting nontrivially on a constant 
number k of qubits. Given two numbers a, b, where b > a and the separation b — a is greater 
than n~" for some constant a, determine whether 

1. H has an eigenvalue not exceeding a, or 

2. all eigenvalues of H are greater than b. 

The two numbers b and a are separated by at least an inverse polynomial in the number 
of qubits in the system, n. LH is a promise problem, therefore we know that either 1 or 2 
are true. Kitaev proved that 

Theorem 1 (LH is QMA-complete) . 5-Local Hamiltonian is QMA-complete. 



36 



One of the directions of the proof is to show that LH is in QMA. The verifier circuit C/, 
saying whether Arthur accepts a given state or not, is constructed as follows. Add an 'an- 
swer' qubit to the system and for each term Hj = X] iV's) (^s|) construct a measurement 
operator 

Wj : \A) ® |0),„, - IV's) ® (V^m + Vl-A.ll)) . (1.44) 

\ / ans 

Taking a general state \rj) = XlsyslV's)) the probability the circuit accepts it (i.e. that I 
measure 1 on the output of C/) is 

p,il) = {v\^{OLnsWjiI<S)\l){lUJW,\v)0\O),,,, (1.45) 
= Y.(^-XMys = l-{r]\Hj\r^). (1.46) 

s 

Add now another 'operator' register to the system. The general circuit for checking all of 
the terms in the Hamiltonian combines the checking operators Wj into 

r 

W = ^\j){jlp^W„ (1.47) 
i=i 

where r is the number of terms in the Hamiltonian. After applying W to the system in the 
state 

Elj)opj ®l^)»|0>,n., (1.48) 

one measures the answer qubit. The probability of getting the outcome 1 is then pi = 
1 — M^M, This is enough to show that LH is in QMA. 

For the other direction, one needs to construct a corresponding instance of the Hamil- 
tonian problem for every instance of a problem in QMA. Each of those has a verifier circuit 
U with the required accept /decline properties associated with it. Given a quantum circuit 
U, Kitaev showed how to construct H with the properties required in the definition of LH 
(Definition [TUD . If the circuit U accepts the state IV'), this Hamiltonian's ground state can 
be constructed from it, and the corresponding expectation value in this state will be smaller 
than a. On the other hand, if the circuit does not accept any state with high probability, 
the ground state of H will have energy higher than b. 



Kitaev's Hamiltonian 

Consider a quantum system consisting of a work register with n qubits qi,. . . ,qn, and a 
clock register. 

= T~(-work ® T~(-clock- (1.49) 

Kitaev's Hamiltonian is a sum of three terms, 

H Kitaev — Hp^op ~\~ Hq^i + Hinput- (1.50) 

The first term checks the proper progression of the circuit U. Taking the idea from 
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Feynman's Hamiltonian computer (see Section ll.3.ip . Kitaev used (jl.4ip to construct a 
Hamiltonian whose ground state is the uniform superposition over the computation called 
the history state: 

1 ^ 

\^)history = ^f=fEl^*)®l*)c' (1-51) 

where 



iV^t) = UtUt^i...U2Ui\iPo), (1.52) 
and is any initial state of the work qubits. This Hamiltonian is 



L 

^ I , 

t=i 



1 ^ 

o ^prop^ (1.53) 



where each term Hp^^p is a projector constructed from Feynman's Hamiltonian ()1.4ip as 
H'p,,p = I (Pt_i + Pt) - (Ut Xt,t^i + (C/t ® Xt,t„i)t) , (1.54) 

V ' 

with Pt = \t) {t\^ the projector onto the state \t)^ of the clock register, and the time- 
increasing operator Xt^t~i changing the clock from t — 1 to t given by (jl.42p . The only 
states that are not obviously annihilated by Hp^^p live in the subspace of states with the 
clock register in the states \t — 1)^ or \t)^. Within that subspace, the states with zero energy 
have the form 



\^) = \a) 0\t-l)^ + Ut\a) (^\t)^, (1.55) 

making it easy to check that the history state ()1.5ip is annihilated by the non-negative 
Hamiltonian Hprop- It can be thought of as checking the correct propagation of the compu- 
tation, as the expectation value of Hprop in any state that is not a history state is greater 
than zero. In fact, the second lowest eigenvalue can be bounded by Ai > for some 
constant ci (see Appendix [A] for a detailed analysis of the Hamiltonian ()1.53p ). 

The last of the work register qubits, qn, is the designated output qubit for the circuit 
U. The second term Hout in (|1.50p adds an energy penalty to states whose g„ at time L is 
|0), meaning that the circuit did not accept them. 

Hout = \0){0\g,,^'^\L){Ll. (1.56) 

This term ensures that H Kitaev (ll-SOp can have a low eigenvalue only if the circuit U outputs 
1 on some input state IV'o) with high probability. 

The work register has n qubits gi, . . . , g„. Out of these, at time f = 0, some are initialized 
in a purported proof state, while others are ancilla qubits required for the verification circuit. 
The last part of H Kitaev, 

H.nput= (1-57) 

k&ancilla 
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checks the correct initiahzation of these ancilla qubits by adding an energy penalty to states 
whose anciUae at time t = (when the clock register is in the state |0)^) are not in the state 
|0). 

These three terms together make up Kitaev's Hamiltonian (jl.SOp . Let me now sketch 
the highlights of the proof that (jl.50p has the properties described in Definition llOi First, 
the expectation value of Hxitaev in the history state (jl.5ip corresponding to an initial state 
{ipyes) that the circuit U accepts can be shown to be upper bounded by Aq*^** < for some 
e. Second, one needs to lower bound the lowest eigenvalue of Hxitaev in the case the circuit 
U doesn't accept any state with high probability, and show that it is separated from Ag^* at 
least as an inverse polynomial in L. All three terms in (|1.5Up are non- negative. The proof 
uses a geometrical argument, showing that the angle between the null spaces of Hprop and 
{Hinput + Hout) is lower bounded by 9, an inverse polynomial in L. The smallest eigenvalue 
of Hxitaev can then be bounded from below by the smallest nonzero eigenvalue of Hprop 
multiplied by a geometric factor coming from 6. The resulting lower bound on the lowest 
eigenvalue of (|1.50p is then \q° > c(l — ^/e)L~^''. Hxitaev thus has the desired properties. 

The final step in the proof that 5-local Hamiltonian is QMA complete is to implement 
the clock register. I describe several clock register constructions in detail in Section 14.1.31 
One of these is Kitaev's domain wall clock. Similarly to Feynman ()1.34p . he encoded his 
clock register in unary using L + 1 clock qubits q . However, the states \t)^ are now encoded 
as 

= |l,Jei...lc,0c,+,...0e,_,0e^> (1.58) 

The Hilbert space of the L + 1 clock qubits cq, ... ,cl is much bigger than the space spanned 
by the proper clock states \t)^ given by (jl.SSp . which I call the legal clock subspace, Hiegai- 
By adding a clock-checking Hamiltonian 

L-l 

Hciock = El01)^01|ct,c,+i (1-59) 
t=o 

to Hxitaev (!1.50p . One can ensure that the low-energy eigenstates of H are close to the 
subspace of states with a proper unary-encoded clock, i.e. whose clock register states are 
close to the legal clock subspace TCiegai ■ The final form of the Hamiltonian is then 

H Kitaev^ — Hprop ~\~ Hq^i -\- Hinput ~l~ H^lock- (1.60) 

In this clock implementation, Pt, the projector onto the state \t), is only a 2-local operator 

= |t)(t|, = |10)(10|,^^,^^^, (1.61) 

with a special case P£ = |1) (l|c^, acting trivially on the rest of the qubits. On the other 
hand, the operator Xt^t-i (|1.42|) in this encoding is 3-local (acting nontrivially on three 
clock qubits): 

Xt,t„i = |t)(t-l|, = |110)(100|,^_^^,^^,^^^, (1.62) 
with a special case at the right end of the clock register, 

Xl,l-i = \L) {L - 11 = 111) (10|,^_^^,^ . (1.63) 



39 



Recalling that 2-qubit gates are universal for quantum computation, the terms in Hprop 
()1.53p . and thus also Hxitaev, become at most 5-local. It remains to show that in the case 
the circuit U rarely outputs 1, the lowest eigenvalue of (jl.60p is also bounded from below 
by an inverse polynomial in n. This proof follows from the geometric arguments outlined 
above without much complication. Therefore, 5-local Hamiltonian is QMA-complete. 

In Chapter HJ I return to this problem, focusing on a special case called Quantum 
Satisfiability, defined by Bravyi in [20j. When all the terms in H are projectors (non- 
negative), one can ask whether there exists a state with energy exactly 0, or whether the 
ground state energy of H is greater than b. I investigate the complexity of several variants 
of Q-SAT and utilize the techniques used in Kitaev's proof I just outlined. 
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Chapter 2 

Adiabatic Quantum Computation 



Farhi et al. [36] introduced Adiabatic Quantum Computing (AQC) as a method to find a 
minimum of a classical cost function by using a quantum system with a time-dependent, 
local Hamiltonian. I describe this model in Section 12.11 and then review the proof of 
Aharonov et al. |7j, who showed that AQC is equivalent to the quantum circuit model, in 
the light a previously unpublished result of Seth Lloyd [63] in Section I2.2[ This motivates 
my definition of a Hamiltonian Computer model in Section 12.2.11 which I later use to prove 
universality of the Quantum 3-SAT Hamiltonian in Section 14.51 and of two Hamiltonian 
Quantum Cellular Automata in Chapter [5j 

There is hope that there may be combinatorial search problems, defined on n bits so 
that = 2", where for certain "interesting" subsets of the instances the run time of the 
Quantum Adiabatic Algorithm grows sub-exponentially in n. A positive result of this kind 
would greatly expand the known power of quantum computers. At the same time it is 
worthwhile to understand the circumstances under which the algorithm is doomed to fail. 
Section \2.3\ based on the paper |34j . 

How to make the Quantum Adiabatic Algorithm Fail 

Edward Farhi, Jeffrey Goldstone, Sam Gutmann and Daniel Nagaj 

The quantum adiabatic algorithm is a Hamiltonian based quantum algo- 
rithm designed to find the minimum of a classical cost function whose do- 
main has size N. We show that poor choices for the Hamiltonian can 
guarantee that the algorithm will not find the minimum if the run time 
grows more slowly than square root of N. These poor choices are nonlocal 
and wash out any structure in the cost function to be minimized and the 
best that can be hoped for is Grover speedup. These failures tell us what 
not to do when designing quantum adiabatic algorithms. 

contains two results about the necessity of local structure for successful AQC algorithms. 
First, in Section [2.3. Il l show a connection between a bad choice of the starting Hamiltonian 
in an adiabatic algorithm and an information theoretical lower bound on the running time 
of any quantum algorithm for unstructured search (see Section ll.2.3p . Second, I show in 
Section [2. 3. 21 that knowing the spectrum of the final Hamiltonian is not enough to determine 
whether finding its ground state via an AQC algorithm is feasible. The local structure 
determining how the final spectrum arises is essential. 
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2.1 Introduction to AQC 



Consider a classical optimization problem with cost function h{z) on n bits, with z = 
0, . . . , 2" — 1. I can construct a Hamiltonian for n qubits, diagonal in the computational (z) 
basis, according to this cost function as 



Af-l 



Hp=^h{z)\z) . (2.1) 



2 = 

The goal is now to find the ground state of Hp, which encodes the solution to the classical 
optimization problem with cost function h{z). If I could prepare the ground state of Hp, 
a measurement of the spins in the z-basis would produce the desired optimal assignment z 
(or one of them, if there are several). 

To prepare the ground state of Hp, a 'beginning' Hamiltonian Hp is introduced with 
a known and easy to construct ground state {gp)- The quantum computer is a system 
governed by the time dependent Hamiltonian 

Hit) = (^l-^^Hp+(^^^Hp = {l-s)Hp + sHp, (2.2) 

where the time t runs from to T. It is also sometimes useful to denote s = t/T, with 

< s < 1. This defines a path in Hamiltonian space between Hp and Hp, with the rate of 
change of H{t) controlled by T. The Schrodinger equation for a state \tp{t)) reads 

|V;(t)) = H{t) \m) ■ (2.3) 

1 now choose the state at t = to be the ground state of Hp 

\m) = \9b) (2.4) 

and run the algorithm for time T (let the system evolve according to H{t)). The power of 
the algorithm comes from the Adiabatic theorem. Suited for our investigation, in this work 
I choose to refer to its formulation due to Jeffrey Goldstone 



Theorem 2 (Adiabatic Theorem). Consider a time- dependent Hamiltonian H{t), and the 
time evolution of the state \4>o) , the ground state of H{0), from time t = to t = T. Then 

e(T) = \\\U{T,0)MO)) - \Mm\ = O (^max^^) , (2.5) 



where U{T,0) is the time evolution operator from time to T corresponding to H{t), while 
ML 

dt 



M{t) = T ll^^ll and g{t) = Ei{t) - Eo{t) is the energy gap of H{t). 



Thus, by the adiabatic theorem, if the gap is large enough, T is large enough, and the 
norm of the derivative of H{t) is bounded, \^l^{T)) will have a large component in the ground 
state subspace of Hp. A measurement of z can then be used to find the minimum of the 
classical cost function h{z) I wanted to find. The algorithm is useful if the required run 
time T does not grow exponentially with the number of qubits n. 

Rescaling the Hamiltonian by a factor E would result in a shorter required running 
time. However, this is a known tradeoff between energy and time. It is thus usual to think 
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the required resources of an AQC algorithm as T ■ \\H\\, or as the required running time for 
a rescaled Hamiltonian whose norm is \\H\\ = 1. 

Note that if the eigenvalue Xo{t) is degenerate, i.e. when we have a subspace Co{t) 
instead of a single state |0o(O) corresponding to the eigenvalue Ao(s), the formulation of 
the theorem can be generalized to involve the projection onto the subspace Cq{T) instead 
of the overlap with the single eigenstate |(/)o(r)). 

It used to be common to state the adiabatic theorem with e(T) = const, for T oc A~^, 
where A = miatg(t) is the minimum energy gap. However, several authors have recently 
rigorously investigated the sufficient conditions for the adiabatic theorem and only proved 
results in which the scaling of T with the gap was worse than [26l[l6l|9]. In this work, 
I choose to use Goldstone's formulation which for the particular Hamiltonian I investigate 
(j2.2p implies constant overlap of the time evolved state with the ground state of Hp if the 
control time depends on the gap like T oc A~^. 

Also note that other time-dependent paths in Hamiltonian space from Hb to Hp are 
possible, each of them defining a specific Quantum Adiabatic Algorithm. In fact, it has 
been shown [31j, that a path change in the Hamiltonian space can be used to beat certain 
counterexamples in which the standard, linear-path QAC requires exponential time to run. 

2.1.1 The power of AQC 

Although AQC was introduced with the hope of solving local optimization problems, the 
verdict on the practical use of AQC for this purpose is so far 'not discouraging'. Several 
numerical studies provided results that did not show exponential behavior in the required 
running time [331I111- Nevertheless, because of computer resource limitations, these could 
not be done for bit numbers greater than n = 24, while behavior change might be expected 
for much larger values of n [19]. Only very recently. Young et al. |113j utilized a different 
numerical method, investigating the gap in a QAC algorithm for the NP-complete Exact 
Cover problem using a Quantum Monte Carlo method, and found polynomial behavior of 
the gap for up to n = 128 qubits. 

On the other hand, soon after AQC was introduced, examples of cases in which a simple 
linear adiabatic algorithm failed were constructed by van Dam et al. [lOlj . However, all of 
these problems were circumvented in Farhi et al. [31] by changing the path in Hamiltonian 
space by adding a random term He to the Hamiltonian as 

H{t) = (1 - s)Hb + sHp + s(l - s)Hc. (2.6) 

In the cleverly constructed example of [lOlj , the gap above the ground state is exponentially 
small only at a single point on the interpolation path in s. Adding the term He avoids 
the bad point on the Hamiltonian path where energy levels cross, making the scaling of the 
required runtime favorable again. 

Recently, more papers investigating cases where certain adiabatic algorithms fail ap- 
peared [1081 I109j . Nevertheless, they share a common design flaw, a bad choice of the 
starting Hamiltonian. In Section 12.31 I show that poor choices for the Hamiltonian can 
guarantee that the algorithm will not find the minimum if the run time grows more slowly 
than ^/N . These poor choices are nonlocal and wash out any structure in the cost function 
to be minimized and the best that can be hoped for is Grover speedup. These failures tell 
us what not to do when designing quantum adiabatic algorithms. 

AQC is interesting on its own as an implementation of a quantum computer, as Aharonov 
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et al. in [7] (see Section I2.2p proved that its computational power is equivalent to that of 
the conventional quantum circuit model. Because of the energy gap, the computation in 
the AQC model is inherently protected against noise. The robustness and fault tolerance 
of the AQC has been recently investigated by Childs et al. [23], Jordan et al. [92], Lidar 
[61j and Lloyd |63j . The challenge today is to find more new, entirely adiabatic algorithms, 
such as the adiabatic state preparation of Aharonov and Ta-Shma [6]. 



2.2 AQC is equivalent to BQP 

Aharonov et al. proved the universality of the AQC model in [7]: 

Theorem 3 (Equivalence of AQC and the quantum circuit model). The model of adiabatic 
computation is polynomially equivalent to the standard model of quantum computation. 

One of the proof directions is straightforward. One can simulate time evolution with an 
AQC Hamiltonian on a standard quantum computer, dividing the time into small intervals 
and approximating the time evolution in each slice. The usual procedure utilizes the finite 
n approximation of the Trotter-Suzuki formula [981 IMl [SH [961 ISS] 

e'iA+B)'''= hm {ei^Mn^BM/ny ^ ^2.7) 

n— >oo V / 

and the resulting unitary transformations are applied using a circuit-based quantum com- 
puter. The precision of this approximation can be improved, as the AQC Hamiltonian 
is changing only linearly in time by integrating the Dyson series for the time-dependent 
Hamiltonian H{t). Assume that H{t) consists of two types of terms, the first contain- 
ing only terms built from trf and the second containing only terms built from (Tf. The 
Hamiltonian can then be written as 

H{t) = (l-^^H,+ {^^^H,, (2.8) 

where and contain only terms built from the respective type of operator. As shown 
in Appendix [Bl the approximation of the time evolution operator U{t + At,t) in the form 
of a product of exponentials, correct to order (At)^, is 

W2 = e-^''2H,g-i62/^.g-ic2i/,^ ^2.9) 

with 

/ t At\ At 
/ t At\ , 

h2 = (^ + ^jAt. (2.11) 

On the other hand, one can encode the evaluation of any quantum circuit U into time 
evolution with a sufficiently slowly changing time-dependent AQC Hamiltonian. The proof 
in [7] uses ideas from Kitaev's construction I described in Section 11.3.21 Consider a system 
with two registers, work and clock (jl.32p . with the work register holding n qubits qi, . . . ,qn 
and with the clock register consisting of L -|- 1 qubits cq, . . . ,cl, where L is the number of 
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gates in the circuit U . Encode the clock states states \t)^ in unary as in (jl.58p . Aharonov 
et. al. choose the initial Hamiltonian 

The term Hinput (|1.57p checks that when the clock is in the state |0)^ = IIcqOci • • • Oc^,), the 
state of the work qubits is initialized in 1 00 ... 0) , by adding an energy penalty to all states 
for which this is not true. Second, Hdock (|1.59p checks whether the clock register states are 
proper unary-clock states. Finally, in the unary clock encoding, the term 

Hclockinit = |1) {Mci (2-13) 

is a projector onto all states whose clock register is not in the state |0)^, the only legal clock 
state whose second clock qubit (ci) is in the state |0). Because the three terms Hinit are 
non-negative, the ground state of Hinit is the state annihilated by all of them: 

\(|)^n^t) = |0o(O)) = |00 . . . 0) 1 1^^ 0^^ . . . 0^^ ) . (2.14) 

The final Hamiltonian is chosen to be a modification of ()1.60p 

H final — Hprop ~\~ Hinput ~\~ Hoiocki (^-l^) 

omitting the term Hout- The term Hprop (|1.53p has a degenerate ground state subspace 
containing all the possible history states (jl.Sip for a given quantum circuit U. The ground 
state of the complete H final is the specific history state (ll.Sip for the quantum circuit U 
starting in the initial state (j2.14p . Using classical Markov chain methods, the authors then 
show that the gap above the ground state of the AQC Hamiltonian 

H{t) = (l - ^) Hinit + ^Hf.nal (2.16) 

is bounded from below by an inverse polynomial in L (the length of the computation U). 
This is possible, because it is enough to look for this gap in the subspace of legal unary 
clock states. Moreover, we know that the eigenvalue gap between the ground state subspace 
and the first excited state of Hprop is polynomially large (see Appendix [A] for details). 

The authors then encode the clock register in unary, using the 3-local construction of 
[54], and prove that 3-local AQC is universal for quantum computation. However, the cost 
of this translation from a quantum circuits to AQC is steep. First, the required runtime 
for the adiabatic algorithm is proportional to at least the inverse gap squared. Second, 
the gap of the final Hamiltonian is small, because of the penalty terms necessary for the 
3-local clock construction of [54]. The required runtime, for a rescaled Hamiltonian with 
\\H\\ = 0(1), is then claimed to be T tx L^^. However, the authors in [7J obtained this 
using a formulation of the Adiabatic theorem which was later found to be wrong, with 
a dependence of the runtime. If, in fact, I use the adiabatic theorem with A^"^^^, 
their result becomes T cx L^^~^^^. The proper adiabatic theorem I stated above has 6 = 1, 
estimating the running time for the simulation of quantum circuits by AQC as T oc L^^. 
Deift, Ruskai and Spitzer [26] later showed that one of the gap estimates in the proof (not 
in the adiabatic theorem) of [7] can be made better, resulting in the total required runtime 
scaling T oc L^"^^^^ (or, with the adiabatic theorem I use, L^^). It is possible that the 
adiabatic theorem with 6 < 1 in the condition T (X A-(2+5) is provable [46]. However, 
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for now, I stay on the safe side and use the adiabatic theorem with 5 = 1, as proved by 
Goldstone [43j. 

Later, Mizel et al. |65j showed that the required running time for simulating quantum 
circuits using AQC scales like n^L^, where L is the number of gates in the circuit, and n is 
the number of work qubits. However, this running time is for a Hamiltonian whose norm is 
ll-ff II = L, so that the real running time for a rescaled Hamiltonian scales like n^L^. Their 
construction is built on 2-particle, 4-site interactions of electrons hopping through an array 
of quantum dots. If we think of their model as implemented using qubits instead of hopping 
particles, it becomes 4-local. Also, they use the knowledge about where in their algorithm 
the gap is small, running the algorithm slower around that spot, allowing them to have 
their runtime scale as good as T tx A~^. 

When I want to keep the interactions 3-local (for qubits), much improvement is nec- 
essary. For this, I suggest using my new 3-local Hamiltonian construction of Section 
14.21 instead of [53]. The required running time for the AQC algorithm then becomes 
T (X L^+^^ = . This scales worse than Mizel's required runtime, but my Hamilto- 
nian is only 3-local for qubits. Also, when implementing my static-qubit construction using 
real hopping particles like Mizel [65], the required interactions couple 2 particles, one of 
which hops between two sites, giving a 2-particle, 3-site interaction. The improvement over 
the previous 3-local constructions is possible because my clock register encoding uses only 
constant norm penalty terms for illegal clock states, as opposed to penalty terms scaling 
like high polynomials in L as in [541 153j . Moreover, the Hamiltonians used can be further 
restricted using the results of Sections 14.4.21 and 14.51 as explained therein. 

This, however, is not the strongest result I want to present here. Using an unpublished 
result of Seth Lloyd [63], in the next Section I show that the resources necessary to perform 
universal quantum computation using 3-local AQC scale only as r • \\H\\ = O(L^log^L), 
where r is the algorithm runtime. However, in this case it is not the adiabatic theorem that 
gives the computational model its power, it is rather the underlying structure of the final 
Hamiltonian whose dynamics is a quantum walk on a line. 

2.2.1 The Hamiltonian Computer model and universality 

In this section, I look at the dynamics of the Hamiltonian ()2.16p in more detail. The result 
presented here was pointed out to me by Seth Lloyd during my Part III doctoral exa nfl and 
I thank him for the access to the write-up of his unpublished work [63j. He showed that 
even when the control time T for the change in the Hamiltonian is short, the result of the 
quantum circuit U can still be obtained from the state of the system with high probability 
after waiting for time r upper bounded by a polynomial in L, the number of gates in the 
circuit U . 

What can happen to the state \(l)init) ()2.14p when we change the Hamiltonian from Hinn 
to H final faster than the required conditions in the statement of the adiabatic theorem? Let 
me consider an extreme scenario, jumping from Hinu to Hf^nai instantly. The system will 
not be in the ground state of H final, but in a superposition of some of its eigenstates, those 
which it could transition to under the action of the Hamiltonian (I2.16|) . To see what these 
are, let me first analyze the time evolution of \(j)init) directly with H final (|2.15p . The term 
determining the dynamics of the system is Hprop, as the terms Hinput and Hdock annihilate 
the initial state and all the states H final brings it to. Consider now the subspace of states 

was made to work it out on the board, and was quite puzzled by it at the time, thinking I must have 
gotten it wrong. 
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Figure 2-1: Padding the circuit U with L gates by 2L extra identity gates. The state of the 
system after L steps stays unchanged. 

which one can transition to from \s) |0)^ using arbitrary powers of Hpyop- There is one 
such subspace for each work register basis state |s), and these subspaces are not connected 
by Hprop- As described in much more detail in Appendix [Aj after a unitary basis change 

\s)<S)\t), ^ \st)0\t)^ = {UtUt-i...U2Ui)\s)<S)\t)^, (2.17) 

the Hamiltonian Hprop becomes block diagonal. The matrix form of Hprop in each of these 
blocks is 

H'prop = ^-lB-^Po-^PL, (2.18) 

where the matrix B (see Appendix [C]) is the adjacency matrix for a line of length L + 1 
and Pq and Pl are projectors onto the endpoints of this line. This gives the system the 
dynamics of a quantum walk on a line. 

Because the blocks for different work register basis states \s) are decoupled, the specific 
initial state {(pmit) = |0) (8> |0)^ time evolved for some time r with Hprop can thus end up 
only in some superposition of the states 

l^t) = IV-t) ® \t)c = {UtUt-1 ■ ■ ■ U2U1) |00 . . . 0) \t), , (2.19) 

where is the state of the work qubits after after t gates of the quantum circuit U have 
been applied to an initial all-zero state |00. . .0). The span of all these states \^t) defines 
the subspace TCq. The state of the system |(/'(t)) after time r has no overlap with states 
outside of TIq, because the Hamiltonian Hprop is constructed so that it doesn't couple TIq 
to the subspaces generated by Hprop from some different state \s / 0) (8) |0)^. 

Note that the state contains the result of the quantum circuit U in its work register. 
I would like to find the system with the computation already done, so I boost the probability 
by padding the initial circuit as in Figure [2^ After the original L gates of the circuit, let 
me add 2L extra identity gates, taking L — > 3L. Now, all the states \^t) with t > L/3 
contain the output of the circuit U in their work register, as the work register of is 
[/ 1 00 ... 0) and stays unchanged for all t > L/3 because I padded the circuit with identity 
gates. 

The time evolution of {(pmit) with Hprop corresponds to a continuous-time quantum 
walk on a line with specific boundary conditions. I can get rid of the boundary conditions 
by wrapping the walk around on a circle, which simplifies the analysis, given in detail in 
Appendix O There I show that this walk mixes so fast, that when I choose a time r 
uniformly at random between and Ti (x L log^ L the probability to find the system in a 
state \t)^ with t > L/3 is close to 2/3. Therefore, when I let the state \4>init) evolve for a 
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time chosen uniformly at random between and 0(L log^L) and measure the clock qubit 
ci^j^, I will obtain the result 1 with probability close to 2/3. This signifies 'success', as it 
implies that the clock register of the state I measured was one of the states with t > L/3. 
The work qubits then contain the result of the quantum computation U. 

To conclude, when I change Hinu to H final quickly, the resulting state of the system will 
be some superposition of states \^t), close to |^o)- When I then let the system evolve for 
extra time r < 0(-L log^ L) with the Hamiltonian Hfi^ai, I can measure the clock register, 
and with high probability I will obtain a result that tells me that the work register now 
contains the result of the quantum circuit U. The overall resources necessary for this 
procedure then scale like r • \\H\\ = 0{L^ log^ L), as the Hamiltonian I use has norm 0{L). 

Let me now formalize the model of computation I just described. 

Definition 11 (Hamiltonian Computer (HC) model). A Hamiltonian Computer aiming to 
simulate a quantum circuit U with L gates is defined by a system consisting of two registers, 
work and clock, a k-local Hamiltonian H, a simple to prepare initial state {ipo), and a bound 
on the required evolution time. The clock register states must be constructed in such a way 
that a measurement distinguishing \t < L/3)^ (failure) and \t > (success) is simple. 

The computation procedure goes as follows: 

1. Prepare the system in the initial state |V'o); ihe ground state of a k-local Hamiltonian 
Ho. 

2. Turn off the initial Hamiltonian Hq, and let the system evolve with a time independent 
Hamiltonian H for a time r chosen uniformly at random between and poly{L). 

3. Make a measurement determining whether the clock register is in a state \t > L/3), 
which signifies success. The probability to obtain such state is greater than p = ^. 
The result of the quantum circuit U is then in the work register of the system. Restart 
otherwise. 

Feynman's Hamiltonian Computer with a unary clock \t)^ = |0 . . . 010 . . . 0))c (see Sec- 
tion [L3T]) is one of the possible implementations of a HC. The initial Hamiltonian Hq is one 
whose ground state is IV'o) = |Ogi • • • Oq„) |0co . . . Ocl), and the time independent Hamilto- 
nian H is given by (jl.SSp . To determine success/failure, I need to measure the first third 
of the clock qubits. When I don't find the spin up there, it means I have succeeded. In the 
case of Feynman's Hamiltonian, the Hamiltonian itself is 4-local, and Lloyd [63] bounds the 
required runtime from above by L^. As seen in the previous Section, Aharonov et al. f7] use 
the 3-local Hamiltonian H final of Kempe and Regev [Mj to simulate quantum circuits using 
AQC, and require resources of the order L^^^^*^. However, when I encode the clock register 
and the clock state transitions as in Section 14.51 the gap of H final whose norm is 0{L) 
scales like L~^. The mixing of the quantum walk is then such that the required running 
time becomes only O(-Llog^L). When I rescale my Hamiltonian to = 1, the resources 
necessary for this HC model are 0(L^ log^L). For more detail on this construction, see 
Section 14. 5[ 

What remains is the question how is one protected against noise in such Hamiltonian 
Computer models. In AQC, it was the energy gap. However, as I have shown, the energy 
gap computed by Aharonov et al. [7] is not relevant, because the transitions to higher 
momenta eigenstates of the quantum walk still contain the result of the computation. The 
relevant energy barrier is the energy cost to transition from the subspace TLiegai to Ti-j-^g^i, 
in which case we would lose the computation. For this purpose, I use a clock-checking 
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Hamiltonian like (jl.59p . whose norm is 0{L). Lloyd [63] showed, that the energy gap for 
such Hamiltonian scales like L~^. In summary, when simulating quantum circuits in a HC 
model, one requires resources r • \\H\\ = O(L^log^-L), and is protected by an energy gap 
scaling like L~^, where L is the number of gates in the circuit. This energy barrier to 
decoherence can be increased, at the cost of increasing the locality of interactions, by using 
error correcting codes like the ones developed for AQC by Jordan et al. [92] . 



2.3 How to make the Quantum Adiabatic Algorithm Fail 

After looking at the power of AQC in general and showing how it can be used to simulate 
any quantum circuit effectively, let me now return to its original intended purpose: solving 
locally constrained optimization problems. 

Recently, a paper announcing the failure of AQC in certain cases appeared [108j , putting 
considerable effort into showing that the gap in a specific AQC algorithm is exponentially 
small. However, the authors make a bad choice of the starting Hamiltonian Hb in ()2.2p . 
In this section I prove some general results which show that with certain choices of Hb or 
Hp the algorithm will not succeed if T is o{y/N), that is T/y/N ^ as ^ oo, so that 
improvement beyond Grover speedup is impossible. I view these failures as due to poor 
choices for Hb and Hp, which teach us what not to do when looking for good algorithms. 
I guarantee failure by removing any structure which might exist in h{z) from either Hp or 
Hp. By structure I mean that z is written as a bit string and both Hp and Hp are sums of 
terms involving only a few of the corresponding qubits. After the preprint [M! was posted 
on the quant-ph archive, two similar papers [109^ [45] concerning this topic also appeared. 

In Section 12.3.11 1 show that regardless of the form of h{z) if Hp is a one dimensional 
projector onto the uniform superposition of all the basis states \z), then the quantum 
adiabatic algorithm fails. Here all the \z) states are treated identically by Hp so any 
structure contained in h{z) is lost in Hp. In Section [2.3.21 1 consider a scrambled Hp that I 
get by replacing the cost function h{z) by h(Tr{z)) where vr is a permutation of to — 1. 
Here the values of h{z) and h{iT{z)) are the same but the relationship between input and 
output is scrambled by the permutation. This effectively destroys any structure in h{z) and 
typically results in algorithmic failure. 

The quantum adiabatic algorithm is a special case of Hamiltonian based continuous 
time quantum algorithms, where the quantum state obeys (j2.3p and the algorithm consists 
of specifying H(t), the initial state \ip{0)), a run time T and the operators to be measured 
at the end of the run. In the Hamiltonian language, the Grover problem can be recast as 
the problem of finding the ground state of 

H^ = EiI-\w){w\), (2.20) 

where w lies between and N — 1. The algorithm designer can apply H^j, but in this 
oracular setting, w is not known. In [37] the following result was proved. Let 

H{t) = Hpit) + H^, (2.21) 

where Hp is any time dependent "driver" Hamiltonian independent of w. Assume also 
that the initial state \ip{0)) is independent of w. For each w we want the algorithm to be 
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successful, that is \^{T)) = \ w). It then follows that 




(2.22) 



The proof of this result is a continuous-time version of the BBBV oracular proof [171 ■ 
proof techniques in this paper are similar to the methods used to prove the result just 
stated. 



2.3.1 General search starting with a one-dimensional projector 

In this section I consider a completely general cost function h{z) with z = 0, . . . , — 1. 
The goal is to use the quantum adiabatic algorithm to find the ground state of Hp given 
by dlT]) with H{t) given by ([22]). Let 

7^ E 1^) (2-23) 

2 = 

be the uniform superposition over all possible values z. If I pick 

Hb = E{I-\s){s\) (2.24) 

and IV'(O)) = \s), then the adiabatic algorithm fails in the following sense: 

Theorem 4. Let Hp be diagonal in the z basis with a ground state subspace of dimension 
k. Let 

H{t) = (1 - t/T)E (I - \s) (s|) + {t/T)Hp. 

Let P be the projector onto the ground state subspace of Hp and let b > be the success 
probability, that is, b = {'4'{T)\ P\i{j{T)). Then 

^ b [n 2Vb 

T > — W —. 

- E \ k E 

Proof. Keeping Hp fixed, introduce — 1 additional beginning Hamiltonians as follows. 
For X = N — 1 lei Vx he a unitary operator diagonal in the z basis with 

{z\ Vx\z)= e^^^^^l^ 

and let 



iV-l 

'^^ 2=0 



so that the form an orthonormal basis. Note also that 

|x = 0) = \s) . 

We now define 

Hxit) = (1 - tlT)Ei^ - \x) (x|) + {tlT)Hp, 

with corresponding evolution operator Ux{t2,ti). Note that H{t) above is HQ{t) with the 
corresponding evolution operator Uq. For each x let the system evolve with Hx{t) from 
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the ground state of Hx{0), which is \x). Note that = VxHqVx and Ux = VxUqVx ■ Let 
\fx) = Ux(T, 0) \x). For each x the success probabihty is {fx \ P\fx), which is equal to b since 
P commutes with Vx- The key point is that if I run the Hamiltonian evolution with Hx 
backwards in time, I would then be finding x, that is, solving the Grover problem. However, 
this should not be possible unless the run time T is of order s/N. 



Let Ur be the evolution operator corresponding to an x-independent reference Hamil- 
tonian 

Hnit) = (1 - t/T)E + {t/T)Hp. 

Let \gx) = '^-^Ifx) be the normalized component of \ fx) in the ground state subspace of 
Hp. Consider the difference in backward evolution from \gx) with Hamiltonians Hx and 
Hji, and sum on x, 

Sit) = ^ ||c/t(r,t) \gx) - Ul,{T,t) \gx)f . 

X 

Clearly S{T) = 0, and 

^(0) = ^\\ui{T,0)\gx) -Ul{T,0)\gx)f . 

X 

Now \gx) = Vb\fx)+VT^\f^) where \f^) is orthogonal to 1/^). Since i7j(T,0) \fx) = \x), 
S{0) = Y,\\^\x) + VT^ 



where for each x, \x^) and \ix) are normalized states with \x^) orthogonal to \x). Since Hr 
commutes with Hp, \ix) = Uj^(T,0) \gx) is an element of the fc-dimensional ground state 
subspace of Hp. Then 

S{0) = 2N -^^Vb{x\ix) + VT^{x^\ix) + C.C. 

X 

> 2N - 2\/hY^Yx\ix) - 2NVl - b. 

x 

Choosing a basis for the k dimensional ground state subspace of Hp and writing 

\ix) = axi \Gi) H \-axk \Gk) gives 



(2.25) 



= VNk. 



Thus 



S{0) > 2N{1 - VT^) - 2VbVNk. 



(2.26) 
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Using the Schrodinger equation, one can find the time derivative of S{t): 



+ c.c. 



Now 



X 

-i^{gx\ U^{T,t)[H^{t) - HR{t)]ul{T,t) \g.)+c.c. 

X 

-2Im^(l - t/T)E{g,\ U^{T,t) \x) {x\ ul{T,t) \g,) . 



< 2i5;(l - t/T) \{9x\ U^{T,t) \x) {x\ ul{T,t) \g^) 

X 

< 2E{l-t/T)Y\{x\ul,iT,t)\g, 



The same technique as in (I2.25P gives 



< 2E{l-t/T)\^. 



Therefore 



T 



dt 



S{t) 



dt < ETVNk. 



Now S{0) < S{T) + \^tS{t)\ dt and S{T) = so 

5(0) < ETVWk. 

Combining this with (j2.26p gives 

ETVNk > 2N{1 - Vl - b) - 2Vb\^, 
which implies what I wanted to prove: 



E\l k E 



□ 



How to interpret Theorem UP The goal is to find the minimum of the cost function h{z) 
using the quantum adiabatic algorithm. It is natural to pick for Hb a Hamiltonian whose 
ground state is |s), the uniform superposition of all \z) states. However if I pick Hb to be 
the one dimensional projector E{I — \s) (s|) the algorithm will not find the ground state 
if T/^/N goes to as goes to infinity. The problem is that Hb has no structure and 
makes no reference to h{z). My hope is that the algorithm might be useful for interesting 
computational problems if Hb has structure that reflects the form of h{z). 

Note that Theorem [5] explains the algorithmic failure discovered by Znidaric and Horvat 
|108j for a particular set of h{z). 
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For a simple but convincing example of the importance of the choice of Hb, suppose I 
take a decoupled n bit problem which consists of n clauses each acting on one bit, say for 
each bit j 

so 



h{z) = zi + Z2 + ■■■ + Zn- (2.27) 
Let me pick a beginning Hamiltonian reflecting the bit structure of the problem, 

The ground state of Hb is |s). The quantum adiabatic algorithm acts on each bit indepen- 
dently, producing a success probability of 

p={l-q{T)r, 

where q{T) — > as T — > oo is the transition probability between the ground state and the 
excited state of a single qubit. As long as nq{T) const. I have a constant probability 
of success. This can be achieved for T of order ^/n, because for a two level system with a 
nonzero gap, the probability of a transition is q{T) = 0(r~^). (For details, see Appendix 
IEI ) However, I know from Theorem [4] that a poor choice of Hb would make the quan- 
tum adiabatic algorithm fail on this simple decoupled n bit problem by destroying the bit 
structure. 

Next, suppose the satisfiability problem I am trying to solve has clauses involving say 3 
bits. If clause c involves bits ic, jc and kc I may define the clause cost function 

^ (^\ — i ^ if -Zfec satisfy clause c, 

\ 1 otherwise. 

The total cost function is then 

h{z)=Y,hc{z). 



To get Hb to reflect the bit and clause structure I may pick 

Hb,c = l[a- 4'^) + (1 - ^^'^) + (1 - 4'^^)" 

with 

Hb = Y,Hb,c. (2.29) 

c 

In this case the ground state of Hb is again \s). With this setup, Theorem|l]does not apply. 
I did a numerical study of a particular satisfiability problem. Exact Cover. For this 
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Figure 2-2: Median required run time T versus bit number. At each bit number there 
are 50 random instances of Exact Cover with a single satisfying assignment. I choose the 
required run time to be the value of T for which quantum adiabatic algorithm has success 
probability between 0.2 and 0.21. For the projector beginning Hamiltonian I use (j2.24p 
with E = n/2. The plot is log-linear. The error bars show the 95% confidence interval for 
the true medians. 
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problem if clause c involves bits ic, jc and kc, the cost function is 




if Zi^ + Zj^ + zk, = l, 
otherwise. 



Some data is presented in FIG. 1. Here we see that with a structured beginning Hamiltonian 
the required run times are substantially lower than with the projector Hb- 



2.3.2 Search with a scrambled problem hamiltonian 

In the previous section I showed that removing all structure from dooms the quantum 
adiabatic algorithm to failure. In this section I remove structure from the problem to be 
solved (Hp) and show that this leads to algorithmic failure. Let h{z) be a cost function 
whose minimum I seek. Let tt be a permutation of 0, 1, . . . , AT — 1 and let 

I will show that no continuous time quantum algorithm (of a very general form) can find the 
minimum of /il'^l for even a small fraction of all tt if T is o{^/N). Classically, this problem 
takes order TV calls to an oracle. 

Without loss of generality let h{0) = 0, and h{l), h{2), . . . , h{N — 1) all be positive. For 
any permutation tt of 0, 1, . . . , A?^ — 1, define a problem Hamiltonian Hp^T^, diagonal in the 
z basis, as 

N-l N-l 

Hp,. = J2 ^'"'(^) \z) {z\ = J2 ^(^) • 

2=0 z=0 

Now consider the Hamiltonian 

H.{t) = HD{t)+c{t)Hp,. (2.30) 

for an arbitrary 7r-independent driving Hamiltonian H£){t) with \c{t)\ < 1 for all t. Using this 
composite Hamiltonian, evolve the 7r-independent starting state |V'(0)) for time T, reaching 
the state \if:^{T)). This setup is more general than the quantum adiabatic algorithm since 
I do not require Hpi{t) or c{t) to be slowly varying. Success is achieved if the overlap of 
|V'7r(r)) with |7r(0)) is large. 
I first show 



Lemma 1. 

5^ II \MT)) - \i^.'{T)) II' < Ah*TN\sfN^, (2.31) 

where the sum is over all pairs of permutations vr, vr' that differ by a single transposition 
involving 7r(0), and h* = h{z)^/ {N — 1). 

Proof. For two different permutations vr and vr' let \tpT^{t)) be the state obtained by evolving 
from \ip{0)) with and let \^p-,^l{t)) be the state obtained by evolving from |V'(0)) with 

Hn'- 
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Now 



dt 



dt 



(MtMn'it)) +C.C. 



< 2 



{Mt)\iHnit)-H,'{t))\i^^it)) 



Consider the case when tt and tt' differ by a single transposition involving 7r(0). Specifically, 
tt' = TT o (a ■f-^ 0) for some a. Now if 7r(0) = i and 7r(a) = j, we have 7r'(0) = j and 7r'(a) = i. 
Therefore, since h{0) = 0, 

Hp,^ - Hpy = c{t)h{a) (li) (il - \i) = c(i)/.(a) (|7r(G)) (7r(a)| - |7r'(a)> <7r'(a)|) , 

so that 

_d 
At 



E l^-W)-IV'T'(i)) < 2\c{t)\Y,Ka) (V^,(i)|(Ka))(7r(a)|-|7r'(a)><7r'(a)|)|^,,(i)) 



7r,7r 



This further simphfies to 



II IV'-(^)) - If ^ 2EM«) (|(V.(i)k(a))| + |<7r'(a)|V.'(i)>|) 



2EE^(") (V'.(t)lvr(a)) +2^^%) (7r'(a)|V.' (0) 

4EE^(«)|(^-(*)i^(«)) 



- 4E JE^(")^ = 



7r V a 



where I used the Cauchy-Schwartz inequahty to obtain the last line. Integrating this in- 
equality for time T, I obtain the result I wanted to prove, 

E||lV'.(r))-KKr))f < 4rriV!Viv^, 

7r,7r' 

where the sum is over tt and tt' differing by a single transposition involving 7r(0). □ 



Next I establish 

Lemma 2. Suppose \2) , \L) are orthonormal vectors and {ipi\i) 
vectors where i = 1, . . . , L. Then for any normalized \(p), 

E II iV'i) - 1^) >bL- 2VI. 

i=l 



> h for normalized 



(2.32) 
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Proof. Write 

i i 

> j;|(i|V.)f -2 5]|(z|^.)||(ib) 

i i 

and use the Cauchy-Schwartz inequality to obtain 



II m - II > &i - |(ii^.)| l^^i^) 

> bL- 2Vl. 



□ 



Let me now state the main result of this Section. 



Theorem 5. Suppose that a continuous time algorithm of the form ()2.30p succeeds with 

2 

probability at least b, i.e. ('i/'7r(r)|7r(0)) > b, for a set of eN\ permutations. Then 



T > 



16/i* 



^JN -I 



e^/7j2 
4h* 



(2.33) 



Proof. For any permutation vr, there are A'^ — 1 permutations vr^ obtained from vr by first 
transposing and a. For each vr let Sj^ be the subset of those — 1 permutations on which 
the algorithm succeeds with probability at least b. Any such permutation appears in exactly 
A^ — 1 of the sets 6-,^ , therefore 

Y,\Sn\ = {N-l)eN\. 
Let M be the number of sets with 15,^1 > f (-^ ~ Now 

|5.|>f(iV-l) |cS.|<f(7V-l) 

E|5^| < M{N-l) + {Nl-M)^{N-l), 

TT 

(A^-l)eA^! < M(A^- 1) + A^!|(A^- 1), 

so M > |A^!, i.e. at least |A^! of the sets S-,^ must contain at least |(A^ — 1) permutations 
on which the algorithm succeeds with probability at least b. For the corresponding vr, I 
have 

II \MT)) - \^.'ST)) f > b'-{N - 1) - 2J'-{N-l). 



by Lemma 2. (Note that the algorithm is not assumed to succeed with probability 6 on vr.) 
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Figure 2-3: The scaled ground state energy E/n for a quantum adiabatic algorithm Hamil- 
tonian of a decoupled problem. The lowest curve corresponds to the original decoupled 
problem. The upper "triangular" curves correspond to single instances of the n-bit decou- 
pled problem, where the problem Hamiltonian was scrambled. 



Since there are at least such vr, 



II \MT)) - If > \n\ (b'-{N - 1) - 2^^(iV-l)) , 

where the sum is over all permutations vr and vr' which differ by a single transposition 
involving 7r(0). Combining this with Lemma 1 I obtain 



16h* Ah* ' 

which is what I wanted to prove. □ 

What I have just shown is that no continuous time algorithm of the form (j2.30p can 
find the minimum of Hp^.,^ with a constant success probability for even a fraction eA^! of all 
permutations vr if T is o{\fN). A typical permutation vr yields an -ffp^jr with no structure 
relevant to any fixed Hjy and the algorithm cannot find the ground state of i^p^-n- efficiently. 

To illustrate the nature of this failure for the quantum adiabatic algorithm for a typical 
permutation, consider again the decoupled n bit problem with h{z) given by (j2.27p and Hb 
given by (|2.28|) . The lowest curve in FIG. 2 shows the ground state energy divided by n 
as a function of t. (Since the system is decoupled this is actually the ground state energy 
of a single qubit.) I then consider the n bit scrambled problem for different values of n. 
At each n I pick a single random permutation tt of 0, . . . , (2" — 1) and apply it to obtain a 
cost function /i(7r~^(z)) while keeping Hb fixed. The ground state energy divided by n is 



58 



now plotted for n = 9, 12, 15 and 18. From these scrambled problems it is clear that if I let 
n get large the typical curves will approach a triangle with a discontinuous first derivative 
at t = r/2. For large n, the ground state changes dramatically as t passes through T/2. 
In order to keep the quantum system in the ground state we need to go very slowly near 
t = T/2 and this results in a long required run time. 

2.3.3 Summary 

I have shown two main results about the performance of the quantum adiabatic algorithm 
when used to find the minimum of a classical cost function h{z) with z = 0, . . . , N — I. 
Theorem |4] says that for any cost function h{z), if the beginning Hamiltonian is a one 
dimensional projector onto the uniform superposition of all the \ z) basis states, the algorithm 
will not find the minimum of /i if T is less then of order This is true regardless of how 
simple it is to classically find the minimum of h{z). 

In Theorem [5] I start with any beginning Hamiltonian and classical cost function h. 
Replacing h{z) by a scrambled version, i.e. /lW(z) = h{TT{z)) with TT a permutation of to 
A'^ — 1, will make it impossible for the algorithm to find the minimum of h^'^'^ in time less 
than order \/iV for a typical permutation vr. For example suppose we have a cost function 
h{z) and have chosen Hb so that the quantum algorithm finds the minimum in time of 
order logA^. Still scrambling the cost function results in algorithmic failure. 

These results do not imply anything about the more interesting case where Hb and Hp 
are structured, i.e., sums of terms each operating only on several qubits. 
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Chapter 3 



Matrix Product States 
on Infinite Trees 

When can one actually simulate local Hamiltonians or find their ground states efficiently? 
It is possible for some Hamiltonians in ID and in tree geometry, by using an approximate 
method based on Matrix Product States. This is what I show in this Chapter, based on the 
paper [69] 

The Quantum Transverse Field Ising Model on an Infinite Tree 
from Matrix Product States 

Daniel Nagaj, Edward Farhi, Jeffrey Goldstone, Peter Shor and Igor Sylvester 

We give a generalization to an infinite tree geometry of Vidal's infinite 
time-evolving block decimation (iTEBD) algorithm [106] for simulating an 
infinite line of quantum spins. We numerically investigate the quantum 
Ising model in a transverse field on the Bethe lattice using the Matrix 
Product State ansatz. We observe a second order phase transition, with 
certain key differences from the transverse field Ising model on an infinite 
spin chain. We also investigate a transverse field Ising model with a specific 
longitudinal field. When the transverse field is turned off, this model has a 
highly degenerate ground state as opposed to the pure Ising model whose 
ground state is only doubly degenerate. 

This Chapter is organized as follows. Section 13.21 is a review of the MPS ansatz and 
contains its generalization to the tree geometry. In Section 13.31 I review the numerical 
procedure for unitary updates and give a recipe for applying imaginary time evolution 
within the MPS ansatz. In Section 13.41 I adapt Vidal's iTEBD method for simulating 
translationally invariant one-dimensional systems to systems with tree geometry. Section 
13.51 contains my numerical results for the quantum Ising model in a transverse field for 
translationally invariant systems. In Section 13.5.11 I first test my method for the infinite 
line, and then in Section [3.5.21 1 present new results for the infinite tree. I turn to the not 
00 model in Section 13.61 and show that my numerics work well for this system even when 
there is a high ground state degeneracy. In Section 13.71 1 investigate the stability of my 
tree results and conjecture that they may be good approximations to a local description far 
from the boundary of a large finite tree system. 
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Figure 3-1: The Bethe lattice (infinite Cayley tree). 



3.1 Introduction 

The matrix product state (MPS) description [29] [75] has brought a new way of approaching 
many-body quantum systems. Several methods of investigating spin systems have been de- 
veloped recently combining state of the art many-body techniques such as White's Density 
Matrix Renormalization Group [110] f83] (DMRG) with quantum information motivated in- 
sights. Vidal's Time Evolving Block Decimation (TEBD) algorithm |105j |104j uses MPS 
and emphasizes entanglement (as measured by the Schmidt number), directing the compu- 
tational resources into that bottleneck of the simulation. It provides the ability to simulate 
time evolution and it was shown that MPS-inspired methods handle periodic boundary 
conditions well in one dimension [103J, areas where the previous use of DMRG was lim- 
ited. TEBD has been recast into the language of DMRG in [lllj |25j and adapted to finite 
systems with tree geometry in [89]. DMRG is especially successful in describing the prop- 
erties of quantum spin chains, the application of basic DMRG-like methods is limited for 
quantum systems with higher dimensional geometry. New methods like PEPS |102j |68j 
generalize MPS to higher dimensions, opening ways to numerically investigate systems that 
were previously inaccessible. 

Here I am interested in investigating infinite translationally invariant systems. Several 
numerical methods to investigate these were developed recently. The iTEBD algorithm 
[106j (see also Sec J3.4p is a generalization of TEBD to infinite one-dimensional systems. A 
combination of PEPS with iTEBD called iPEPS [50] provides a possibility of investigating 
infinite translationally invariant systems in higher dimensions. My contribution is a method 
to investigate the ground state properties of infinite translationally invariant quantum sys- 
tems on the Bethe lattice using imaginary time evolution with Matrix Product States. The 
Bethe lattice is an infinite tree with each node having three neighbors, as depicted in Fig l3-ll 
It is translationally invariant in that it looks the same at every vertex. This geometry is 
interesting, because of the following connection to large random graphs with fixed valence. 
Moving out from any vertex in such a random graph, you need to go a distance of order 
logn, where n is the number of vertices in the graph, before you detect that you are not on 
the Bethe lattice, that is, before you see a loop. 

I choose to investigate the quantum transverse field Ising model on the Bethe lattice. 
Note that I work directly on the infinite system, never taking a limit. First I test the iTEBD 
method on a system with a known exact solution, the infinite line. Then I turn to the Bethe 
lattice with the new method I provide. In both cases, the Hamiltonian is given by 




(3.1) 
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where the sum over i is over all sites, and the sum over is over all bonds (nearest 

neighbors). I show that imaginary time evolution within the MPS ansatz provides a very 
good approximation for the exact ground state on an infinite line, resulting in nearly correct 
critical exponents for the magnetization and correlation length as one approaches the phase 
transition. I obtain new results for the quantum Ising model in transverse field on the 
infinite tree. Similarly to the infinite line, I observe a second order phase transition and 
obtain the critical exponent for the magnetization, /3t ~ 0.41 (different than the mean-field 
result). However, the correlation length does not diverge at the phase transition for this 
system and I conjecture that it has the value l/ln2. 

I also investigate a model where besides an antiferromagnetic interaction of spins I add 
a specific longitudinal field |cj* for each spin: 

I choose the longitudinal field in such a way that the interaction term in the computational 
basis takes a simple form, |00) (00| j^, giving an energy penalty to the |00) state of neighboring 
spins. (I follow the usual convention that spin up in the z-direction is called 0.) I call it the 
NOT 00 model accordingly. This model is interesting from a computational viewpoint. The 
degeneracy of the ground state of //not oo at /i = is high for both infinite line and infinite 
tree geometry of interactions. I am interested in how my numerical method deals with this 
case, as opposed to the double degeneracy of the ground state of (13. 1[) at /i = 0. I do not 
see a phase transition in this system as I vary J and h. 

3.2 Matrix Product States 

If one's goal is to numerically investigate a system governed by a local Hamiltonian, it is 
convenient to find a local description and update rules for the system. A Matrix Product 
State description is particularly suited to spin systems for which the connections do not 
form any loops. Given a state of this system, I will first show how to obtain its MPS 
description, and then how to utilize this description in a numerical method for obtaining 
the time evolution and approximating the ground state (using imaginary time evolution). I 
begin with matrix product states on a line (a spin chain), and then generalize the description 
to a tree geometry. In 13.31 I give a numerical method of updating the MPS description for 
both real and imaginary time simulations. 

3.2.1 MPS for a spin chain 

Given a state {ip) of a chain of n spins 

IV') = Yl c.. |si)i...|Si)JSi+i)i+i...|Sn)„, (3.3) 

...SiSi + l... 

I wish to rewrite the coefficients Cs^^...^s„ as a matrix product (see [77] for a review of MPS) 

C...,s,,s,+i,... - 2^ ...A^ ''a,b \ ^b,c -^c ••• W-^) 

...abc... 
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Figure 3-2: Two successive Schmidt decompositions on a line allow me to find the T tensor 
for the marked site and the two A vectors for the bonds coming out of it. 



using n tensors F^*) and n — 1 vectors A^*-*. The range of the indices a,b, . . . will be addressed 
later. After decomposing the chain into two subsystems, one can rewrite the state of the 
whole system in terms of orthonormal bases of the subsystems. A*^*) is the vector of Schmidt 
coefficients for the decomposition of the state of the chain onto the subsystems 1 . . .i and 
i + 1 . . . n. 

In order to obtain the A's and the F's for a given state \^|J), one has to perform the 
following steps. First, perform the Schmidt decomposition of the chain between sites i — 1 
and i as 

IV') = i;Va)w-lAr^)|<^a),.„,„, (3.5) 
a=l 

where the states on the left and on the right of the division form orthonormal bases required 
to describe the respective subsystems of the state The number Xi-i (the Schmidt 

number) is the minimum number of terms required in this decomposition. The Schmidt 
decomposition for a split between sites i and i + 1 gives 

IV') = f;|e.)wA«|06).+i,...,.. (3.6) 

b=l 

These two decompositions (see FIG J3-2P describe the same state, allowing me to combine 
them to express the basis of the subsystem i, . . . ,n using the spin at site i and the basis of 
the subsystem i + 1, . . . ,n as 

s=0,l 6=1 

where I inserted the a[*'' for convenience. This gives me the tensor F^*). It carries an index s 
corresponding to the state \s) of the i-th spin, and indices a and b, corresponding to the two 
consecutive divisions of the system (see FIG I3-2p . Because {(pa) (and \6b)) are orthonormal 
states, the vectors A and tensors F obey the following normalization conditions. From (j3.6p 
I have 

EAp = l, (3.8) 

6=1 
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Figure 3-3: The three Schmidt decompositions on a tree required to obtain the T tensor for 
the marked site and the three A vectors for the bonds emanating from it. 



while (13. 7p imphes 

{4>a'\<Pak...,n = ^ E rj^f Af^r^'^A^ = 5.,, , (3.9) 
s=0,l b=l 

and 

moth.... = E E Ar^)r«r Ari)r«/ = 6,,, . (3.10) 

s=0,l a=l 



3.2.2 MPS on Trees 

Matrix Product States are natural not just on chains, but also on trees, because these 
can also be split into two subsystems by cutting a single bond, allowing for the Schmidt- 
decomposition interpretation as described in the previous section. The Matrix Product 
State description of a state of a spin system on a tree, i.e. such that the bonds do not form 
loops, is a generalization of the above procedure. Tree-tensor-network descriptions such as 
the one given here have been previously described in [89j . 

(k) 

Specifically, for the Bethe lattice with 3 neighbors per spin, I introduce a vector Xa^. for 
each bond k and a four-index (one for spin, three for bonds) tensor ral',l],am for each site i. 
I can then rewrite the state \ip) analogously to (|3.3p . (|3.4p as 

IV')=f n EAi?)( n Eri?,al,a„)|--->k.>|...>, (3.11) 

where a/, Um, CLn are indices corresponding to the three bonds /, m and n coming out of site 
i. Each index ai appears in two F tensors and one A vector. To obtain this description, one 
needs to perform a Schmidt decomposition across each bond. This produces the vectors 
A^'^. To obtain the tensor F^*^ for site i, one needs to combine the three decompositions 
corresponding to the bonds of site i as depicted in Fig l3-31 Analogously to (13. 7p . expressing 
the orthonormal basis for the first subsystem marked in Fig l3-3l in terms of the state of the 
spin |sj) and the orthonormal bases for the latter two subsystems in Fig J3-31 one obtains 
the tensor Fi*^;a^,a„ for site i. 

The normalization conditions for a MPS description of a state on a tree are analogous 
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to (I3:8]1- (I3TTO . I have 



E^if = l' (3-12) 

Xk Xl 

VVVrW'^* \(fc)2\(/)2p{i),s Co|o^) 

^ au^aua^i'^ak ^ai ak,ai,am "am,a^, , yo.i^oj 

s=0,l afe=l a;=l 

and two other variations of (j3.13p with A;, / and m interchanged. 

3.3 Simulating Quantum Systems with MPS 

I choose to first describe the numerical procedures for a chain of spins. Then, at the end of 
the respective subsections, I note how to generahze these to tree geometry. 

3.3.1 Unitary Update Rules 

The strength of the MPS description of the state hes in the efficient apphcation of local 
unitary update rules such as U = e~*"^^* (where A is an operator acting only on a few 
qubits). First, I describe the numerical procedure in some detail, and then, in the next 
Section, discuss how to modify the procedure to also implement imaginary time evolution. 

Given a state IV') as a Matrix Product State, I want to know what happens after an 
application of a local unitary. In particular, for a 1-local U acting on the i-th spin, it suffices 
to update the local tensor 

vff^U!,T^^f. (3.14) 

The update rule for an application of a 2-local unitary V acting on neighboring spins i and 
i + 1, requires several steps. First, using a larger tensor 

e^:* = Ai-i) (rS^A?r(;:^)-*) x^t'\ (3.i5) 
b 

I rewrite the state IV') as 

a,c s,t 

After the application of V , the tensor in the description of changes as 

©afc^EKv'®^/. (3.17) 

s't' 

One now needs to decompose the updated tensor to obtain the updated tensors F^*), 
F*-*^^^ and the vector A*-*^. I use the indices a, s and c, t of to introduce combined indices 
(as) and (ct) and form a matrix T(^as),{ct) with dimensions 2xi^i x 2xi+i as 

T{as),{ct) = 0a,c- (3-18) 
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Using the singular value decomposition (SVD), this matrix can be decomposed into T = 
QAW, where Q and W are unitary and A is a diagonal matrix. In terms of matrix elements, 
this reads 

T{as),{ct) = Q{as),6-f^fe,f)^6,{ci)- (3.19) 
b 

The diagonal matrix D = diag(A^*^) gives me the updated Schmidt vector A^*\ The updated 
tensors F^*) and r^*^^) can be obtained from the matrices Q, W and the definition of (|3.15p 
using the old vectors A^*"^-* and A*-*^^^ which do not change with the application of the local 
unitary V. After these update procedures, the conditions (j3.8p - (j3.10p are maintained. 

The usefulness/succintness of this description depends crucially on the amount of entan- 
glement across the bipartite divisions of the system as measured by the Schmidt numbers 
Xi- To exactly describe a general quantum state {ip) of a chain of n spins, the Schmidt 
number for the split through the middle of the chain is necessarily Xn/2 = 2""/^. Suppose I 
start my numerical simulation in a state that is exactly described by a MPS with only low 
Xi's. The update step described above involves an interaction of two sites, and thus could 
generate more entanglement across the i,i + 1 division. After the update, the index b in 
a[*^ would need to run from 1 to 2xi to keep the description exact (unless Xi already is at 
its maximum required value Xi = 2™"^^*'"~*J'). This makes the number of parameters in the 
MPS description grow exponentially with the number of update steps. 

So far, this description and update rules have been exact. Let me now make the de- 
scription an approximate one (use a block-decimation step) instead. First, introduce the 
parameter x, which is the maximum number of Schmidt terms I keep after each update 

(i) 

step. If the amount of entanglement in the system is low, the Schmidt coefficients A^ 
decrease rapidly with b (I always take the elements of A sorted in decreasing order). A MPS 
ansatz with restricted Xi = X '^iH hopefully be a good approximation to the exact state {ip). 
However, I also need to keep the restricted x throughout the simulation. After a two-local 
unitary update step, the vector A^*^ can have 2% entries. However, if the b > x entries in 
a[*^ after the update are small, I am justified to truncate A^*^ to have only x entries and 
multiply it by a number so that it satisfies (j3.8p . I also truncate the F tensors so that 
they keep dimensions 2 x x x X- The normalization condition (|3.10p for F^*) will be still 
satisfied exactly, while the error in the normalization condition ()3.9p will be small. This 
normalization error can be corrected as discussed in the next section. This procedure keeps 
the state within the MPS ansatz with restricted x- 

The procedure described above allows me to efficiently approximately implement local 
unitary evolution. To simulate time evolution 

\m) = e-'"'\m), (3.20) 
with a local Hamiltonian like (j3.ip . I first divide the time t into small slices At and split 

(x) (z) 

the Hamiltonian into two groups of commuting terms H^^ and Hm ■ Each time evolution 
step e~*-^^* can then be implemented as a product of local unitaries using the second order 
Trotter-Suzuki formula 




(3.21) 
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The application of the product of the local unitaries within each group can be done almost 
in parallel (in two steps, as described in Section [3.4p . as they commute with each other. 

These update rules allow me to efficiently approximately simulate the real time evolution 
(j3.20p with a local Hamiltonian H for a state within the MPS ansatz with parameter 
X- The number of parameters in this MPS description with restricted x is then n(2x^) for 
the tensors P^*) and (n — l)x for the vectors A*^*-*. The simulation cost of each local update 
step scales like 0{x^), coming from the SVD decomposition of the matrix Q. For a system 
of n spins, I thus need to store 0{2nx^ + nx) numbers and each update will take O(nx^) 
steps. 

The update procedure generalizes to tree geometry by taking the tensors P with dimen- 
sions 2xxxxxxasin Section [3.2.2[ For a local update (on two neighboring spins i and 
i + 1 with bonds labeled by /, m, n and n, o,p) I rewrite the state \ip) analogously to (|3.16p 

as 

IV') = E E®lia,)M\^c.,)\^a,)\s)At).+^\^a.)\K)- (3-22) 

af^ ,ai ,ao ftip s,t 

using the tensor 

{a^.ai),{aoap) ^^"-k a.1 Z-^ \ ak,ai,am,^^am ^ am,ao,ap J ^^ao ^^ap ) yu.^uj 

with combined indices {atai) and (aoflp). One then needs to update the tensor Q as described 
above (|3.17p - (|3.19p . The decomposition procedure to get the updated vector A*^™) and the 
new tensors P*-*^ and p(*+^) now requires O(x^) computational steps. The cost of a simulation 
on n spins thus scales like 0{nx^). 

3.3.2 Imaginary Time Evolution 

Using the MPS ansatz, I can also use imaginary time evolution with e~^* instead of (j3.20p 
to look for the ground state of systems governed by local Hamiltonians. One needs to 
replace each unitary term e~*^^* in the Trotter expansion (|3.2ip of the time evolution with 
g-AAt followed by a normalization procedure. However, the usual normalization procedure 
for imaginary time evolution (multiplying the state by a number to keep {ip\''P) = 1) is now 
not enough to satisfy the MPS normalization conditions (j3.8p - (l3.10p for the tensors P and 
vectors A I use to describe the state 

The unitarity of the real time evolution automatically implied that the normalization 
conditions (I3.8p . ()3.10p were satisfied after an exact unitary update. While there already 
was an error in (j3.9p introduced by the truncation of the x + 1 • • • 2x entries in P'^*^ the 
non-unitarity of imaginary time evolution update steps introduces further normalization 
errors. It is thus important to properly normalize the state after every application of terms 
like e~^^* to keep it within the MPS ansatz. 

In |106j , Vidal dealt with this problem by taking progressively shorter and shorter steps 
At during the imaginary time evolution. This procedure results in a properly normalized 
state only at the end of the evolution, after the time step decreases to zero (and not 
necessarily during the evolution). I propose a different scheme in which I follow each local 
update e~^^* by a normalization procedure (based on Vidal's observation) to bring the 
state back to the MPS ansatz at all times. The simulation I run (evolution for time t) 
thus consists of many short time step updates e~^^*, each of which is implemented using 
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a Trotter expansion as a product of local updates e~ . Each of these local updates is 
followed by my normalization procedure. 

I now describe the iterative normalization procedure in detail for the case of an infinite 
chain, where it can be applied efficiently, as the description of the state requires only 
two different tensors T (see Section r3.4.ip . One needs to apply the following steps over and 
over, until the normalization conditions are met with chosen accuracy. 

First, for each nearest neighbor pair i,i + 1 with even i, combine the MPS description 
of these two spins (I3.15p - (l3.16p . forming the matrix T ()3.18p . Do a SVD decomposition of 
T (jS.lOp to obtain a new vector A^*). The decomposition does not increase the number of 
nonzero elements of A'-*^ as the rank of the 2x x 2% matrix T ()3.19p was only x (coming 
from (jS.lSp ). I thus take only the first x values of A^*) and rescale the vector to obey 
'^a=i -^a^^ ~ 1- using this new A^*\ I obtain tensors r^*\ r(*+i) from (j3.19p . and truncate 
them to have dimensions x x x x 2. Second, I repeat the previous steps for all nearest 
neighbor pairs of spins i,i + 1 with i odd. 

I observe that repeating the above steps over and over results in exponential decrease 
in the error in the normalization of the T tensors. Note though, that the rate of decrease 
in normalization errors becomes much slower near the phase transition for the transverse 
field Ising model on an infinite line (see Section 13.5. ip . 

In practice, I apply this normalization procedure by using the same subroutine for the 
local updates e~^^^, except that I skip the step (|3.17p . which is equivalent to applying the 
local update with At = 0. The normalization procedure is thus equivalent to evolving the 
state repeatedly with zero time step (composing two tensors T and decomposing them again) 
and imposing the normalization condition on the vectors A. Note though, following from the 
definition of the SVD, that each decomposition assures that one of the conditions (j3.9p . (13.10p 
is retained exactly for the updated tensors T. The errors in the other normalization condition 
for the r tensors are decreased in each iteration step. 

The numerical update rules for a system with tree geometry are a simple analogue of 
the update rules for MPS on spin chains. Every interaction couples two sites, with tensors 
^Ifb 'c ^Ifd 'e ^ with the three lower indices corresponding to the bonds emanating from 
the sites. One only needs to reshape the tensors into r|^^''j'*^ and ^^^^^ and proceed as 
described in p.lSp and below. 

3.4 MPS and Translationally Invariant Systems 
3.4.1 An Infinite Line 

For systems with translational symmetry such as an infinite line all the sites are equivalent. 
I assume that the ground state is translationally invariant, and furthermore pick the tensors 
r^*) and vectors A^*^ to be site independent. For fixed x the number of complex parameters 
in the translationally invariant MPS ansatz on the infinite line scales as 2^^. 

When using imaginary time evolution to look for the ground state of this system, within 
this ansatz, it is technically hard to keep the translational symmetry and the normalization 
conditions after each update. Numerical instabilities plagued my efforts to impose the 
symmetry in the procedures described above. In |106j . Vidal devised a method to deal with 
this problem. Let us break the translational symmetry of the ansatz by labeling the sites 
A and B as in FIG ]3-4i This doubles the number of parameters in the ansatz. The state 
update now proceeds in two steps. Let the site pairs AB interact and update the tensors 
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Figure 3-4: The parametrization and update rules for the infinite hne. 




Figure 3-5: The two-layer, directed labeling of the tree. 



r(^), r^-^) and the vector A^^-^). Then let the neighbor pairs BA interact, after which I 
update the tensors F^'^^ F*^"^^ and the vector A^^"^^. What I observe is that after many state 
updates the elements of the resulting F^"^) and the F^^^ tensors differ at a level which is 
way below my numerical accuracy (governed by the normalization errors) and I am indeed 
obtaining a translationally invariant description of the system. 

One of the systems easily investigated with this method (iTEBD) is the Ising model 
in a transverse field (15.371) on an infinite line. Vidal's numerical results for the real time 
evolution and imaginary time evolution [106\ of this system show remarkable agreement with 
the exact solution. I take a step further and also numerically obtain the critical exponents 
for this system. Further details can be found in Section [3. 5^ where I compare these results 
for the infinite line to the results I obtain for the Ising model in transverse field on the 
Bethe lattice. 



3.4.2 An Infinite Tree 

For the infinite Bethe lattice. My approach is a modification of the above procedure in- 
troduced by Vidal. In order to avoid the numerical instabilities associated with imposing 
site-independent F and A after the update steps, I break the translational symmetry by 
labeling the "layers" of the tree A and B (denoted by half-circles and triangles), as in 
FIGiESl 

The Bethe lattice is also symmetric under the permutation of directions. Tensors F with 
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Figure 3-6: The two-step interactions for the infinite tree. 



fuh directional symmetry obey Ta,b,c = T^^ca = 'rc,a,b = '^c,b,a = '^b,a,c = Ta.cfe- However, for 
the purpose of simple organization of interactions, I will also partially break this symmetry 
by consistently labeling an 'inward' bond for each node, as denoted by the flat sides of the 
semi-circles and the longer edges of the triangles in Fig l3-51 This makes the first of the three 
indices of Ta,b,c special. However, I keep the residual symmetry Ta,b,c = ^a,c,b- I can enforce 
this by interacting a spin with both of the spins from the next layer at the same time. The 
update procedure for the interaction between the spins now splits into two steps, interacting 
the layers in the AB order first, and then in the BA order as in Fig l3-61 Similarly to what 
I discovered for the line, the differences in the elements of the final F^^^^ and F^^^ are well 
below the numerical accuracy of my procedure. 

The scaling of this procedure is more demanding than the O(x^) simulation for a line. 
The number of entries in the matrix used in each update step is 2^^ x 4^^, therefore 
the SVD decomposition requires O(x^) steps. The scaling of my numerical method is thus 
0{x^) for each update step. 

3.4.3 Expectation Values 

A nice property of the MPS state description is that it allows efficient computation of 
expectation values of local operators. First, for a translationally invariant system on a line 
(with only one tensor F and one vector A), I have for an operator O^*) acting only on the 
i-th spin 

X X 

(^io«iv)= o;;;,,x:i:(AaF:;:A,)(A„F:;,A,), (3.24) 

<ij,s'=0,l a=l 6=1 

where O^*^ , = (s^|0^*^|sj). Similarly, for the expectation values of O^^^O^-'^ (assuming j > i), 

Si,s'...,Si,s'. a,,e b,b',... 

x(A„F5A,,F^;+rA,,...A,,Fj;Ae) 
x(AaF:;,A,F^;rA.---A,Fj,Ae). 
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Defining a x matrix B (where one should think of (66') as one combined index ranging 
from 1 to x^) as 

B(bb'),{cc') = y^rg^^rg*^,AcAe', (3.26) 

s 

and vectors v and w with elements again denoted by a combined index (66') = 1 ... as 

= Eol;!.^E(^«)'r:>r:>A,A,,, (3.27) 

-(^^') = Eo!^l.'.E^^e^^:(Ae)^ (3.28) 

I can rewrite p.25p as 

O^^O^^^ IV') = v^ BB - j-B w, (3.29) 



There is a relationship between the eigenvalues of the matrix B and the correlation func- 
tion (O^'^^O^^^) — {0^^^){0^^^). One of the eigenvalues of S is //i = 1, with the corresponding 
right eigenvector 

f^P)=T.T.Vb',ci>^c?, (3.30) 

s c 

and left eigenvector 

Plli] = Sb,b'Xl (3.31) 

which can be verified using the normalization conditions (j3.9p and (j3.10p . I numerically 
observe that /^i = 1 is also the largest eigenvalue. (Note that > 1 would result in 
correlations unphysically growing with distance.) Denote the second largest eigenvalue of 
B as fi2- Using the eigenvectors of B, I can express B^~^~^ in (j3.29p . as 

When computing the correlation function, the term that gets subtracted exactly cancels the 
leading term involving fii = 1. Therefore, if \fi2\ is less than 1, (j3.32p implies 

_ (C)W)(C)(i)) oc /i'^'-^l. (3.33) 

The correlation function necessarily falls of exponentially in this case, and the correlation 
length ^ is related to ^2 as ^ = — l/ln^2- 



The computation of expectation values for a MPS state on a system with a tree geometry 
can be again done efficiently. For single-site operators 0^^\ the formula is an analogue of 
(I3.24P with three A vectors for each T tensors which now have three lower indices. For 
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two-site operators, the terms in ()3.29p now become 



B(cc'),{dd') - y^y^rgerfrg*e^^,(Ae)^ArfArf/, (3.34) 



s e 



V(cc') 



E E(^'^)'(^^)'r:>,cr3,c'AcAe', (3.35) 



, a.b 



^i'id') = E0Z'E^2e,f^J,eji^ef{Xf?. (3.36) 



The correlation length is again related to the second eigenvalue of the B matrix as in (|3.33p . 

3.5 Quantum Transverse Field Ising Model 

My goal is to investigate the phase transition for the Ising model in transverse magnetic field 
p.ip on the infinite line and on the Bethe lattice. I choose to parametrize the Hamiltonian 
as 

^ = I E (1 - + E (1 - <) ' (3-37) 

where < s < 1 and b is the number of bonds for each site (5 = 2 for the line, 6 = 3 for 
the tree). I will investigate the ground state properties of p.37p as we vary s. The point 
s = corresponds to a spin system in transverse magnetic field, while s = 1 corresponds to 
a purely ferromagnetic interaction between the spins. 

3.5.1 The Infinite Line. 

I present the results for the case of an infinite line and compare them to exact results 
obtained via fermionization (see e.g. [82], Ch.4). Vidal has shown |106j that imaginary time 
evolution within the MPS ansatz is capable of providing a very accurate approximation for 
the ground state energy and correlation function. I show that even using x smaller than used 
in [106j . I obtain the essential information about the nature of the phase transition in the 
infinite one-dimensional system. I also obtain the critical exponents for the magnetization 
and the correlation length. 

In FIG 13-71 I show the how the ground state energy obtained using imaginary time 
evolution with MPS converges to the exact energy as x increases. The exact solution for a 
line has a second order phase transition at the critical value of s, sl = |, and the ground 
state energy and its first derivative are continuous, while the second derivative diverges at 
s = Sl- I plot the first and second derivative of E with respect to s obtained numerically 
and compare them to the exact values in FIG. \'6-8\ observing the expected behavior already 
for low X- 

The derivative of the exact magnetization Af = (c^) is discontinuous at sl, with the 
magnetization starting to rise steeply from zero as 

Moc(x-xl)^, (3.38) 

with the critical exponent Pl = h- Here x is the ratio of the ferromagnetic interaction 
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Figure 3-7: Transverse Ising model on an infinite line. Fractional difference of the ground 
state energy obtained using MPS and the exact ground state, near the phase transition at 
sl = ^- The energy scale is logarithmic. 
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Figure 3-8: Transverse Ising model on an infinite line. The first and second derivative with 
respect to s of the ground state energy obtained via MPS compared with the exact result. 
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Figure 3-9: Transverse Ising model on an infinite line. Magnetization obtained using MPS 
vs s, with = 10"^. 



strength to the transverse field strength in ()3.37p and is 

X = — ^ — -, (3.39) 
2(1 -s)' ^ ^ 

with the value x = = 1 at the phase transition (s^, = |). I plot the magnetization 
obtained with my method in FIG I3-91 To obtain the magnetization depicted in the plot, I 
used a small symmetry breaking longitudinal field with magnitude = 10~®. In FIG I3-10| 

1 plot M vs. X — XL on a log-log scale. I also plot a line with slope 0.125. Observe that as x 
increases, the data is better represented by a line down to smaller values of x — xl. For the 
largest x displayed, a straight line fit of the data between 4 x 10~^ < x — xl < 10~^ gives me 
a slope of 0.120. Note that the mean field value of the critical exponent for magnetization 
is 0.5. 

The correlation function (cri^Vi"''') — {(J^^){az''^) can be computed efficiently using (j3.25p . 
Away from criticality, it falls off exponentially as e"'*"-^'/''. The falloff of the correlation 
function is necessarily exponential as long as fi2, the second eigenvalue of B, is less than 1. 
The exact solution for the correlation length ^ near the critical point has the form 

^(x\x-xl\~^ (3.40) 

as X approaches xl. Already at low x the iTEBD method captures the divergence of the 
correlation length. In FIG. 13-lH I plot ^ vs. \x — xl\ on a log-log plot, together with a 
line with slope —1. Again, as x increases, the data is better represented by a line closer 
to the phase transition. For the highest x displayed, a straight line fit of the data between 

2 X 10~^ < XL — x < 4 X 10~^ gives me a slope of —0.92. Note that the mean-field value 
of the critical exponent for the correlation length is 0.5 (corresponding to slope —0.5 in the 
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Figure 3-10: Transverse Ising model on an infinite line. Log-log plot of magnetization vs. 
X — XL. I also plot a line with slope = 0.125. 
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Figure 3-11: Transverse Ising model on an infinite line. Log-log plot of the correlation 
length vs. \x — xl\. I also plot a line with slope —i'l = —1- For each x in the plot, I choose 
a numerical value of the critical point as the point at which the correlation length is 
maximal. 
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Figure 3-12: Transverse Ising model on an infinite tree. The first and second derivative 
with respect to s of the ground state energy obtained via MPS. 



graph) . 

3.5.2 The Infinite Tree (Bethe Lattice). 

The computational cost of the tree simulation is more expensive with growing x than the 
line simulation, so I give my results only up to x = 8. I run the imaginary time evolution 
with 10000 iterations (each iteration followed by several normalization steps) for each point 
s, taking a lower x result as the starting point for the procedure. I also add a small 
symmetry-breaking longitudinal field with magnitude = W~^. 

I see that the energy and its first derivative with respect to s are continuous. However, 
I now observe a finite discontinuity in the second derivative of the ground state energy 
(see FIG. I3-I2p . as opposed to the divergence on the infinite line. This happens near 
s = ST^ 0.5733. 

Similarly to the one-dimensional case, the magnetization quickly grows for s > st, while 
it has (nearly) zero value for s < st (see FIG. I3-I3p . In FIG I3-1'I1 I plot the magnetization 
vs. X — XT on a log-log scale, where x is 

(3.41) 



3(1 -s)' 

with the value x = xy ~ 0.451 at the phase transition (where s = st ^ 0.5733). I want to 
test whether the magnetization behaves like 

M oc (x - xrf (3.42) 

for X close to xt, which would appear as a line on the log-log plot. As x grows, the data is 
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Figure 3-13: Transverse Ising model on an infinite tree. Magnetization vs. s, with = 
10-8. 



better represented by a straight line closer to the phase transition. If I fit the x = 8 data 
for 4 X 10"^ < X — xl4x < lO^'^, I get /3 = 0.41. I add a line with this slope to the plot. 
Note that the mean-field value for the exponent f3 is 0.5, just as it is for the infinite line. 

I observe that the correlation length now rises up only to a finite value (see FIG. I3-I5p . 
As I increase Xj the second eigenvalue of the B matrix, ^2, approaches a maximum value 
close to ^. I conjecture that the limiting value of ^2 is indeed ^, which corresponds to 
a finite correlation length with value (In 2)"^. Note that for the infinite line, the second 
eigenvalue of B approaches 1, and so the correlation length is seen to diverge at the phase 
transition. 



I now look at a model with a different interaction term. Starting with an antiferromagnetic 
interaction, I add a specific longitudinal field at each site. As in the previous section, I 
parametrize the Hamiltonian p.2p with a single parameter s: 



with 6 = 2 on the line and 6 = 3 on the tree. I choose the longitudinal field in such a 
way that the nearest-neighbor interaction term Hij becomes a projector, expressed in the 



3.6 The Not 00 Model 



NOT 00 




(3.43) 
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Figure 3-14: Transverse Ising model on an infinite tree. Log-log plot of magnetization vs. 
X — XT, with Bz = 10~®. I also plot a line with slope 0.41. 
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Figure 3-15: Transverse Ising model on an infinite tree. A linear plot of the correlation 
length vs. s. 
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Figure 3-16: The not 00 model on an infinite line. The ground state energy and its first 
two derivatives with respect to s. 



computational basis as 

Hij = m{00\^^, (3.44) 

thus penalizing only the |00) configuration of neighboring spins. Accordingly, I call this 
model NOT 00. The ground state of the transverse Ising model (I3.37P at s = 1 has degeneracy 
2. For (j3.43p on the infinite line or the Bethe lattice, the degeneracy of the ground state at 
s = 1 is infinite, as any state that does not have two neighboring spins in state |0) has zero 
energy. 

3.6.1 Infinite Line 

I use my numerics to investigate the properties of (|3.43p on the infinite line as a function 
of s. My numerical results show continuous first and second derivatives of the energy 
with respect to s (see FIG J3-16"]) . The magnetization M = (az) decreases continuously and 
monotonically from at s = to a final value of —0.606 at s = 1 (see FIG I3-171) . The 
second eigenvalue of the B matrix (j3.33p rises continuously from at s = 0, approaching 
0.603 at s = 1 (see FIG J3-17p . Because /i2 < 1, the correlation length ^ is finite for all 
values of s in this case. These results imply that there is no phase transition for this model 
as I vary s. 

As a test of my results, I compute the magnetization at s = 1 exactly for this model 
on a finite chain (and ring) of up to n = 17 spins. I maximize the expectation value of 
Hb = Yli within the subspace of all allowed states at s = 1 (with no two zeros on 
neighboring spins), thus minimizing the expectation value of the second term in (j3.43p for 
s approaching 1. I compute the magnetization M = (al) for the middle i = [^J spin 
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Figure 3-17: The NOT 00 model on an infinite line. Magnetization as a function of s and 
correlation length cLS Si function of s. 



for the ground state of the not 00 model exactly for a finite chain and ring of up to 17 
spins at s = 1. As I increase n, the value of M converges to —0.603 (much faster for the 
ring, as the values of M for n = 14,17 differ by less than 10^"^). Recall that I obtained 
M = —0.606 from my MPS numerics with x = 16 for the NOT 00 model on an infinite 
line. I also compare the values of the Schmidt coefficients across the central division of the 
finite chain (n = 16) to the elements of the A vector obtained using my MPS numerics with 
X = 32. I observe very good agreement for the 11 largest values of A^, with the difference 
that my MPS values keep decreasing (exponentially), while the finite-chain values flatten 
out at around Xk>i4 ~ 10~^ (see FIG I3-18]) . The behavior of the components of A from 
MPS doesn't change with increasing x- 

In the ground state of (|3.43p at s = 1, the overlap with the |00) state of any two 
neighboring spins is exactly 0. If the T tensors are the same at every site, the component of 
the state \^) that has overlap with the state |00) on nearest neighbors can be expressed as 

^(A„rO,A6rO,A,) |OO)|0e). (3.45) 

a,b,c 

Furthermore, when the F tensors are symmetric, the elements of the A vectors must be 
allowed to take negative values to make this expression equal to zero. Note that until now, 
I used only positive A vectors, knowing that they come from Schmidt decompositions, which 
give me the freedom to choose the components of A to be positive and decreasing. 

The negative signs in the A vector can be absorbed into every other F tensor, resulting 
in a state with two different F (for the even and odd-numbered sites) and only positive 
A's. In fact, this is what I observe in my numerics, which assume positive A, but allow 
two different F tensors (see 13. 4.1]) . If I allow the elements of A to take negative values, my 
numerically obtained F tensors are identical. 
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Figure 3-18: Ground state of the NOT 00 model on a line at s = 1. Comparison of the exact 
Schmidt coefficients for a division across the middle of a finite chain and of the MPS values 
(x = 32) for an infinite hne. 



3.6.2 Infinite Tree 

Here, I numericaUy investigate the not 00 model (j3.43p on the Bethe lattice. As on the 
line, the numerics show continuous first and second derivatives of the energy with respect 
to s (see FIG J3-19]) and a continuous decrease in the magnetization from at s = to 
—0.671 at s = 1 (see FIG I3-20|) . The correlation length behaves similarly as on the line, 
increasing with s, but it reaches a maximum at s = 0.96 for x = 8. The maximum value of 
^ is apparently lower than l/ln2, (see FIG J3-20]) . meaning that on the tree, the correlation 
function {a^z^ai^^) — {(7't^){cr^P) falls off with distance faster than 2~l*~-'l for all s. 



3.7 Stability and Correlation Lengths on the Bethe Lattice 

I have found that on the Bethe lattice, for both my models, the second eigenvalue 112 of 
the matrix B ()3.33p . which determines the correlation length, apparently is never greater 
than i. In this section I argue that this is a model-independent, and calculation method 
independent, consequence of assuming that a translation-invariant ground state is the stable 
limit of a sequence of ground states of finite Cayley trees as the size of the tree grows. For 
a related problem, the stability of recursions for Valence Bond States on Cayley trees has 
been investigated by Fannes et.al. in [30j . 

The Hamiltonians (j3.ip and ()3.2p each consist of sums of terms , Hm"^ , as in (j3.2ip , 

where each term depends on a single ax and each on a neighboring pair oi (Tz- I 
calculate the quantum partition function 

Z(/3) = tre"'^^ (3.46) 
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Figure 3-19: The not 00 model on an infinite tree. The ground state energy and its first 
two derivatives with respect to s. 
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Figure 3-20: The not 00 model on an infinite tree. Magnetization as a function of s and 
correlation length cLS Si function of s. 
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Figure 3-21: A system of classical spins on a lattice whose layers are Cayley trees. A and 
B denote the Boltzmann factors. 



as the limit of 



Z{N, At) = tr 



TV 



(3.47) 



as At — > 0, ^ oo with NAt = [i. To find the properties of the ground state, I take 
/3 ^ oo so that I need Z{N, At) as At ^ 0, iV ^ oo with iVAt oo and N{At)^ (to 
make the error in using the Trotter-Suzuki formula go to zero). I interpret ()3.47l) as giving 
the classical partition function of a system of Ising spins (s = ±1) on a lattice consisting 
of N horizontal layers, each of which is a Cayley tree of radius M (i.e with a central node 
and concentric rings of 3, 3 x 2, 3 x 2^, . . . , 3 x 2^~^ nodes). I can write ()3.47p as 

Z{N,At) = Y,\{A\{B, (3.48) 

where the sum is over all configurations of A^x (3x2^ — 2) spins s = ±1 and the products are 
of a Boltzmann factor A for each horizontal link in the Cayley trees, and a Boltzmann factor 
B for each vertical link between corresponding nodes in neighboring layers (see FIG I3-21]) 
(layer N is linked to layer 1 to give the trace). The factor A for the link between nodes i,j 
in the same horizontal layer is given by 

A{si,Sj) = e"^*-^W)(""^^\ (3.49) 
The factor B for the link between nodes i' in the same vertical column is given by 

B{sus[) = {a, = s'i\e-^'"^) \a, = Si). (3.50) 

(z) (x) 

When all the terms -f^(jj) are of the same form, as are all the terms H^^-^ , the form of the 
factors A and B does not depend on which particular links they belong to. 

Each term in the sum, divided by Z, can be thought of as the probability of a config- 
uration {s}. In what follows I will keep A^ and At fixed and consider the limit M — > oo, 
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Figure 3-22: The M = 2 Cayley tree. 



i.e. finite Cayley tree — > Bethe lattice. I will then suppose that my results, which are 
independent of the form of A and B (provided A,B>{)) will also hold after the — > oo 
limit is taken, i.e. for the quantum ground state. 

I will think of the lattice as a single tree, with each node being a vertical column of 
spins. (A recent use of this technique to investigate the quantum spin glass on the Bethe 
lattice is in [59].) Denote by s the vector of N values of s along a column. Let K{s) be 
the product of the N factors B{s,s') along a column and let L{s,s') be the product of 
N factors A(s, s') on the horizontal links between nearest neighbor columns. Let Zm be 
the partition function for a tree of radius M. I can calculate Zm by a recursion on M as 
follows: 

Zm = Y.K{s)[FM(s)f, (3.51) 

s 

FMis) = Y.L{s,s')Kis')[FM-i{s')f, (3.52) 

s' 

Fo{s) = 1. (3.53) 

It is easy to see that this recursion gives the correct Zm (the case M = 2 is shown in 
FIG|332|). I also see that 

PMis) = -^Kis)[FM{s)f (3.54) 

is the probability of the configuration s along the central column. For the case = 1, 
i.e. classical statistical mechanics on a tree, this is the well-known method to find an exact 
solution [H]. 

In order to have a well-defined translationally invariant limit as M ^ oo I would like 
the recursion (j3.52p for Fm to have an attractive fixed point F which Fm approaches as 
M — > oo. 'Attractive' means that if I start the recursion with a different Fo(s), sufficiently 
close to Fq{s) = 1, the limiting value of FMis) will be the same fixed point. This in turn 
implies that on a Cayley tree with large M, small changes in the Hamiltonian on the outer 
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edge will have small effects on the properties of the central region. 

First however I need to fix the overall normalization of Fm{s), since if Fm{s) satisfies 
(I3.52p . so does a^'^ Fm{s) which rules out an attractive fixed point. 

Let 

Fm{s) = ZIFm{s), (3.55) 

so that 

Y,Kis)[FM{s)f = 1, (3.56) 

s 

and 

Pm{s) = K{s)[FM{s)f. (3.57) 

The recursion relation becomes 

Fm{s) = XMY,L{s,s')K{s')[FM-i{s')f, (3.58) 

s' 

with Xm determined by the normalization condition (j3.56p . I can now suppose that 

Fm{s) F{s) as M ^ oo, (3.59) 

with 

F{s) = XY,Lis,s')K{s')[Fis')f, (3.60) 

s' 

and 

Y,K{s)[F{s)f = 1. (3.61) 

s 

To determine whether F is an attractive fixed point, let 

Fm{s) = F{s) + fM{s), (3.62) 

Am = A(l + eAf), (3.63) 

with Jm and em ^ as M — > oo. To first order in Jm, ^m, (|3.58p and (|3.56p become 

fM{s) = eMF{s) + 2Y,T{s,s')fM-i{s'), (3.64) 

s' 

Y,K{s)[Fis)]'fM{s) = 0, (3.65) 

s 

where 

T{s, s') = XL{s, s')K{s')F{s'). (3.66) 
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From l^m\i . 

Y,ns,s')F{s') = F{s), (3.67) 

s' 

and since L{s, s') = L{s', s), 

Y,K{s')[Fis')fT{s',s) = K{s)[F{s)?, (3.68) 

i.e. the linear operator T has an eigenvalue one, with right eigenvector F and left eigenvector 
KF^ (which from (j3.6ip have scalar product one). (j3.64p now gives 

Y,K{s)[F{s)ffM{s) = eM + 2Y,K{s)[F{s)ffM-i{s), (3.69) 

s s 

so from (I3.65p . €m = 0. Let 

T^{s,s') = Tis,s')-Hs)Kis')[Fis')?, (3.70) 

so that 

J2T^is,s')F{s') = 0, (3.71) 



and 



(j3.64p now becomes 



Y,K{s')[F{s')?T^{s',s) = 0. (3.72) 



fM{s) = 2 ^ T^is, s')fM-iis'). (3.73) 



(j3.73p shows that F(s) is an attractive fixed point if and only if 



|r^||<i (3.74) 



(for a tree with valence p + 1 at each vertex, ^ is replaced by ^). 



I can in fact prove that there does exist an F{s) satisfying (j3.60p and (j3.6ip . for which 
the corresponding T"*- has a maximum eigenvalue less than ^. Define a function <I> of F{s) 

by 



m = y\F{s)fK{s)L{s,s')Kis')[Fin?- (3.75) 
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Figure 3-23: Computing the probability distribution of sq and Sd- 



I look for a maximum of <I> with F{s) restricted to the region 

Y,K{s)[F{s)f = 1, (3.76) 

s 

F{s) > 0. (3.77) 

A maximum must exist, but it might be on the boundary of the region, i.e. it might have 
F{s) = for some values of s. Elementary calculations (omitted here) establish that 
stationary values of $ on the boundary cannot be maxima. At stationary points in the 
interior of the region, i.e. with F{s) > for all s, ()3.60p and (I3.6ip must be satisfied. 
If such a stationary point is a maximum, all the eigenvalues of 1 — 2T-^ are > 0, i.e. all 
the eigenvalues of T"*" are < ^. This is weaker than the attractive fixed point condition, 
which also requires that no eigenvalue is less than — i, but does correspond to the observed 
property of ^2- 

I now examine the joint probability distribution of sq and s^, where denotes the central 
column and d a column distance d from the center. One can see from FIG 13-23] that 

PMiso,Sd) = ^ V [FM{so)fK{so)Liso,si)FM^iisi)... (3.78) 

si,— ,Sd-l 

. . . L{sd-i, Sd)K{sd)[FM-d{sd)]'^ . 
As M — > oo (with d fixed) this becomes (using (|3.66p ) 

P{so,Sd) = C[F{so)]^K{so)T''{so,Sd)Fisd), (3.79) 
where the normalization C is determined by 

Y,Piso,Sd) = l. (3.80) 

Using ([3:671) I find 

P{so) = J2 P(^o, Sd) = CK{so)[F{so)f, (3.81) 

Sd 

so from (|3.6ip . C = 1 (and P{so) agrees with the limit of (|3.57p ). Expressing ()3.79p in 
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terms of T-^ using dSTTO]) . (I3TT|) and (l3T2]l . 



P(so,Sd)-P(so)P(sd) = K(sl,)[F(so)]2(r^)'^(so,Srf)F(sd). (3.82) 

Thus the correlation between sq and falls off as fi'^, where /i is the eigenvalue of T"*" with 
maximum modulus, and so from (j3.74p . faster than 1/2^^. 

If this conclusion is correct (and clearly the argument is less than rigorous), it estab- 
lishes more than my experimental observation that fi2 < 5- The quantum limit of -P(so, s^) 
encodes not only the static correlation {ipol SQSd iV'o) in the ground state, but also the imag- 
inary time dependent correlation {ipo] e^*soe~^*Sd lipo) which in turn determines the linear 
response as measured by Sd to a time-dependent perturbation proportional to sq. If this 
indeed falls off faster than 1/2*^, then there is some hope that the Bethe lattice can be used 
as a starting point for investigation of fixed valence random lattices. 
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Chapter 4 



Quantum 



Satisfiability 



The topic of this chapter is Quantum Satisfiability, the quantum analog of classical Satisfi- 
ability. I discussed the locally constrained classical problem in Section 11.1.31 where I also 
summarized the well known complexity of its many variants (see Table [LT] in Section [LL2]) • 
On the other hand, there are still several open questions about Quantum Satifiability, and 
I choose to investigate some of them here. 

First, in Section r4.1.1I I introduce Quantum Satisfiability and summarize what is presently 
known about its complexity. Then I review the proof techniques I use throughout this chap- 
ter, focusing on several clock constructions in Section [4.1.3[ Section is based on the paper 



A New Construction for a QMA Complete 3-local Hamiltonian 

Daniel Nagaj, Shay Mozes 

We present a new way of encoding a quantum computation into a 3-local 
Hamiltonian. Our construction is novel in that it does not include any 
terms that induce legal-illegal clock transitions. Therefore, the weights of 
the terms in the Hamiltonian do not scale with the size of the problem as 
in previous constructions. This improves the construction by Kempe and 
Regev [5T , who were the first to prove that 3-local Hamiltonian is complete 
for the complexity class QMA, the quantum analogue of NP. 
Quantum k-SAT, a restricted version of the local Hamiltonian problem 
using only projector terms, was introduced by Bravyi 20J as an analogue 
of the classical k-SAT problem. Bravyi proved that quantum 4-SAT is 
complete for the class QMA with one-sided error (QMAi) and that quantum 
2-SAT is in P. We give an encoding of a quantum circuit into a quantum 
4-SAT Hamiltonian using only 3-local terms. As an intermediate step to 
this 3-local construction, we show that quantum 3-SAT for particles with 
dimensions 3 x 2 x 2 (a qutrit and two qubits) is QMAi complete. The 
complexity of quantum 3-SAT with qubits remains an open question. 



One of the implications of this construction is that Adiabatic Quantum Computation can 
be used for simulating quantum circuits much more effectively than Aharonov et al. showed 
in [7] (see Section [2?2]) . 

Next, in Section [4.31 I review the result of Aharonov et al. [5] who proved that Quantum 
2-SAT on a line with 12-dimensional particles is QMAi complete. I decrease the required 
dimensionality of particles in their construction to d = 11 in Section [4.3.21 



m 
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In Section 14.41 I present a very recent result by Eldar and Regev, who proved that 
Quantum 2-SAT for a cinquit-qutrit pair (particles with d = 5 and d = 3) is QMAi 
complete, using a novel triangle clock construction. In Section 14.4.21 I use their result to 
obtain another QMA-complete 3-local Hamiltonian construction built from constant norm, 
3-local terms of a restricted type. 

Finally, I take yet another step towards figuring out the complexity of Quantum 3-SAT in 
Section [4.51 Using my train switch clock construction, I show that a simple Quantum 3-SAT 
Hamiltonian is universal for quantum computation in the Adiabatic Quantum Computing 
model of Chapter [2] and in the more general Hamiltonian Computer model of Section 12.2.11 
Moreover, the required running time for this model scales only slightly worse than linearly 
with the length of the computation. 

4.1 Introduction 
4.1.1 Quantum k-SAT 

The Quantum A;-SAT (Q-A:-SAT) promise problem was introduced by Bravyi [20] as an 
analogue of classical /c-SAT. In classical Satisfiability, one needs to determine whether there 
exists a bit string satisfying all the boolean clauses of a given problem instance. In Q-k- 
SAT, the problem is to determine whether a Hamiltonian acting on n qubits has a zero 
eigenvalue, or whether all its eigenvalues are higher than e > n~° for some constant a. 
Moreover, the Hamiltonian 

^Q-fe-SAT = J^Pi, (4.1) 

is composed of /c-local projector terms 

Pi = Pl=I^'-''-''>^mm{,.,,,,.^y, (4.2) 

where each Pj acts nontrivially on k qubits {q\...ql.}. For a 'yes' answer to a problem 
instance, the ground state of -f^Q-fc-SAT rnust be exactly zero, i.e. a state which is annihilated 
by all the projectors Pi must exist. This shows the analogy to classical SAT, where all the 
boolean clauses have to be satisfied. If the states are computational basis states and 
e = 1, the problem reduces to classical A:-SAT. The special case of commuting projectors Pi 
has been analyzed in [21J. 

Local Hamiltonian, the problem I introduced in Section ll.3.2| is more general than 
Q-A;-SAT. It is the quantum analogue of classical MAX-A:-SAT, where one is interested 
in the properties of the ground state of the total Hamiltonian (the sum of the terms). 
Viewing Quantum A;-SAT and fc-local Hamiltonian as quantum analogues of classical fc-SAT 
and MAX-fc-SAT, it is interesting to investigate and compare their complexities. Known 
complexity results about the classical problems are summarized in Section 11.1.21 Table 11.11 
In Table 14.11 I present the currently known results about Quantum Satisfiability and Local 
Hamiltonian problems. In this Chapter, I prove those marked with *. 

In |2Qj, Bravyi gave a classical polynomial algorithm for Quantum 2-SAT showing that 
Q-2-SAT belongs to P. He then proved that Quantum fc-SAT is QMAi complete for k > A. 
The class QMAi is a special case of QMA, where only one sided error is allowed (see 
Definition [8] in Section ri.2.2p . When the answer to the problem instance is 'yes', the verifier 
circuit U can output 'yes' on some state with certainty. Thus, the ground state energy 
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Quantum 


qubits 


qudits 


Q-A;-SAT 


/fc = 2 : in P 

A; = 3 : contains NP, 

universal for BQP* 
A; > 4 : QMAi-complete 


Q-(5,3)-SAT : QMAi-complete 
Q-(ll, 11)-SAT in ID : 

QMAi- complete* 

Q-(3,2,2)-SAT : QMAi-complete* 


fc-local 
Hamiltonian 


k>2: QMA-complete 

even on a 2D grid 
A; = 3 with constant norm terms : 

QMA-complete* 





Table 4.1: Known complexity of Quantum Satisfiability and Local Hamiltonian problems. 
My own results, presented later in this Chapter, are marked with *. 



of a Quantum A;-SAT Hamiltonian must be exactly zero in the 'yes' case. When I release 
the certainty requirement on the verifier circuit C/, it becomes a verifier circuit for a QMA 
problem instead. Therefore, when I prove that a certain Q-SAT problem is QMAi complete, 
the QMA-completeness of the mother Local-Hamiltonian problem immediately follows as 
a side result. I present a result of this type in Section 14.21 where I give a new 3- local 
Hamiltonian construction with constant norm terms. 

Bravyi's original definition also required all of the terms in the Hamiltonian to be pro- 
jectors. However, using A;-local positive semidefinite operator terms Hf with zero ground 
state and constant norm instead of projectors Pi in (j4.ip is an equivalent problem. Quan- 
tum A:-SAT with positive semidefinite operators contains Q-/c-SAT with projectors. On the 
other hand, if one is able to solve Q-/c-SAT with projectors, one can solve Q-A:-SAT with 
positive semidefinite operators as well. For each positive semidefinite operator Hf, define 
a projector P-^ with the same ground state subspace. If Hp = ^ has a zero ground 
state, so does ff+ = "^Hf. If the ground state energy of Hp = "^P^ is greater than e, 
the ground state energy of Hj^ is greater than ce, where c is a constant. Therefore, instead 
of only projector terms, I can equivalently use positive semidefinite operators in Q-fe-SAT 
Hamiltonians. 

4.1.2 Constructing the Hamiltonians 

The two types of results about Quantum k-SAT Hamiltonians that I present here are proofs 
of QMA (or QMAi) completeness and proofs of universality of using a Hamiltonian for 
simulating quantum circuits. In both cases, the quest is to prove these results for ever 
simpler Hamiltonians. Better locality (lower k, the number of particles involved in each 
term), simpler geometry of interactions (2D grid, ID), lower dimensionality of particles 
involved, and using only terms of a restricted type are what I am aiming at. 

To prove that a Hamiltonian problem is QMA-complete, one first needs to show that it 
belongs to QMA, and then to show that it also contains QMA. I have presented Kitaev's 
proof of QMA-completeness for the 5-Local Hamiltonian problem in Section 11.3.21 and 
because a lot of the results I present in the next two Chapters build on it, I now quickly 
summarize it. Before considering the special case of 5-local Hamiltonian, Kitaev proved 
that /c-local Hamiltonian is in QMA for any constant k. Therefore, the question whether 
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a Hamiltonian problem belongs to QMA is usually answered quickly by a reduction to 
fc-Local Hamiltonian. For the other direction, one needs to reduce every instance of a 
problem in QMA into an instance of the Hamiltonian problem in question. Every instance 
of a problem in QMA is a 'yes/no' question. Moreover, every problem instance has a 
verifier circuit U which verifies claims (supposed proofs) about the answer. This circuit 
must have the completeness and soundness properties. Completeness means that if the 
answer to the problem instance is 'y^s', there exists a state that the circuit U will accept 
with high probability. The soundness property means that if the answer to the problem is 
'no', there is no state that fools the verifier circuit U with high probability. I then need to 
construct a Hamiltonian whose ground state would have low energy (or exactly zero energy 
for Quantum fc-SAT) if the answer is 'yes', and whose ground state would have considerably 
higher energy if the answer is 'no'. 

All the results in this Chapter use the same general form of this Hamiltonian, but differ 
in the implementation of some of the terms. Let me first look at what they have in common 
and focus on the differences in Section [4.1.31 The verifier quantum circuit U with L gates 
acts on n qubits. Consider now a quantum system whose Hilbert space 7i is larger than 
T~^work = C®". A way to do this is to append a clock register to the system as 

'H = Hwork 'Hdock- (4.3) 

Then, encode the progression of the circuit U 

\i>t) = UtUt^i . . . U2U1 \^o) , t = 0, . . . , L (4.4) 

into states 



m = \^t)(S)\t), (4.5) 

of the larger system. How exactly the states \t)^ and transitions between them are imple- 
mented constitutes a 'clock construction' and I explain the possibilities in detail in Section 
14.1.31 Following Kitaev, the Hamiltonian I then construct for the larger system is 

H — Hprop ~\~ HquI + Hinput ~l~ -f^ciocfc ~l~ Hclockinit- (^'6) 

It is built in such a way that the history state for the quantum circuit U corresponding to 
some initial state {tpo) 

1 ^ 1 ^ 

\^h^story) = -J=== ^ l^t) = -7=^ 1^*) ® (4-7) 

+ ^ t=o + 1 

is close to the ground state of H . In fact, for a 'yes' instance of Quantum A;-SAT, it is the 
ground state of H. 

The first term, Hprop is a Hamiltonian whose ground state has the form {^history)'- 

L 

Hprop = Hpj.gp (4-8) 

t=i 

H'prop = \{l®{Pt + Pt-i)-Ut®Xt,t-i-{Ut®Xt,t-i)^) , (4.9) 
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where Pt = \t) {t\^ is a projector onto the state \t)^ of the clock register, the operator 
Xt,t-i = 1^) ~ l|c increases the time register from t — 1 to t and the unitary gate Ut acts 
on the corresponding work qubits. Observe that it is built from terms which check whether 
the amplitude of each l^'f-i) and {"^t) is the same. Therefore, the history state (I4.7p is its 
eigenvector with eigenvalue 0. 

The last of the n work qubits, Qn, is the output qubit for the circuit U. The term 

Hout = \0){0\g^^PL (4.10) 

takes care of whether the circuit U accepts some state or not. It gives an energy penalty 
to any history state built from a state IV'no) that is not accepted by U, i.e. whose output 
qubit after L gates is not in the state |1). The term 

Hinput= Il)(l|,,®^0 (4.11) 

k&ancilla 

takes care of proper initialization of ancilla qubits, in the case they are required for U. The 
final two terms Hdock and Hdockinit add an energy penalty to all states of the clock register 
Ticiock that are not in the subspace spanned by the states \t)^. The implementation of these 
two terms and how the states \t)^ are actually encoded is the topic of the next Section. 

For the other type of result, universality of a Hamiltonian H for simulating a quantum 
circuit, one needs to show that time evolution (for not too long a time) from a simple initial 
state with the Hamiltonian H can produce the output of the circuit U with high probability. 
This type of proof again requires the encoding of the progression of a quantum circuit U 
into states \^t) of a larger system. However, proving this type of result is different in that 
I need to analyze the dynamics of the system, not just the ground state of H. However, 
what simplifies the problem is that I can choose the initial state myself. I present one such 
result in Section 14.51 Later, the whole Chapter [5] is dedicated to showing universality for 
Hamiltonians which are both time and translation-invariant. 

4.1.3 Constructing Clocks 

In this Section, I review several methods of implementing the clock register in (14. 3p . The 
goal of a clock construction is to encode the progression ()4.4p of a quantum circuit U with 
L gates on n work qubits into a sequence of orthogonal states \^t) of a larger system Ti. It 
is desirable for a clock construction to have: 

• Easily implementable unique transitions 

One needs to construct a Hamiltonian inducing transitions between the states |^t)- It is 
important to ensure that each state \^t) can transition directly only to the states or 
i-e. 

H\^t) = at-i\^t-i) + (3t\^t) +lt+i\^t+i) , (4.12) 

and not into any other state. Undesirable transitions can be suppressed by adding large 
penalty terms as in the work of Kempe and Regev [5l], but such constructions can not be 
used for proving results about Quantum fe-SAT. In Sections I4.2M.5I I give ways to construct 
the clock register with unique transition rules, allowing me to prove new results about 
Q-fe-SAT and strengthen old results about the Local Hamiltonian problem. 
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• Locally checkable encoding, initialization and final detection 

The full Hilbert space TL can be much larger than the subspace Tiiegai spanned by the states 
\^t)- I niust be able to detect whether a state belongs to the subspace TLiegai using only local 
terms. In practice, this is done by adding an energy penalty to any state outside TLiegai- 
Also, I need to be able to make a local projection measurement onto \^ l) and 1^'^)"'". 

• Local transition rules involving only a few (clock) particles, 

• Low dimensionality of the (clock) particles, 

Simple local interactions in the system 7i and low dimensionality of the particles involved 
are essential for my goal of obtaining stronger results about ever simpler local Hamiltonians. 

• Simple geometry of interactions. 

Finally, for practical implementation, the /c-particle interactions in H need to involve only 
spatially close particles. The ID construction of [5] presented in Section [4.31 is one example 
of this. 

Let me now present a few of clock register encodings and discuss their advantages and 
disadvantages. 

Adding a clock register 

The first option is to directly add a clock register to the system as in Feynman's construction 
in Section 11.3.11 

H = Hwork ® 'Hdock- (4-13) 

A set of orthogonal states \t)^ of the clock register is then used to label the progress of the 
computation, while the contents of the work register holds the state of the work qubits after 
t gates have been applied to an initial state {tpo) 

\^t) = {UtUt-i...Ui\il,o))®\t)^. (4.14) 

^ V ' 

The space spanned by the clock states \t)^ then defines the legal clock subspace TLiegai- 
An essential ingredient of the clock constructions I present here is a local way of checking 
whether a state of the clock register belongs to the legal clock subspace. It is usually 
implemented by projector terms adding an energy penalty to states outside of TLiegai- 

One could think of implementing the clock register by simply labeling the basis states 
of TLciock by \t)^. The minimum required number of qubits in TLciock would then be logL. 
However, such a construction is not local, as the transition rules between different clock 
states necessarily involve (logL) qubits. 

The simplest local implementation of the clock register is the pulse clock used by Feyn- 
man [38]. A state corresponding to time t is the state of L + 1 qubits cq, . . . , cl, with a 
single clock qubit cj in the state |1) 

= |00... 00100... 00) . (4.15) 
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This single up-spin denotes the active site in the clock register. To transition between clock 
states \t)^ and \t + 1)^ (and vice versa), I need to look only at 2 clock qubits and apply the 
operator 

Xt+,,t + xl,^, = |01) (10|,^^,^^^ + |10) {0%^^^^^ . (4.16) 

This low locality of transitions is a big advantage of the pulse clock. Proper encoding of 
the clock register is assured by operators |11) (11| which give an energy penalty to incorrect 
clock states with two active spots (two spin up particles) in the system. These are also 2- 
local. On the other hand, the big disadvantage of the pulse clock is that it is impossible to 
check whether the clock is not in the state \D) = |0 . . . 0)^, without a single spin up, using 
only local, projector terms. This 'dead' state is annihilated by the transition operators 
(j4.16p . so it does not 'move' anywhere. For this reason, the pulse clock is suited only for 
constructing Hamiltonian computers (see Section 12.2. ip , where one can pick the initial state 
of the clock and thus rule out the state \D). It is unusable for proving QMA-completeness 
results, unless one rules out the state \D) by adding terms that energetically favor a single 
spin up in the system over the state \D). 

I already presented Kitaev's unary domain wall clock in Section 11.3.21 

|t)^ = I11...1100...00). (4.17) 

The active spot in the clock register is now determined by the position of the domain wall 
10. The unique transition rules for this clock are now 3-local, |. . . 100 . . . ) |. . . 110 . . . ), 
but the advantage one gains are simple check operators for the clock. Illegal clock states 
are detected by the projector |01) (01 1 on consecutive clock qubits, and the dead state 
\D) = |0 . . . 0)^ can be easily ruled out using the projector |0) (0|^^ on the first clock qubit. 
Also, one can use a 1-local transition rule |0) <-> |1) between clock states as in |54j or [7], 
with the caveat that the illegal transitions such as 1 11100000) 1 11100100) can occur and 
need to be taken care of by using large penalty terms. I work around this problem while 
keeping the locality low in the following constructions. 

It is desirable to decrease the locality of transition rules as much as possible. In Bravyi 
[20] introduced a clock made out of particles with d = 4. The upside of his construction 
is that some of the transition rules involve only one clock particle (albeit 4-dimensional) . 
The interaction required to increment the clock and to apply a two-qubit gate at the same 
time thus involves particles with dimension 4x2x2, and can be thought of as a 4-qubit 
interaction. I now present a combined domain wall and pulse clock made out of particles 
with d = 2 and d = 3, building on the idea of [20]. Consider first a pulse clock made out of 
particles with dimension d = 3, with the progression of states 



|0), = !>□□□□), (4.18) 

|1), = !<□□□□), (4.19) 

|2)^ = !□>□□□), (4.20) 

|3)^ = !□<□□□), (4.21) 

|4)^ = !□□[>□□), (4.22) 

|5)^ = !□□<□□). (4.23) 



Some of the transitions are 1-local ( > < ) while others are 2-local ( < □ □ 1> ). There- 
fore, one can engineer a 3-local interaction for a qutrit and two qubits which increments 
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Figure 4-1: The 'triangle' clock transition rules allowing to apply a controlled gate to two 
work qubits using interactions that involve one work qubit at a time. 

the clock and concurrently applies a unitary gate to two work qubits. The initialization 
problem of this pulse clock is then solved by combining it with a domain wall clock 1 111000) 
as follows: 



The transition rules are 3-local for the transition of the active spot from the domain wall to 
the pulse clock 10 □ <-> 11 < , while they stay 1-local for the < <-> > transition. Also, note 
that 3-local operators suffice to detect whether the clock register state has a single active 
spot (either 10 in neighboring domain wall qubits or 1 < or 1 [> when there is a pulse 
clock particle at the domain wall) . I implement this clock using qubits rather than qutrits in 
Section 1121 which results in a new QMA-complete 3-local Hamiltonian construction without 
the requirement for large norm terms penalizing illegal clock states. 

Very recently, Eldar and Regev [28] had a revolutionary idea and abandoned the linear 
progression of states \^t) to create a 'triangle' clock as in Figure In this construction, 
the progression of states \^t) can take two different paths, depending on the state of one of 
the work qubits. Eldar and Regev use it to prove interesting results about Quantum 2-SAT 
with qudits, shattering the misconception that to apply a two-qubit gate, one has to interact 
with both qubits in the same step. I review their work in Section [4.4.11 and present my own 
results using this construction in Section r4.4.2[ I then take this idea further, creating a more 
complex, but more symmetric, train switch clock in Section [4.51 There I use it to prove the 
final result of this Chapter: one can simulate BQP in the Hamiltonian Computer model of 
Section 12.2.11 using a Quantum 3-SAT Hamiltonian. 

Geometric clocks 

Besides directly adding a separate clock register to the system and leaving the work qubits 
at a particular location in the system, one has the option of having the n work qubits move 
around by using higher dimensional particles. Consider a simple 2D construction using 
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qutrits, where the state space of each particle is 



TLn = 



O 



(4.31) 



with TCq a 2-dimensional subspace (can hold the state of a qubit) and TC^ a one-dimensional 
subspace. When a particle is in the state | Kl), it indicates that no qubit is present at the 
site. The location of the qubits along the strip then determines the states 



lO 

lO 



(4.32) 



Observe that the state \^t) is necessarily orthogonal to any state \^t') for which the lo- 
cation of the qubits is different. However, synchronization in this naive geometric clock 
construction is a problem, resulting in transition rules which are not unique and need to be 
at least 4-local, if one wants to apply two qubit gates Ut while moving the qubits forward. 
The challenge is to make a geometric clock encoding with unique transition rules which 
would involve only a few particles with low dimension. In other words, the goal is to get the 
work qubits to move along in a coordinated fashion and to have the corresponding gates Ut 
applied to them along the way. 

In the papers of Mizel [65], a geometric clock like this is planned to be implemented 
in a system of electrons hopping between quantum dots. The interactions in [65j involve 2 
particles hopping through 4 sites. 

Another geometric clock construction appeared in [7j , where the authors used 6-dimensional 
particles on a 2D grid. However, it doesn't have unique transition rules and thus requires 
large penalty terms for illegal clock states. Their progression of states |^'t) is 



]® □ □ 



]® □ □ 
]0Q □ 



® □ □ 
® □ □ 
(><)□□ 



]0 □ □ 
]® □ □ 



1® □ □ 
1 Q 
1 ^0 Q 



M □ 



where the state space of each 6-dimensional particle is a direct sum © Hq © W gi © W □ , 
with the 2-dimensional subspaces Wg) and Ti.Q able to hold the state of a qubit. The 
transitions in this type of clock can be made unique and checkable without large penalty 
terms if one makes the particles 12-dimensional. 

The recent construction of Aharonov, Gottesman, Irani and Kempe [5] shows that a 
geometric clock can be constructed even in ID, on a chain of 12-dimensional particles. 
What is even better, the legal clock states are locally checkable and the transition rules 
are unique and 2-local. The underlying idea in getting the qubits to move is a 'pass the 
hat forward and then jump over it' technique illustrated in Figure l4-2i I explain this 
construction in detail in Section 14.31 where I also decrease the required dimensionality of 
particles to d = 11. 

Janzing and Wocjan [48j also use a geometric clock in their complicated Hamiltonian 
Computer construction. Recently, myself and Wocjan [7l] simplified their model greatly, 
and use a novel geometric clock construction with transition rules reminiscent of diffusion 
instead of the conventional quantum walk on a line. I present this result in Chapter [5] on 
Hamiltonian Quantum Cellular Automata. 
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Figure 4-2: The 'pass the hat forward and then jump over it' technique for a geometric 
clock (qubit transport) on a line from [5j. The progression of 20 states is shown. 

4.2 A new 3-local QMA complete Local Hamiltonian 
(Quantum 5-SAT from 3-local terms) 

Quantum fc-SAT, the special case where the Hamiltonian is a sum of local projectors was 
defined and studied by Bravyi [20j as a natural analogue of classical /c-SAT. It is in P for 
k = 2, and it is QMAi complete for k > 4. However, the classification of quantum 3-SAT 
is still an open question. Its mother problem, 3-local Hamiltonian was shown to be QMA 
complete by Kempe and Regev in |54j. This result was further improved, showing that 2- 
local Hamiltonian is QMA complete in [53^ [74] . These constructions use the Hamiltonian of 
the form (j4.6p . and encode the clock register using the unary domain wall clock (j4.17p . Since 
the terms in their Hamiltonian are no longer 5-local as in Kitaev's proof in Section [L321 the 
corresponding terms in H^^^p do not only verify proper application of Ut, but also induce 
transitions from legal clock states into illegal ones. The subspace TLwork ® '^legai is thus no 
longer invariant under the action of H. To fix this, the penalty associated with illegal clock 
states is made high (scaling as a high polynomial in n), effectively forcing the ground state 
of the Hamiltonian to reside in the subspace of legal clock states. The Hamiltonian then 
contains 0{L) terms with weights that scale as 0(L^^) and has thus norm ||i?|| = 0(-L^^), 
where L is the number of gates in the computation [54] • This was later improved to weights 
of order O(L^), resulting in = 0{L'^) in [53]. 

Note that there are two energy scales in this problem. The norm of the Hamiltonian 
\\H\\, and the energy difference b — a > \/poly{L) in the definition of local-Hamiltonian. It 
is important to keep track of both of these. One can always rescale the Hamiltonian such 
that ||-ff|| is a constant, but this also shrinks the energy difference b — a. The quantity |^ 
is an indicator of the strength of the result. Physically speaking, it relates the precision 
required in obtaining the ground state energy (finding whether it is below a or above b) 
to the strength of the interactions (as given by Proving QMA completeness for the 

local-Hamiltonian where = const, would give a quantum analogue of the famous POP 
theorem [HI [11] ■ 

My goals are much simpler. I show in sect ion [4 . 2 . 5 1 that the 3-local Hamiltonian problem 
is QMA complete using a Hamiltonian with 0{L) terms with only constant operator norms. 
This makes the norm of my Hamiltonian scale as \\H\\ = 0{L), while I keep the same energy 
difference b — a = l/poly{L) as in the previous constructions. The O(L^) increase in the 
ratio of the two energy scales of the problem thus makes my construction physically much 
more interesting. 
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In this work, I show a new reduction from a verifier quantum circuit to a 3-local Hamil- 
tonian. The novelty of my construction is that it leaves the space of legal clock-register 
states invariant. Therefore, the weights of the terms in my Hamiltonian do not scale with 
the size of of the input problem. Such terms do appear in the constructions of [54] and 
[53j . As an intermediate step in the construction, I prove that quantum 3-SAT for qutrits 
is QMAi-complete. 

The rest of Section 14.21 is organized as follows. After reviewing Bravyi's proof that 
Quantum 4-SAT is QMAi-complete in Section 14. 2. 1^ I present a qutrit-clock construction 
in Section 14.2.21 and show that quantum 3-SAT for particles with dimensions 3x2x2 
(the interaction terms in the Hamiltonian couple one qutrit and two qubits) is QMAi- 
complete. The existence of such a construction was previously mentioned but not specified 
by Bravyi and DiVincenzo in [20] as [27]. In Section [4.2.51 1 show how to encode the qutrit 
clock particles from Section 14.2.21 into a pair of qubits in such a way that the Hamiltonian 
remains 3-local, obtaining a new construction of a QMA complete 3-local Hamiltonian. This 
Hamiltonian is composed of 4-local positive semidefinite operator terms. However, each of 
these 4-local operators is composed of only 3-local interaction terms. It is not a Quantum 
3-SAT Hamiltonian, since the 3-local terms by themselves are not positive semidefinite 
operators. I discuss the complexity of Quantum 3-SAT and further directions in Section 

MM 



4.2.1 Bravyi's Quantum 4-SAT 



In [20j, Bravyi proved that quantum fc-SAT belongs to QMAi for any constant k. Further- 
more, he showed that quantum 4-SAT is QMAi complete using a new realization of the 
clock. In his construction Bravyi uses L + 1 clock particles with 4 states: unborn, active 1 
(ai, input for a gate), active 2 (02, output of a gate), and dead. These 4 states of a clock 
particle are easily realized by two qubits per clock particle. There are 2L legal clock states: 

\C2k-1) = \d_^_^^^aiu^^^^^ji), and \C2k) = \ d_^_^^^a2 , (4.33) 

k-l L-k k-l L-k 

for 1 < < L. A clock Hamiltonian Hdock = Hdockinit + Y^k=i ^dLk i"equired to check 
whether the states of the clock are legal. 

H^Lk = \d){d\,^j\u){u\,^, (4.34) 



+ 



+ 



fc+1 



(""Ifc ^ ( \d) {d\ + |ai) (ai| -|- \a2) (02! )^ 
\ai) (ail + \a2) («2| )^ i?) \ai) (aij + \a2) (02! + \d) 



fc+1 

Hdockinit = \u) (till + \d) {d\^. (4.35) 
The Hamiltonian checking the correct application of gates is Hp^op = X]fc=i ^prop^ with 

Hjyrlp = ^ ® ("ilfc + ^ ^ l«2) {a2\k -Uk® \a2) {ai\t^ -Ul® \ai) (02!^ ) . (4.36) 

Each such term verifies the correct application of the gate Uk between the states \ai) and 
[02) of the A;-th clock particle. This only requires interactions of the A:-th clock particle 
(qubit pair) and the two work qubits the gate C/fc is applied to. Each of the terms is thus a 
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4-local projector. 

I need another Hamiltonian term to propagate the clock state \C2k) into |C2fc4-i) while 
leaving the work qubits untouched (that is, for the ground state \il^2k)work = \'4^2k+i) work) ■ 
This is done by the 4-local clock-propagation Hamiltonian Hdockprop = I]fc=i ^dlckprop^ 



with 



Hfockprop = ^(l«2)(a2U^|n)Hfc_,, + |d)(d|,0|ai)(ai|,+i) (4.37) 



2 

The final ingredients in this construction are 

^nit = 5];|1)(1U® |ai)(ai|i, (4.38) 

n=l 

Hout = |0)(0|„„,®|a2)(a2L. (4.39) 

Applying Kitaev's methods [SS] to this construction, Bravyi shows that the quantum 4-SAT 
Hamiltonian (a sum of 4-local projectors) 

H — HqIqi-i^ -\- H (.[Qf^if^p^Qp -\- Hinit ~\~ Hq^i -\- Hprop (4.40) 

is QMAi complete. 

Bravyi's original definition required all of the terms in the Hamiltonian to be projectors. 
However, as I have shown at the end of Section [4.1.H using positive semidefinite operator 
terms in H instead of just restricting ourselves to projectors is an equivalent problem. 



4.2.2 A qutrit clock implementation 

In this section I present a new realization of the clock which builds on Bravyi's quantum 
4-SAT realization described above. Using this clock construction, I prove that quantum 
3-SAT for qutrits is QMAi-complete. 

First, I need to show that quantum 3-SAT with qutrits is in QMA. I can use Bravyi's 
proof that quantum /c-SAT for qubits is in QMAi for any constant k. Given an instance of 
quantum 3-SAT for qutrits, I convert it into an instance of quantum 6-SAT for qubits by 
encoding each qutrit in two qubits and projecting out one of the four states. According to 
Bravyi, this problem is in QMAi and therefore so is the original quantum 3-SAT problem 
with qutrits. 

For the other direction in the proof, I need to construct a quantum 3-SAT Hamiltonian 
for qutrits, corresponding to a given quantum verifier circuit U for a problem in QMA. The 
terms in the Hamiltonian I will construct act on the space of one qutrit and two qubits 
(particles with dimensions 3x2x2). 



4.2.3 Clock register construction 

The clock-register construction in the previous section required 4 states for each clock 
particle: \u) , |ai) , jaa) and \d). Let me first explain why Bravyi's construction requires two 
"inactive" states: \d) and \u). If I only use \d) (i.e., have legal clock states of the form 
\d. . . daid . . .d) and \d . . . da2d . . . d)), I immediately get a 3-local Hamiltonian for qutrits. 
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Figure 4-3: Clock register consisting of 2L qubits and L qutrits. 



However, in Bravyi's construction, the first clock particle is never in the state \u), and the 
last one is never in the state \d) (see (14.35P ). This ensures that at least one clock particle 
is in an active state. When not using the state I can no longer exclude the state with 
no active particles \dd . . . d) in a simple local fashion. 

I fix this by modifying the clock register as shown in in Fig |4-3i The clock register now 
consists of 2L qubits and L qutrits. The 2L qubits ci,di, . . . ,Ck,dk play the role of the 
usual unary |1 . . . 1100 ... 0) clock representation, while the L qutrits ti,t2, ■ ■ ■ play the 
role of Bravyi's clock with just three states {d, ai, 02). 

I define the legal clock space TLiegai as the space spanned by the 2>L states \Cm)- These 
states are defined for \ < k < L as follows: 

|C73fc_2) = \{lld){lld)...{lld){md){md){md)...{md)), (4.41) 

■V ^ v 

k—l times L—k times 

|C3fc-i) = |(lld)(lld)...(ll(i)(llai)(00(i)(00d)...(00d)), 

k—l times L—k times 

ICgfc) = |(lld)(ll(i)...(ll(i)(lla2)(00d)(00(i)...(00d)). 

V V 

k—\ times L—k times 

The first state (|C3fc_2)) corresponds to the time when the qubits are "in transport" from 
the previous gate to the current (/cth) gate. The second one corresponds to the time right 
before application of gate C/^ and the third corresponds to the time right after the gate 
was applied. The structure of such clock register can be understood as two coupled "unary" 
clocks, the qubit one (c^, d^) of the 11 . . . 11100 ... 00 type and the qutrit one (the t^'s) of 
the 00 . . . 00100 ... 00 type. Formally, the legal clock states satisfy the following constraints: 

1. if dfc is 1, then is 1. 

2. if Cfc+i is 1, then is 1. 

3. if tfc is active (01/02), then dfc is 1. 

4. if ifc is active (01/02), then c^+i is 0. 

5. if dfc is 1 and Cfc_(_i is 0, then is not dead (d). 

6. ci is 1. 

7. if is 1, then is not dead. 

The last two conditions are required to exclude the clock states |(00(i)(00d) . . . (OOd)) and 
|(ll(i)(lld) . . . (lid)) that have no active clock terms. The clock Hamiltonian Unlock = 
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Figure 4-4: Illustration of the two-step clock pointer propagation. 



Hdockinit + Y.k=i ^dLki + I]fc=i ^dLk2 verifies the above constraints. 



rrik) 
^clockl 



rrik) 



|01) (01|,^^,^ + |0) (0|,^ (|ai) (ail + la^) (asj) , (4-42) 
|01) (01|,„e,^., + ( (ail + |a2> (aal |1> 

Hclock^nit = |0)(0|,^+|l)(l|,^^|<i)((i|,^. 
(k) 

Only the last term in H^i^^^2 ^ 3-local projector, acting on the space of two qubits and 
one qutrit. The rest of the terms are 2-local projectors on two qubits, or a qubit and a 
qutrit. The space of legal clock states Tliegai is the kernel of the clock Hamiltonian Hdock- 



4.2.4 Checking correct application of gates and clock propagation 

The gate-checking Hamiltonian Hprop = J2k=i ^prop is an analogue of (|4.36|) . with 
H^rlp = ^(Wfc® {ai\t^ + \a2) (a2lt, ) - Uk ^ jaa) (oil^^ - ul |ai) (aal^^ ) . (4.43) 

The clock propagation proceeds in two steps. First, the "active" spot in the clock 
register moves from the state |a2) of the qutrit tk to the |10) state of the next two qubits 
Cfe+i, dk+i- After this, it moves into the state \ai) of the next qutrit t^+i, as in Fig l4-4[ The 
Hamiltonian checking whether this happened, while the work qubits were left untouched, is 

Ak) . rrik) 



^clockprop — l^k=l ^dockpropl Z^fc=l ^dockprop2^ 

^( (lOlc.A ^ l'^) (d\t, + 111) (ll|c,A ® ) (4-44) 

^( 111) (lOlc.A ^ l«i) (d\t, + |10) (ll|e,A ® 1^) ("il^O' 



rrik) 

clockpropl 



rrik) 

clockprop2 



i( |a2) (a^l,^ 100) (001,^^^,,^^^ + \d) {d\,^ |10) (10|,^^^,,^^^ 
l[\d) (a^l,^ |10) (001,^^^,,^^^ + |a,) {d\,^ 100) (10|,^^^,,^^^ 



The input Hamiltonian checks whether the computation has properly initialized ancilla 
qubits. 

Na 

Hinit = 5^|1)(1|„® |ai)(ai|,^. (4.45) 

n=l 
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Finally, the output Hamiltonian checks whether the result of the computation was 1. 



out 



|0) {^\out ® l«2) (a2 



(4.46) 



All of the terms coming from (j4.42p - (j4.46p in the Hamiltonian 




prop- 



(4.47) 



are projectors. Therefore, the ground state has energy zero if and only if there exists a zero 
energy eigenstate of all of the terms. If there exists a witness \ip) on which the computation 
U gives the result 1 with probability 1, I can construct a computational history state (14. 7p 
for a modified circuit U = [7^ • I • I • C/l-i • I • I • • • C/i • I, where the "identity" gates correspond 
to the clock propagation in my construction, with nothing happening to the work qubits. 
This state is a zero eigenvector of all of the terms in the Hamiltonian (j4.47p . 

I now need to prove that if no witness exists (the answer to the problem is "no" ) , then 
the ground state energy of (j4.47p is lower bounded by l/poly{L). Let me decompose the 
Hilbert space into 



where TLiegai is the space of legal clock states (on which Hdock |a) = 0). The Hamiltonian 
(|4.47p leaves this decomposition invariant, because it does not induce transitions between 
legal and illegal clock states. Since any state in Ti^g^i violates at least one term in Hdock, 
the lowest eigenvalue of the restriction of ()4.47p to Tiwork ^ '^'tgai least 1. On the 

other hand, the restriction of H to the legal clock space is identical to the legal clock space 
restriction of Bravyi's Hamiltonian ()4.40p from the previous section. Therefore, his proof 
using the methods of Kitaev [56j applies to this case as well. He shows that if a no witness 
state for the quantum circuit U exists, then the ground state energy of the restriction of 
(j4.40p to 'Hwork^'^iegai lower bounded by l/poly{L). This means that if there is no witness 
state for the verifier circuit U, the ground state of (j4.47p is lower bounded by l/poly{L). 
This concludes the proof that quantum 3-SAT with qutrits (in fact, quantum 3-SAT on 
particles with dimensions 3 x 2 x 2, a qutrit and two qubits) is QMAi complete. 

The existence of another 3x2x2 construction for quantum 3-SAT (i.e., a Hamiltonian 
with terms acting on one qutrit and two qubits) was already mentioned in [20] as [27j . 
though that construction was not specified. I choose to write out my result explicitly as 
it serves as a natural intermediate step towards the new 3-local Hamiltonian construction 
described in the following section. 



4.2.5 The new 3-local QMA complete construction (for qubits) 

In Bravyi's Quantum 4-SAT construction [20], the clock particles (qubit pairs) can be in 
4 states. In the previous section, I required only 3 states of the clock particles and used 
qutrits as particles with these three states. I start with the clock-register construction (see 
Fig l4-3p from the previous section, with legal states as in (j4.4ip . However, I now encode 
the three states of every clock qutrit tk using a pair of qubits rk,Sk- The new clock register 



W — iJ~(-work ® Ti-legal) © {'Hwork ® '^legal)- 



(4.48) 
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Figure 4-5: Clock register construction with 2L + 2L qubits. 



is depicted in Fig |4-5[ 

^ ;7f l*^)*. ^100).,,.,, (4.49) 

This encoding allows me to obtain a new 3-local Hamiltonian construction (for the QMA- 
complete 3-local Hamiltonian problem) . This Hamiltonian is a quantum 4-S AT Hamiltonian 
whose 4-local positive semidefinite operator terms consist of just 3-local interactions. 

I am looking for 3-local terms that flip between the clock states |ai) |a2), while 
simultaneously (un)applying a 2-qubit gate Uk on two work qubits. I encode the active 
states of a clock particle into entangled states, and thus I am able to flip between these 
clock states with a term like Zi involving only one of the clock particles. Thus my 2-qubit 
gate checking Hamiltonian involves only 3-local terms (acting on one clock qubit or Sk 
and the two work qubits on which the gate acts). 

First, I define the legal clock space TLiegai as the space spanned by the 3L states \Cm)- 
These states are defined for 1 < A; < L as follows (compare to (|4.4ip ): 

\C^k-2) = I (11)(00) . . . (ll)(00)j ® |10)e,A ® I00)r.,s. ^ I (00)(00 ) • • • (00)> (4.50) 

fc— 1 times L—k times 

\Czk-i) = I (11)(00) . . . (11)(00)) ® |ll),„rf, ® (|01) - |io)),,,,, ® I (oo)(oo) . . . (00) ), 

k—l times L—k times 

\C^k) = I (11)(00) . . . (11)(00)) ® |ll),,,rf, ® -)= (|01) + 110)),^^,^ ® I (00)(00) . . . (00)). 

^ V ' v2 ^ V ' 

fc— 1 times L—k times 

Similarly to the construction of the previous section, the first state {\C2k-2)) corresponds 
to the time when the qubits are "in transport" from the previous gate to the current (/cth) 
gate. The second one corresponds to the time right before application of gate and the 
third corresponds to the time right after the gate Uk was applied. 

Formally, the legal clock states for this construction satisfy the following constraints: 

1. if dk is 1, then is 1. 

2. if Ck+i is 1, then dj. is 1. 

3. the pair rfc,Sfc is not in the state |11). 

4. if the pair r^, Sk is active (in the state (|01) ± |10))/\/2), then dk is 1. 

5. if the pair r^, Sk is active (in the state (|01) ± |10))/-v/2), then Ck+i is 0. 

6. if dk is 1 and Ck+i is 0, then the pair rfc, Sk is not dead (in the state |00)). 
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7. ci is 1. 



8. if is 1, then the pair vl, sl is not dead (in the state |00)). 

The last two conditions are required to make the clock states |(00)(00)... (00) (00)) and 
I (11) (00) . . . (11) (00)) with no active spots illegal. The clock Hamiltonian Hdock = Hciockinit+ 
ELi ^dicfei + Efe=i^ ^docfc2 verifies the above constraints. 

HiiLki = |01)(01L„,, + |0)(0|,^® (|1)(1|,^ + |1)(1|,J + |11)(11|,^_,^, (4.51) 



hSL, = |oi)(oi|,„,,^^ + (|i)(i|,, + |i)(iLj®|i)(i|,,^^ + /.f , 

hf^ = |l)(lk®^(^r. + ^^J®|0)(OU+, + |ll)(llL..> 
Hclockinit = |0)(0|,^ + |1)(1|^^®|00)(00|,^_,^. 



(k) (fc) 

All of the terms involve only 3-local interactions. All terms in -ffc/ocfci' -^docfc2 ^^'^ Hdockinit, 

(fc) 

are are projectors. The term h\ corresponds to the sixth legal state condition. It is a 4- 
local projector onto the space spanned by (illegal clock) states \ ldi^{00)r^,,sk^ck+i) ^ 
|0(11)1) and 1 1(1 1)1). Note that even though ^4^^^ is a 4-local projector, it is only constructed 
of 3-local terms. 

4.2.6 Checking gate application with 3-local terms 

Let me start by writing out a Hamiltonian that checks the correct application of a single- 
qubit gate Uk- 



rr(k).one-,um ^ 1 / I ® |01 " 10) (01 - 10|,^,,^ - ® |01 + 10) (01 - 10|,^,,^ 

2 i I ® |01 + 10) (01 + 10|,^ - ul ® |01 - 10) (01 + 10|,^ 



,(4.52) 



where |01 it 10) is a shortcut notation for the normalized entangled states (|01) it |10))/-\/2. 
This Hamiltonian is a 3-local projector. Note that in the case Uk = I, this Hamiltonian 
becomes the projector (I — X)/2 on the space of active clock states {|01 — 10) , |01 + 10)}. 

For a two-qubit gate C/fc, the above construction would be 4-local. However, I am be able 
to construct this 4-local projector using only 3-local terms. To do this, I require the 2-qubit 
gate to be symmetric"' Uk = u\. This is a universal construction, since the symmetric gate 
CNOT (or C<^) is universal. Now I can write 

HW,tn.o-,ulnt ^ 1 ^ 1 (I _ z,^z,,) -Uk^l {Zr, - Z.j) . (4.53) 

The first term in this Hamiltonian, (I — Z^-^Zs,^) /2, is a projector onto the space of ac- 
tive clock states, |01 it 10)^^ as I needed. The second term contains (Z,.^ — Zg^)/2^ 
which has zero eigenvalues for the states lOO)^ ^ and 1 11)^ „ , and fiips between the states 
|01 — 10) |01 -I- 10) . Altogether, this is a 4-local projector made out of only 

3-local terms. 

4.2.7 Clock propagation 

After a gate Uk is applied, I need to "propagate" the pointer (the active state of the qubit 
pair rk,Sk) to the next pair of qubits rk+i,Sk+i- This is done in two steps, as shown in 
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|oi+iq| (o) (o) < — ► 00 (i) (o) 



00 < ► (1) (1) 01-10 



Figure 4-6: Illustration of the two-step clock pointer propagation. 

FiglMl 

For each step, I want to write a 3- local positive semidefinite Hamiltonian with terms 
acting on 4 consecutive qubits rk,Sk,Ck+i,dk+i (for the second step of the clock pointer 
propagation, the four qubits in play are Ck, dk,rk, Sk), with zero eigenvalue for the legal 
clock-propagation states, and perhaps also some illegal clock states, which will be disallowed 
by other terms in the Hamiltonian {Hdock)- For the first step, these desired eigenvectors 
with zero eigenvalues are 

l«i).„.„c,+,A+i = |00)r„sJOO),,^,,,,^, , (4.54) 
l«3>.„.„,,^„,,^, = ^(101) + 110)).,,., 100),^^^,,^^, + i= 100),^,,^^ |io),,^,,,,^, , 

The state that I want to exclude (make it a nonzero eigenvector) is the legal clock state 
with incorrect pointer propagation: 

\»^)r.,-4.,. = ^(101) + |10)k,«J00>.,,„.,+i-^|00),^,,J10),^^^,,^^^. (4.55) 

Let me present the Hamiltonian penalizing this state. 

H^Lpropi = Il0>(10|,^^^,,^^^ (4.56) 
+ i(|01) + |10))((01| + (10|),^^,^ 

1 



V2 
1 

+ 2|11)(11 



|o)(iU + |o)(iLj0|io) (001,^^^,,^^^ 



^(|i)(0U + |i) (olJ^Ioo) (10|,^^^,,^^^ 



It is positive semidefinite, with eigenvalues (x7), 1 (x4), 2 (x3) and 3 (x2). Its zero 
energy eigenvectors are |ai), |a2), jas), |a4) expressed above, and three illegal clock states, 
|00)^^^^^ |01>c,+i,d, 1, (|01)-|10)) |01) and (|01)-|10)) |11). The state \a-^) is an eigenvector 

(k) (k) 

of H^iockpropi ^'^^^ eigenvalue 2. This means that Hl^ockpropi fulfills its job in punishing the 
legal states of the clock register ()4.55p . which do not correctly propagate the clock. This 
Hamiltonian term is a positive semidefinite operator, while Bravyi's original definition of 
quantum fc-SAT requires the terms in the Hamiltonian to be projectors. However, as I have 
shown at the end of Section [4.1.1l quantum fe-SAT with positive semidefinite operator terms 
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is equivalent to quantum A;-SAT with only projector terms. 

For the second step, the desired zero energy eigenvectors are 



1 



and the state I want to exclude is 



l/3^)c,A,r..,.. = ^|10),^^,J00),^,,^,-|11),^,,^-(|01)-|10)),„,,. (4.58) 



The Hamiltonian with these properties is a simple analogue of (j4.56p : 

HSLprop2 = Il0)(10|,^,,^ (4.59) 

+ ^(ioi)-|io))((oi|-(io|),^,,^ 



^2^ 

+ 2 111) (111.,,.,. 



This is again a positive semidefinite operator with eigenvalues (x7), 1 (x4), 2 (x3) and 3 
(x2). Its zero energy eigenvectors are |/3i), \(32), IPs), lA) expressed above, and three illegal 

clock states, lOl)^^^^^ \^^)r^„Sk' l^^) d^l) + l^^)) 1^^) d^l) + l^^))' ^^^^^ penalized 

I I \ (ft) 
by H clock- The state |/5 ) is an eigenvector of Hd^ckpropi '^ith eigenvalue 2, which is what 

(k) 

I intended H^i^^^^^^ to do. 

Just as in the previous section, the total Hamiltonian leaves the decomposition into 
T~iiegai ® 'H'tgai invariant while all illegal clock states violate at least one term in Hdock- 
Again, up to a constant prefactor, the restriction of H to the legal clock space is the same 
as that of the Hamiltonian in ()4.40p . The proof of the necessary separation between positive 
and negative instances then follows the proof in the previous section. This concludes the 
proof that quantum 4-SAT with positive semidefinite operators made out of 3-local terms 
is QMAi complete. 



4.2.8 Discussion and further directions 

In the above, I proved that quantum 3-SAT for particles with dimensions 3 x 2 x 2 is QMAi 
complete. I have shown in Section 14.1.11 that quantum A;-SAT with positive semidefinite 
operator terms (not just projectors) is equivalent to quantum A;-SAT with projector terms. 
I presented a new 3-local construction of a quantum 4-SAT Hamiltonian with positive 
semidefinite operator terms, proving that quantum 4-SAT with 3-local interactions is QMAi 
complete. 

The currently known complexities of classical and quantum satisfiability problems are 
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shown in Table 14.11 Quantum 3-SAT contains classical 3-SAT and therefore is NP-hard. 
Unlike classical /c-SAT, which is known to be NP-complete for k > 3, quantum /c-SAT is 
only known to be QMAi complete for k > 4 [20]. In my opinion, it is unlikely that one can 
show that quantum 3-SAT (k = 3) is also complete for QMAi. One indication for this arises 
in my numerical explorations, where random instances of quantum 3-SAT for a reasonable 
number of clauses generally have no solutions, unless the clauses exclude non-entangled 
states. This may suggest that the hardness of quantum 3-SAT actually lies only in the 
classical instances (3-SAT) and a classical verifier circuit for quantum 3-SAT might exist. 

Another reason comes from dimension counting. This argument, however, is only valid 
for the specific encoding of a circuit into the Hamiltonian I used. I worked with a tensor 
product space 'Hwork^'^dock-, encoding the computation in the history state (j4.7p . Encoding 
an interaction of two qubits requires at least an 4 -|- 4 = 8 dimensional space (4 for the two 
qubits before the interaction and 4 for the qubits after the interaction). On a first glance, 
a three-local projector on the space of two work qubits and one clock qubit (2x2x2 = 8) 
seems to suffice. However, I must ensure that this interaction only occurs at a specific 
clock time. When the two work qubits interact with just a single clock qubit (ffipping it 
between states before/after interaction) ambiguities and legal-illegal clock state transitions 
are unavoidable. This transforms the problem from the SAT-type (determining whether a 
simultaneous ground state of all terms in the Hamiltonian exists) to the MAX-SAT type 
problem (determining the properties of the ground state energy of the sum of terms in the 
Hamiltonian). A single clock qubit cannot both determine the exact time of an interaction 
and distinguish between the states before and after the interaction. I managed to overcome 
this obstacle by using three-dimensional clock particles with states d, oi and 02 and a 
2 X 2 X 3 = 12 dimensional space for encoding the interactions. I believe that this can 
not be further improved with more clever clock-register realizations within the usual 7i = 
T~(-work ® T~(-ciock framework. However, a recent novel idea by Eldar and Regev [28] presented 
in Section 14.4.11 could be a step in this direction. They introduce a novel 'triangle' clock 
construction, and show that Quantum 2-SAT for particles with dimension 5 x 3 is QMAi- 
complete. 

A different approach to encoding a quantum computation into the ground state of a 
Hamiltonian is the geometric clock (see Section 14.1.31 found in the work of Aharonov et.al 
[7]. Their idea is to lay out the qubits in space in such a way, that the shape of the state (the 
locations of the work qubits moving around in the system) uniquely corresponds to a clock 
time. Their motivation was to show that adiabatic quantum computation [7] is polynomially 
equivalent to the circuit model. As a side result, they showed that the nearest-neighbor 
2-local Hamiltonian problem with 6-dimensional particles is QMA complete. Actually, 
as Kempe et al. proved, even 2-local Hamiltonian with 2-dimensional particles (qubits) 
is QMA complete ^3]. I did not succeed to improve or reproduce the construction of 
quantum 3-SAT for particles with dimensions 3x2x2 using the geometric clock framework. 
However, one can use the idea of a geometric clock to construct quantum 2-SAT for higher 
dimensional particles (qudits). I know that classical 2-SAT for particles with dimensions 
3x3 contains graph coloring, and is thus NP-complete. Using a rather straightforward 
modification of the Aharonov et. al. construction, one can prove that quantum 2-SAT for 
1 2-dimensional particles is QMA complete. Combining the work/clock and the geometric 
clock constructions, a much tighter result can be shown. Specifically, one can construct 
quantum 2-SAT for particles with dimensions 9x4 and prove that it is QMAi complete. 
However, this is now made obsolete by the result of Eldar and Regev presented in Section 
14.4.11 which does not use a geometric clock. It remains to be seen what are the minimal 
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dimensions of particles for which quantum 2-SAT is QMAi complete. 



4.3 Quantum 2-SAT on a line 

In 2007, Aharonov, Gottesman, Irani and Kempe |5j proved a surprising result about the 
power of quantum systems on a line. One might have thought before that finding ground 
states or simulating time evolution of spin chains is not too hard a task with today's modern 
DMRG and MPS methods (see Chapter[3]). However, [5j shows that finding ground states of 
2-local Hamiltonians on a line of d = 12 dimensional particles is QMA complete. Utilizing 
this construction, recently Schuch, Cirac and Verstraete [84] proved that even finding ground 
states of 2-local Hamiltonians on a line whose ground states are known to be Matrix Product 
States with constant dimension is NP-hard. 

In this section I revisit [5] and show that the problem remains QMA complete even 
for particles with dimension d = 11. Moreover, [5] shows that one can use a simplified 
nearest-neighbor Hamiltonian on a line of d = 9-dimensional particles to perform universal 
Quantum Adiabatic Computation. This result was recently improved to a system of 8- 
dimensional particles by Chase and Landahl [22]. Concurrently with their work, myself 
with P. Wocjan in [71] have proved that a translationally invariant Hamiltonian on a line 
of 10-dimensional particles can be used for this purpose. When I release the translational 
invariance requirement, I get a Hamiltonian for a chain of 8-dimensional particles as well. 
I present this result later in Section [5. 2[ 



4.3.1 A Line of (i = 12 Dimensional Particles (Quantum (12,12)-SAT) 

As discussed in Section [4. 1.31 to show that a Local Hamiltonian (or Quantum-SAT) problem 
is QMA complete, one usually proceeds as follows. First, encode a quantum circuit into a 
progression of orthogonal states \^t) of a larger system, either by adding a clock register 
to the register holding the work qubits, or by using a geometric clock construction. The 
corresponding Hamiltonian is then constructed to check the proper transition rules for the 
progression of states \^t)- What remains to be shown are the properties of the low- lying 
spectrum of H. 



The d = 12 Geometric clock on a line 

On a line, only nearest neighbor interactions are allowed. This rules out the use of a separate 
clock register, which means a geometric clock (see Section [4. 1.3p is a necessity. This is also 
a reason for the dimensionality {d = 12) of the particles. The state space of each particle 
consists of four two-dimensional subspaces 



o 



a qubit marked as 'active', 

a qubit labeled 'messenger', 

a qubit to the left of the active site, 

a qubit to the right of the active site 



(4.60) 
(4.61) 
(4.62) 
(4.63) 
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Figure 4-7: A quantum circuit for n = 3 qubits with K = 4 rounds of nearest neighbor 
gates Uk,s- 



and four more states, each of them 1-dimensional 

M : to the left of the qubits ('done/dead'), (4.64) 

□ : to the right of the qubits ('ready/unborn'), (4.65) 

O : a 'turn' state at the qubit sequence boundary, (4.66) 

<l : a 'push' state used to move the qubits to the right. (4.67) 

Altogether, the state space of each particle is 12-dimensional, with the structure 

ni2=Qeh = {q(S)l)®h = ( O X,-}) ©{K, 0,0, <}, (4.68) 



where Q is the subspace for the 4 types of qubits with q the internal state space of a qubit 
and / its four possible labels. The state space for the four other states is labeled h. 

Without loss of generality, take a quantum circuit on n qubits which consists of K 
rounds of nearest neighbor gates as in Figure 14-71 

U = {UK,n-lUK,n-2 ■ ■ ■ Uk,i) • • • {U2,n^lU2,n~2 • • • f/2,1) (C/l,n-lC/l,n-2 • • • C/l,l) , (4.69) 

where the gate Uk^s acts on the pair of qubits s, s + 1. I now encode the progression of U 
into a set of states of a line of qudits with length nK. The initial state IV'i) has a sequence 
of n qubits on the left, and the rest of the chain is in the 'ready' state □. The leftmost 
qubit is active . Here I write out the initial state corresponding to a quantum circuit on 
n = 3 qubits with (7 = 4x2 nearest neighbor gates (here K = 4, as in Figure . 

= (gOOH H H H H H H H H (4-70) 

Using a few rules, a progression of states \xpt) is now constructed. The rules are prepared 
in such a way that there is always only one possible state \ipt+i) for a given state {tpt)- 

The first rule says that when the active spot is at position kn + m where 1 < k < K 
and 1 < m < n, it can pass to the right. When this happens, the gate Uk,m is first applied 
to the state of the corresponding qubits {kn + m,kn + m + 1), while their labels change 
from ® □ to 0®. 

la : ®0 ^ Uk,m{®^) (4.71) 
on particles {kn + m,kn + m + 1). 

On the other hand, the active spot @ just moves to the right, changing the qubit labels 
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without modifying their internal states. 

16 : ©O ^ O© (4.72) 

Here is the corresponding progression of states (only la applies here): 

(g)OOQ H H H H H H H H (4.73) 
®(g)OQ H H H H H H H H (4.74) 
©©(§)□□□□□□□□□ (4.75) 

After these n — 1 applications of rule la, the internal state of the n qubits holds the state 
of the quantum circuit after the first round of gates. 

The second set of rules involves the active spot reaching the front of the chain of qubits. 
Rule 2a applies to particle pairs (kn — 1, kn) for integer k, while rule 26 applies everywhere 
else. The reason for this is to ensure that when working out the progression of states 
backwards, there is always only one of the rules that applies in a given situation. 

2a : (g) □ — > ®0 on particles (/cn - 1, fcn), (4.76) 
26 : ©□ — > ©0 everywhere else, (4.77) 

Thus, after applying rule 2a on particles (1,2), I get 

©©© □□□□□□□□□ (4.78) 
©©©O □□□□□□□□ (4.79) 

The next three rules facilitate the sending of a message back to the left end of the chain. 

(4.80) 
(4.81) 
(4.82) 

Also note that rule 3 is simply O 
is <\ — > (3 for the leftmost particle of the line. This produces the state progression 

©©©O □□□□□□□□ (4.83) 

©©© (4.84) 

©©<©□□□□□□□□ (4.85) 

©<©©□□□□□□□□ (4.86) 

<©©©□□□□□□□□ (4.87) 

(^©©©□□□□□□□a (4.88) 

Observe now that the qubits have moved one step to the right from where they were in 
(j4.73p . Finally, the last set of rules activates the leftmost of the qubits. Depending on the 
position this qubit is at, it will become either (g or @. 

6a : O© — > for particles (A;n,A;n + 1), (4.89) 

66 : O© — > (4.90) 



3 


□ 


— . < □ 


4 


©< 


<© 


5 


M < 


— > MO 


- 




for the ris 
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The progression of states thus continues as 



and 



OOOOQ H □□□□□□ 


(4.91) 


^©OOQ H □□□□□□ 


(4.92) 


K®@OQ □ □□□□□□ 


(4.93) 


K®®(g) □□□□□□□□ 


(4.94) 


□□□□□□□ 


(4.95) 


^00® 


(4.96) 


^0® <OQ H □ □ □ □ □ 


(4.97) 


^® <OOQ Q □ □ □ □ □ 


(4.98) 


^ <OOOQ Q □ □ □ □ □ 


(4.99) 


^OOOOQ Q □ □ □ □ □ 


(4.100) 


K K@OOQ Q □ □ □ □ □ 


(4.101) 


K K(@OOQ □ □ □ □ □ □ 


(4.102) 


K K®(g)OQ H □ □ □ □ □ 


(4.103) 


K K0®(gQ □ □ □ □ □ □ 


(4.104) 


m □□□□□□ 


(4.105) 


m ^00® 


(4.106) 


K K®® <OQ H □ □ □ □ 


(4.107) 


M <OOBB □ □ □ □ 


(4.108) 


M M <000 □□□□□□ 


(4.109) 


K KOOOQH H □ □ □ □ 


(4.110) 


M M K(g)00 □□□□□□ 


(4.111) 



Note that in steps (j4.92p and (j4.10ip the last qubit is marked as a 'messenger' @ by rule 
6b. This rule is used at the left end of the qubit sequence until the qubits have moved n 
steps to the right. Then, in (j4.11ip . the leftmost qubit is at position kn + 1, which makes 
it become ® by rule 6a. The second round (A; = 2) of gate applications now begins. 



The Hamiltonian 

As in Section [4.21 the Hamiltonian is an implementation of (j4.6p . 

The part checking the propagation of the computation is constructed from the above 
rules as 

riK-l 6b 
s=l r=la 

where each term R^^^ is a projector acting on one or two neighboring qubits corresponding 
to the rule r. 

As an example, I write out the term R^^""^ which checks the proper application of rule la. 
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This projector acts nontrivially on an 8-dimensional subspace of two neighboring particles 
spanned by 0@- Using the notation of (I4.68p . this subspace can be thought of 

as 



(gi fe® ^2), (4.113) 

where qi and q2 are states of the two qubits, while their labels (Zi, I2) are in the state (►, •) or 
(x , ►). The corresponding projector (only for clock particles in position kn + m, kn + m+1) 
is 

= ® (!►•)(► -I + |x^)(x^|U (4-114) 

- liUk,m)^^^^^0\x>){>-\i^i^ (4.115) 



91,92 

where the gate Uk,m is acting on the internal state of the two qubits. 

The rest of the projectors are much simpler, therefore I write out only one of them. 
Rule 2a involves the 4-dimensional subspace of two particles spanned by □ and M O ■ 

Ptn-iM) = (4-117) 



The low-lying spectrum of H 

What remains in the Hamiltonian construction are the legal clock state checking operators. 
First, analogously to the |01) (01 1 check operators of the domain wall clock, 1 need 

Hcheckoi = 1^) (^1 ® 1^) (^1 (4-118) 

A&Sa B&Sd 

+ E E 1^) ® ' (4-119) 



where 



Sa = {©,®,0,<}, (4.120) 

Sd = {K,®}, (4.121) 

Sr = {B,Q}. (4.122) 

Next, I need to check that there is an active site in the system at all by 

HchecklaO = E E ^ - (4-123) 

Next, 1 need to check that there is not more than one active site in the system by adding a 
term 

Hchecku = E E 1^) (^1 ® 1^) (^1 ■ (4-124) 
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In fact, this is not sufficient to rule out that the two active sites wouldn't be at two places in 
the system that are spatially separated. However, Aharonov et al. [5j prove a clairvoyance 
lemma, which says that the expectation value of the complete Hamiltonian in a state with 
two active states is lower bounded by an inverse polynomial in n. This is because the two 
active sites have to move in the system, and thus will end up on neighboring sites. Also, 
one needs to check whether the number of qubits in the system is right. Whether there are 
not too few of them can be checked by 

Hf,^ = \B){Bl^ (4.125) 

on the n-th particle from the left. On the other hand, whether there are not too many of 
them can be checked by 

^„^an, = |®0)(®□lc„,c„+l■ (4-126) 

Finally, the initial state of the ancillae is checked by 

n 

k=m+l 

and the very first state in the sequence is initialized by 

Hciockinit = |m)(m|,^ + |©)(@|,^. (4.128) 

Altogether, for any quantum circuit U , Aharonov et. al. [5j constructed a Hamiltonian for 
d = 12 dimensional particles on a line with nearest neighbor interactions, whose ground state 
energy encodes whether the quantum circuit U can output 1 with high probability on some 
state. Thus, Quantum (12,12)-SAT on a line is QMAi complete and Local Hamiltonian for 
d = 12 particles on a line is QMA complete. 

4.3.2 Q-(11,11)-SAT on a line is QMAi complete. 

Here I present a modification of the above clock construction which requires only 11- 
dimensional qudits. The complexity of Q-(11,11)-SAT on the line is thus still QMAi com- 
plete (and Local Hamiltonian with d = 11 qudits on a line is QMA complete). This result 
is my own previously unpublished work. 

There are 8 dimensions required for the qubits in the previous clock construction. I 
needed two types of active qubits, one to send a message to the right @ and the other one 

to facilitate the application of gates. I also required two types of inactive qubits, one in 
the state 'done' (x) (to the left of the active spot), and the other one 'ready' Q (to the right 
of the active spot). This can be reduced by having only one type of inactive qubit Q, whose 
property ('ready' or 'done') is determined by the parity of its position. To compensate for 
this, I introduce two more non-qubit states ffl and B . The initial state is 

KKKKKK (g) BO BOB BBBBBBBBBBBB (4.129) 

^ V ' 

3 qubits separated by B 

where I now follow each qubit by the non-qubit state B. The progression of the states \^t) 
becomes somewhat more involved. To move the qubits two spaces to the right, one goes 
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through 



K K K K K ^(gBOBOB H H H H H H H H H H H H (4.130) 

K K K K K K ffl(g)OBOB H H H H H H H H H H H H (4.131) 

KKKKKKfflQ^BOBHHHHHHHHHHHH (4.132) 

K K K K K K fflQS(g)QB □□□□□□□□□□□ □ (4.133) 

KKKKKKBOfflO®BHHHHHHHHHHHH (4.134) 

K K K K K K BOfflOffl® □ H H H H H H H H H H H (4.135) 

K K ^ K K K BOfflOfflQO B B B B B B B B B B B (4.136) 

K K ^ K K K BQBOfflOfflO B B B B B B B B B B (4.137) 

KKKKKKfflQBOfflOffl<BBBBBBBBBB (4.138) 

K K K K K K fflOfflOBO < B H H H H H H H H H H (4.139) 

KKKKKKfflOfflOB<OBBBBBBBBBBB (4.140) 

^ K K K K K BQBQ < BOB B B B B B B B B B B (4.141) 

^ ^ ^ K ^ K BQB <]QBOB B B B B B B B B B B (4.142) 

K K K K K K BQ < BQBQB B B B B B B B B B B (4.143) 

K K K K K K ffl <OBOBOB B B B B B B B B B B (4.144) 

KKKKKK<BOBOBOBBBBBBBBBBB (4.145) 

K K K K K BQBQBOB B B B B B B B B B B (4.146) 

KKKKKKKOOBOBOBBBBBBBBBBB (4.147) 

K K K K K K K K@BOBOB B B B B B B B B B B (4.148) 



Observe that in ()4.148p . the qubits have moved two spots to the right with respect to their 
starting position (|4.130p . This is an analogue of obtaining the state ()4.92p from the state 
()4.73p in the previous Section. After the qubits have moved 2n spots to the right, the 
leftmost of them becomes again, and the next round of gate applications commences. 

The formal rules governing this progression of states are: 

la : ®B — ^ B® (4.149) 
lb : ©B — > B© (4.150) 

shifting the active qubit position. Next, when a qubit marked 'apply gates' meets another 
qubit, depending on their position, a unitary gate is applied to the internal states of the 
two qubits as the second qubit gets the 'apply gates' label. 

2a : ®0 ^ t/fc,„^(®®) (4.151) 
on particles (2kn + 2m, 2kn + 2m + 1), 

where k and m are integers. Meanwhile, the transition rule for a 'messenger' qubit meeting 
another qubit is just 

26 : ©0 ^ 0@. (4.152) 
When the qubit with the 'apply gates' label is at the right end of the qubit sequence, it 
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changes according to 



2a : ®Q — > (g)0 (4.153) 
on particles {kn, kn + 1), 

while everywhere else it changes according to 

2b : ©□ — y ®0. (4.154) 

The reason for this distinction between rules 2a and 2b is that the rules applied backwards 
also have to be unique. To send the active spot to the left now takes several steps: 

3a : OH — > SO (4.155) 
on particles {2k + 1,2k + 2), 

on neighboring particles with (odd, even) locations, while on the other pairs of particles the 
transition rule is 

36 : OH — > <B (4.156) 
on particles {2k + 2, 2A; + 3). 

Next, the qubits are pushed to the right as the active symbol <l moves to the left. Then 
the active spot bounces back as it reaches the left end of the qubit sequence. 

(4.157) 
(4.158) 
(4.159) 
(4.160) 

Finally, the leftmost qubit becomes active. When the particle pair involved is at position 
{kn, kn+ 1), it means that the qubits have moved 2n spots to the right and another round 
of gate applications can start. This activation proceeds according to the rule 

6a : CO — > (4.161) 
on particles {2kn — 1, 2kn), 

while everywhere else the last of the qubits is activated into the 'messenger' state @ as 

66 : — ^ ^© (4.162) 



4a 


ffl < — 


< B 


46 


o< - 


- <o 


5a 


K < — 


MO 


56 


B - 


MO 



4.4 Quantum 2-SAT in general geometry (for qudits) 
4.4.1 The Triangle Clock and Q-(5,3)-SAT 

In a recent paper, Eldar and Regev j^Hj, present a novel idea. Their 'triangle' clock con- 
struction abandons a simple linear progression of states |^'f) encoding the progression of 
the quantum circuit U . This allows them to prove that certain variant of Quantum 2-SAT 
for higher spins is QMAi complete. 

The complexity of classical 2- Satisfiability with higher dimensional particles is known. 
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apply gate U to 
target qubit 




applies if the control 
qubit f/j is 1 

applies if the control 
qubit is 



Figure 4-8: The 'triangle' clock transition rules allowing to apply a controlled gate to two 
work qubits using interactions that involve only one work qubit at a time. 



Already the simplest variant, (3,2)-SAT, i.e. 2-SAT for a trit and a bit is NP-complete, as 
shown in Section 11.1.31 On the other hand, the complexity of Q-SAT for higher spins is an 
open question. In Section 14.21 1 proved that Quantum (3,2,2)-SAT, i.e. Quantum 3-SAT 
for a qutrit and two qubits, is QMAi complete. Here I focus on the higher spin version of 
Quantum 2-SAT. Eldar and Regev's result which I now review is that Quantum (5,3)-SAT 
is QMAi-complete. Here (5, 3) means that each projector in their Hamiltonian involves one 
cinquit (a particle with dimension d = 5) and one qutrit {d = 3). 

The proof of QMAi-hardness relies on the 'triangle' clock construction depicted in Figure 
The state 



(4.163) 



has an allowed transition to the state 

\%+i) = {4'"Ui^t))<S)\t + i),, 



(4.164) 



which has the operator ax applied to a target work qubit q2- Next, the transition from 
the state |^'() to the state |^'t+2) is conditioned on a control (work) qubit qi to be in the 
state |0). On the other hand, if the control qubit qi is in the state |1), the transition 
l^'f+i) — > \^t+2) applies. In this fashion, a two-qubit controlled gate, in this case a 
CNOT, is applied to two work qubits, while interacting with only one of them at a time. 
The corresponding 'triangle' Hamiltonian consists of three terms: 



+ 



+ 



1 r. 



0" 



(92) 



t 

t+l,t 



(4.165) 



|0) (0| 



-^2 



1 r 
2 



Pt + Pt+2 - 
Pt+1 + Pt+2 



Xt+2,t + X 



t+2,t 



Xt+2,t+i + X^ 



+2,t+l 



The first term checks that for the transition between states t and t + 1, a ax gate is applied 
to the target qubit 52- The second and third terms check the transitions \^t) l^t+2) 
and l^'f+i) — > |^'f+2)i conditioned on the state of the control work qubit qi. The reader is 
welcome to check that the ground state energy of is exactly zero, and the states in the 
ground state subspace of iif^ are of the form 



Ca (• • • + l^t) + l^t+l) + |*t+2) + ■ • • ) 



(4.166) 
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where 

l^t) = l^)®|i)c> (4-167) 

l^m) = (ai^^^lV)) ®|t + l)e> (4.168) 

\%+2) = (CNOT,,,,, IV)) 55 |t + 2)^ (4.169) 

and is a state of the work register. 



The goal of Eldar and Regev in |28j is to construct a Quantum 2-SAT Hamiltonian (for 
a qutrit and cinquit), and for that one needs to implement the clock in a 2-local fashion. 
The clock register is made of alternating particles with dimension d = 3 and d = 5. The 
3-dimensional clock particle states are □ , ■ and M , while the 5-dimensional clock particle 
states are o , -^^ A , ► and x . The progression of clock states is similar to the domain 
wall combined with the pulse clock, as discussed in Section [4. 1.3t 



|l)c = 


■ 


□ 


o □ o □ o 


(4.170) 


|2). = 


M < 


□ 


o □ o □ o 


(4.171) 


|3). = 


M ▲ 


□ 


o □ o □ o 


(4.172) 


|4). = 


M ► 


□ 


o □ o □ o 


(4.173) 


= 


M X 


■ 


o □ o □ o 


(4.174) 


|6). = 


m X 




□ o □ o 


(4.175) 


= 


m X 




A □ o □ o 


(4.176) 


|8)c = 


m X 




► □ o □ o 


(4.177) 


|9). = 


M X 




X ■ o □ o 


(4.178) 


|io), = 


m X 




X m MB o 


(4.179) 


|ll)c = 


M X 




X M A B o 


(4.180) 


|12). = 


m X 




X Kl ► □ o 


(4.181) 


|13)c = 


M X 




X Kl X ■ o 


(4.182) 


|14). = 


m X 


m 


X m X m M 


(4.183) 



with a single active particle for a given state \t)^ (denoted by a black filled symbol). The 
three clock states \t)^ , \t + 1)^ and \t + 2)^ involved in the triangle construction in Figure 
14-81 correspond to three consecutive states of the clock register in whose the active site is 
on one 5-dimensional particle, such as |2)^ , |3)^ and |4)^. 



The propagation Hamiltonian Hprop in 14.61 now consists of four types of terms. For a 
neighboring clock qutrit and cinquit I have 

^^U(3,5) = ^(|")("|(3)®I(5)+I(3)®H)(^I(5)) (4.184) 
- ^ (|K ^) (■ o 1(3^5) + |Ho)(K ^1(3^5)), (4.185) 
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while for a neighboring clock cinquit and qutrit, the transition is checked by 



TTt 

prop{5,3) 



o l^)(^l(5)»%)+I(5)»l")( 



1(5) 



(5,3) 



□ )(X 



■1(3) 

1(5,3) 



(4.186) 
(4.187) 



The third type of term in Hp^op is the above described 'triangle' Hamiltonian (j4.165p . 
implemented using the following projectors and transition operators: 

(4.188) 
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This term checks the application of a controlled-not (CNOT) gate on two work qubits in a 
2-local fashion, as each of the three terms in ff^ is a projector acting nontrivially on the 
Hilbert space of one cinquit (5) and one qubit (2). 

So far, I have shown how Eldar and Regev check the application of a 2-local controlled 
gate on two work qubits by using projectors on one cinquit and one qutrit at a time. 
However, for universality one also needs to apply single qubit unitary gates. The final, 
fourth type of term in Hprop is then 



Tjt 

prop{single) 



(Pt + Pt 



t+ll 



+ I - ( Pt+l + Pt+2 



Xt+i,t -{U(S) Xt+i,t)^^ (4.189) 

(4.190) 



t+2,t+l 



+ x, 



t 

t+2,t+l 



and it is implemented using (j4.188p . The first line is a projector which checks the application 
of a unitary gate U on a single work qubit between states \^t) and |^'(+i). The second line 
is a projector checking a transition from |^t+i) to l^f+2) without changing the work qubits. 

Finally, the legal clock subspace checking operators are the analogue of the operators 
|01) (01| for the domain wall clock. They read 



Hclocki3,5) = \M){M\(^\o){o\ 

+ |H)(H|®(I-|o)(o|) 

+ |m)(Q|®(I-|o)(o|) 



(4.191) 
(4.192) 
(4.193) 



for a neighboring qutrit and cinquit, while for a neighboring cinquit and qutrit they are 



Hciock{5,3) = |x)(x|®|Q)(H| 

+ (I-|x)(x|)0|H)(H| 

+ (I- I x)(x|)0|K)(K| 

Whether there is an active site in the clock register at all is taken care of by 

Hclockinit = \a){^\c,+\^){^\ci^, 



(4.194) 
(4.195) 
(4.196) 



(4.197) 
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where ci and cl are the leftmost and the rightmost qutrits in the clock register. 

Eldar and Regev then proceed to show that if the circuit U accepts a state IV'j/es); 
the history state corresponding to IV'j/es) with extra terms for the states above the linear 
progression of states such as in Figure is the zero-energy ground state of H. Also, 

they show that in the 'no' case, the ground state energy of H is bounded from below by an 
inverse polynomial in n. Thus, Quantum (5,3)-SAT is QMAi-complete. In the following 
sections, I build on this result. 



4.4.2 A QMAi-complete 3-local Hamiltonian with restricted terms 

Here I implement the 2-local triangle clock which originally uses cinquits and qutrits, using 
qubits. This way I obtain another QMAi complete 3-local Hamiltonian construction which 
does not require any large penalty terms. Moreover, all the terms involved are constructed 
from a restricted set of operators. This is an unpublished result obtained with Peter Love 
during the Computational Complexity of Quantum Hamiltonian Systems workshop in Lei- 
den, Netherlands (July 2007). 

As in Section 14. 2^ I use a combined domain wall and pulse clock. First, replace each 
d = 3 particle in the construction of the previous section with two qubits u^i , W2 as 

□ — > 00, (4.198) 
■ — > 10, (4.199) 
K — > 11, (4.200) 

as in the domain wall clock in Section [4.1.31 Next, replace each d = 5 particle in the above 
construction with three qubits pi,P2,P3 of a pulse clock as 

o — > 000, (4.201) 

< — > 100, (4.202) 

A — > 010, (4.203) 

► — > 001, (4.204) 

X — > 000. (4.205) 

Note that the transitions between the active (pulse) states of the d = 5 clock involve only 
two of the three qubits pi. Therefore, the clock register transitions between \t)^, \t + 1)^ and 
\t + 2)^ in the triangle construction of the previous section are now implemented 2-locally 
as 



< > A 

< > ► 

▲ > ► 
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|01 




10 




|01 




10 




|01 





(4-206) 
(4-207) 
(4-208) 



Adding an interaction with a single work qubit at a time, this becomes 3-local. The propa- 
gation Hamiltonian is thus made from projectors on three qubits. This is true also for the 
transition from an active state of a qutrit clock particle to an active state of the neighboring 
cinquit clock particle. It is implemented 3-locally as 

mo ^ m< : |lOL„^JO)p^ |ll)^^^^Jl)p^, (4.209) 
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while moving the active state from a cinquit clock particle to a qutrit is done as 

►Q— X" ■■ |l)p3l00U,-2 ^ |0)p3|l0)-i,-2- (4-210) 



The legal clock state progression is now 
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(4.223) 



However, this is not a proof that Quantum 3-SAT is QMAi complete. I have not found a 
way to make the terms checking the legal clock states 3-local projectors. If I want the term 
which ensures proper coupling of the domain wall and the pulse clock to be a non-negative 
operator, it necessarily becomes 5-local. On the other hand, I have found a 5-local positive 
operator made out of 2-local operators, which does the job. This is reminiscent of the idea 
presented in Section I4.2[ 

The operators checking whether a state is in a legal clock subspace are the following. 
First, for the domain-wall clock particles, the operator 

Hciock(u,aii) = El01)(01|-.,-^ (4.224) 
{k,l} 

where w^jWi are consecutive domain wall qubits, checks the domain wall clock is correctly 
encoded. Next, I check whether the pulse clocks have at most one pulse in each of them by 

Hclock{pulse) = (|11) (ll|pW,p^'=) + |11) (ll|pW,pW + |11) (ll|pW,pW (4.225) 

k 

where Pi'\p2^\p^'^ is the k-th. pulse clock triplet. Finally, I need to check that the domain 
wall and the pulse clocks are correctly coupled. In other words, I need to check that the 
state in which the pulse clock was not initialized at the domain wall 

|...ll 000 00...) (4.226) 
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is illegal, while the states 
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(4.233) 



are perfectly fine. For this task, I now construct two versions of a 5-local positive semidefinite 
operator Hdockpuise, the first made from 3-local and the second made from 2-local terms. 
In both cases, Hdockpuise acts nontrivially on the space of three pulse qubits pi and two 
neighboring domain wall qubits Wi. 

In the first case, Hdockpuise uses at most 3-local interactions: 

It checks whether there is a single 1 in the pulse clock, when the neighboring domain wall 
qubits are 1 and 0. 

The second variant uses only 2-local interactions: 

4o"i;le = |l)(lL,®(cT^)+a^)+a^)+^^)-2) 

+ ^|ll)(ll|p^_p^,+ j;|l)(l|p^®|l)(lL^, (4.235) 

checking that when a domain wall qubit to the left of a pulse clock is 1, exactly one of the 
following 4 qubits (3 pulse and 1 domain wall) must be 1, as seen in (|4.229p - (|4.233p . 

Altogether, I constructed a Quantum 5-SAT Hamiltonian from at most 3-local interac- 
tion terms. Each of the 0(n) terms in this Hamiltonian has constant norm. Moreover, I 
retained the properties of the previous construction, as I only implemented the Hamiltonian 
(|4.6p in a different system. Therefore, Quantum 5-SAT made from 3-local interactions of 
a restricted type is QMAi complete. Moreover, 3-local Hamiltonian with constant norm 
terms of a restricted type is QMA-complete. 

4.5 Train Switch Construction: Universal Quantum Compu- 
tation using a Quantum 3-SAT Hamiltonian 

In this Section, building on Regev's triangle construction idea, I prove that one can perform 
universal quantum computation using a Quantum 3-SAT Hamiltonian H^s (for qubits). It 
is suitable for an Adiabatic Quantum Algorithm where I slowly change the Hamiltonian 
from a simple starting Hamiltonian Hb to -^35, or for a Hamiltonian Computer model 
where I let an easily prepared starting state evolve with H^s for not too long a time. 
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Figure 4-9: The 'pulse' clock construction consists of T + 1 clock qubits on a line, where 
only one of them is in the state |1) at a time. 
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Figure 4-10: The 'train switch' clock construction for the application of a CNOT gate using 
only 3-local terms. There are four extra qubits inserted between q and Q+i. The transitions 
from the original clock qubits to the gadget qubits are controlled by one of the work qubits. 
During the transition from ui to U2, the target work qubit is flipped. This gadget checks 
the application of a CNOT gate using 3-local projectors on qubits. 

As I want to perform a quantum computation, as opposed to constructing a Hamiltonian 
with a unique ground state, I have the advantage that I can choose the initial state of the 
system. This allows me to use the pulse clock (see Section 14.1.31 and Figure I4-9P 

\t)^=\0,,...Oct.^lctOct+^...Oc,) (4.236) 

as my starting point. The main innovation I add to it is the following train switch clock 
construction. My goal is to use a 3-local projector to check whether a 2-qubit controlled 
gate CNOTqj^qj properly applied to the work register when transitioning from the state 
l^f) to For this, I replace two clock qubits q and ct+i with a six-qubit gadget as 

in Figure 14-101 If the control qubit qi is |1), the active site in the clock register follows 
the upper rail, and takes the lower rail otherwise. Moreover, when the active site in the 
clock register moves from ui and U2, the target qubit q2 is flipped. This effectively applies a 
CNOT gate between qi and q2, between the states and {"^t+i)- In detail, the transition 
checking operators are 

Ht,u, = ^{l(S){Pt + Pu,)-\l){%^{Xu„t + xl^^)), (4.237) 

Huuu, = l{^^(.Pu,+Pu,)-cTi^'^®[Xu„u^+Xl^^^)), (4.238) 

Hu,,t+i = ^(l®(P„, + Pi+i)-|l)(l|g,55(Xi+i,„,+x/+,^„J) (4.239) 

for the upper rail, with the projectors P and time-increase operators X implemented as 
before for the pulse clock (see Section l4.1.3p : 

Pt = |l>(l|c.> (4-240) 
Xh,a = |10)(01|,_,. (4.241) 
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Analogously, for the lower rail I write 

Ht,di = l{^^{Pt + Pu,)-\0){0\^^(^[Xa,,t + xi;)), (4.242) 

Hd„d2 = l{^®{Pdr+Pd,)-l®{Xa„d,+xl,^)), (4.243) 
Hd„t+i = ^(l®(Pd,+Pt+i)-|0)(0|,^C5(Xi+i,d,+X^,J). (4.244) 
Let me examine the action of this Hamiltonian on the state 

l**> = l^i) ® \t)c = (« |0>,, |a> + & m<i2,..) ® \t)c • (4.245) 
When I consider only the forward moving terms in 

Htrain(t,t+1) = Ht,ui + -f^Mi,«2 + ^u^^t+l + Hf^di + -f^<ii,<i2 + Hd2,t+1 (4.246) 

(forgetting the identity and time-decreasing terms), I obtain the following three states: 

= a|0)^Ja)^^,...®|di), + 6|l)^J/3),,,...®K),, (4.247) 

1*^''^^)) = a|O)^Ja)^^,...0|d2), + Hl),J/3%2,...®Me, (4.248) 

l**+i> = («|0),J«),„... + 6|l),J/?%2,...)c^|t + l)c> (4.249) 

where 

W),,,..=4''¥),.,... (4.250) 

is the state obtained of the work qubits (besides qi) after the work qubit q2 was flipped. 
Observe that the work qubits of the state \^t+i) have the CNOT gate apphed to qi,q2, as 
I wanted, therefore 

\^t+i) = (CNOT,,,,jVt))®|i + l)c. (4.251) 

For a single qubit gate application, the transition checking is much simpler. At places 
where I do not insert the train switch gadget, I use the usual single qubit gate checking 
operator 

Hsingleit) = ^ (l ® (P* + Pt+l) -U® Xt+i,t - {U ^ ^t+l,t)^) . (4.252) 

The Hamiltonian 

H3S = ^ -f^smflie(t) + XI ^train{t,t+l)^ (4.253) 

corresponding to a quantum circuit ?7 is a Quantum 3-SAT Hamiltonian, because it is 

composed of 3-local projectors on qubits. 

The progression of states with the occasional intermediate states and 
\^^)f'^^'^^^), when a CNOT gate is on order, has the geometry of a line. Each of the states 
l^t) is connected by a transition only to two neighboring states. Let me call Tiiegal the 
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Figure 4-11: a) Computation on a line of states \^t)- b) The line corresponding to a circuit 
U padded with extra identity gates. The region where the computation is done is marked, 
c) Computation on a cycle of states \^t) coming from two copies of the computation on a 
line, d) The cycle corresponding to a padded circuit U. 



subspace of Ti spanned by the states Because I constructed it so, the Hamiltonian 

Hss (I4.253P does not induce transitions between the subspace Hiegai and 'H'^^gai- The time 
evolution of an initial state in TLiegai is thus governed solely by the restriction of H^s to 
TCiegai- Morcovcr, in the basis this restriction has the form 



H3S 



(4.254) 



same as Kitaev's propagation Hamiltonian (see Section 11.3.21 and Appendix [A]) . It is a 
Hamiltonian for a quantum walk on a line of length L' = poly{n), where n is the number of 
qubits the quantum circuit U acts on. Here L' = I+L1+3L2 where Li is the number of single 
qubit gates in U and L2 is the number of CNOT's in the circuit. To avoid complications 
in the analysis of the required running time coming from the endpoints, I can change the 
'line' of states |^f) into a circle of length 2L' as in Figure B-lll First, double the number 
of qubits in the clock register as 



|CoCi ...CL) 



|CoCl 



■ cl) Icqc'i ...4) , 



(4.255) 
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with the second 'hne' |cgc'^...c^) having analogous transition rules. Second, identify the 
endpoint qubits, i.e. cq = Cg and cl = c'^. This gives the clock register the geometry of a 
cycle, with a unique state |0)^ = IIcqO • • • 0) denoting time t = 0. The active site (spin up) 
in the clock register can proceed towards cl both ways, as in Figure 14-111 This new set of 
states \^t) with the geometry of a circle then defines the subspace "Hl^g^i- The Hamiltonian 
H-is restricted to this subspace is 



= \{l-B^), (4.256) 

legal 



where Bo is the adjacency matrix for a cycle of length 2L. The dynamics of this system is 
the quantum walk on a cycle, and I analyze it in detail in Appendix O There I show that 
when starting from the state on site A of the line and letting the system evolve for a time 
chosen uniformly at random between zero and a number not larger than O(Llog^L), the 
probability to find the state farther than L/3 from the starting point is close to |. This 
corresponds to finding a state with a spin up on clock qubit Ci-^l/^ (lower part of the circle) 
or Cti>L/^ (upper part of the circle). In both cases, because the circuit U was padded with 
identity gates, the state of the work qubits is then and contains the output of the 
quantum circuit U . Therefore, one can simulate a quantum computation U with Li single 
qubit gates and L2 CNOT's with this system as follows: 

1. Construct a clock register with C = 1 + Li + 4L2 qubits and transition rules as 
described by the train switch construction. 

2. Pad the clock register with extra 2C qubits (equivalent to padding the circuit U with 
identity gates). 

3. Double the clock register and make it into a cycle with length L = 6C. 

4. Initialize the system in the state 

1^0) = I00...0) 0|lcoO...O). (4.257) 



5. Let it evolve with H-^s (|4.253p for a time chosen uniformly at random between zero 
and r < 0(L log^ L). 

6. Measure the clock register qubits in the success area (farther than L/3 from the initial 
site, see also Figure [4-lip . With probability close to | you will find the active site 
(spin up) there. The work register now contains the output of the quantum circuit U . 
Restart otherwise. 

With this I have shown that Quantum 3-SAT Hamiltonians are powerful enough to 
perform universal quantum computation in the Hamiltonian Computer model based on a 
quantum walk. When the number of qubits and gates involved in the circuit is L, what I need 
for my Quantum 3-SAT Hamiltonian computer are thus 0{L) qubits, a Hamiltonian with 
0{L) projector terms with norm 0(1), resulting in \\H\\ = 0{L), and a running time of order 
0{L log^ L). The rescaled required resources (time x energy) scale as r- \\H\\ = ©(L^log^L). 
The initial part of this Hamiltonian, turning off the initial Hamiltonian pinning the initial 
state in {"^^q) can be turned off and H^g turned on adiabatically or quickly, as a quick change 
will not destroy the computation as discussed in Section 12.2.11 Moreover, Seth Lloyd [63] 
showed that this Hamiltonian Computer model is protected by an energy gap between the 
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legal subspace Hl^g^^i and the subspace orthogonal to it. Noise in the system will induce 
transitions away from 'H°gg„;, but their energy cost is going to be scaling like 0{L~^). The 
computation is thus protected by an energy barrier of the order 0{L~^). 
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Chapter 5 

Hamiltonian Quantum Cellular 
Automata 

Can universal quantum computation be performed by time evolving a simple system by a 
translationally invariant, time-independent Hamiltonian? I present ways how to do it in 
this chapter, largely based on the paper [TT], so far posted on the arXiv. 

Hamiltonian Quantum Cellular Automata in ID 

Daniel Nagaj, Pawel Wocjan 

We construct a simple translationally invariant, nearest-neighbor Hamil- 
tonian on a chain of 10-dimcnsional qudits that makes it possible to real- 
ize universal quantum computing without any external control during the 
computational process. We only require the ability to prepare an initial 
computational basis state which encodes both the quantum circuit and its 
input. The computational process is then carried out by the autonomous 
Hamiltonian time evolution. After a time polynomially long in the size of 
the quantum circuit has passed, the result of the computation is obtained 
with high probability by measuring a few qudits in the computational basis. 
This result also implies that there cannot exist efficient classical simulation 
methods for generic translationally invariant nearest-neighbor Hamiltoni- 
ans on qudit chains, unless quantum computers can be efficiently simu- 
lated by classical computers (or, put in complexity theoretic terms, unless 
BPP=BQP). 

This chapter is organized as follows. First, in Section [5.11 1 introduce classical cellular 
automata, quantum cellular automata and the Hamiltonian Quantum Cellular Automaton 
model. Then, in Section 15.2.11 1 present a HQCA in ID with cell size d = 10. I give 
another HQCA construction in ID with cell size d = 20 in Section 15.3.11 Throughout the 
required runtime analysis. Appendix [C] on the quantum walk on a line and Appendix [Dl 
about diffusion of free fermions on a line is referenced. 

5.1 Introduction 

One of the most important challenges in quantum information science is to identify quan- 
tum systems that can be controlled in such a way that they can be used to realize universal 
quantum computing. The quantum circuit model abstracts from the details of concrete 



131 



physical systems and states that the required elementary control operations are: (i) initial- 
ization in basis states, (ii) implementation of one and two-qubit gates, and (iii) measurement 
of single qubits in basis states. Meanwhile, many other models have been proposed such 
as measurement-based quantum computing [8ll [72l [60l [24] , adiabatic quantum computing 
[331 [7], or topological quantum computing [55l [79l [THl [66l EH HH HO] that reduce or modify 
the set of elementary control operations. However, the common principle underlying all 
these models is that the computation process is always driven by applying a sequence of 
control operations. 

Instead, I consider a model that does not require any control during the computational 
process. This model consists of a quantum system with a Hamiltonian that makes it possible 
to realize universal quantum computing by the following protocol: (1) prepare an initial 
state in the computational basis that encodes both the program and input, (2) let the 
Hamiltonian time evolution act undisturbed for a sufficiently long time, and (3) measure a 
small subsystem in the computational basis to obtain the result of the computation with high 
probability. I refer to this model as a Hamiltonian quantum computer and more specifically 
as a Hamiltonian quantum cellular automaton (HQCA) provided that the Hamiltonian acts 
on qudits that are arranged on some lattice, is invariant with respect to translations along 
the symmetry axis of the lattice, and contains only finite range interactions. Most natural 
Hamiltonians have these properties, so it is important to construct HQCA that are as close 
as possible to natural interactions. 

Hamiltonian QCA are related to the more usual discrete-time QCA (for further review 
of the different types of quantum cellular automata I refer the reader to |76] ) . However, 
while the evolution of discrete-time QCA proceeds in discrete update steps (corresponding to 
tensor products of local unitary operations, see e.g. [80^ 187]). the states of Hamiltonian QCA 
change in a continuous way according to the Schrodinger equation (with a time-independent 
Hamiltonian). For this reason, Hamiltonian QCA are also called continuous-time QCA [76] . 
Also, in the HQCA model, all the couplings (interactions) are present all the time, while for 
the the discrete-time QCA, the execution of updates on overlapping cells is synchronized by 
external control. Therefore, the nearest-neighbor interactions of a HQCA have to include 
a mechanism that ensures that the logical transformations are carried out in the correct 
order. 

The motivation to consider Hamiltonian computers is threefold. First, it is a fundamen- 
tal question in the thermodynamics of computation how to realize computational processes 
within a closed physical system. Such Hamiltonian computers were presented and discussed 
by Benioff [TB], Feynman [39|, and Margolus ^64j- Second, Hamiltonian quantum cellular au- 
tomata could lead to new ideas for reducing the set of necessary control operations in current 
proposals for quantum computing by using the inherent computational power of the inter- 
actions. HQCA are at one end of the spectrum of possible implementations; more realistic 
perspectives for quantum computing could arise by combining this model with more con- 
ventional models involving external control operations throughout the computation. Third, 
this model can show the limitations of current and future methods in condensed matter 
physics for simulating the time evolution of translationally invariant systems. If evolving 
with a certain Hamiltonian can realize universal quantum computing, then there cannot 
exist any classical method for efficiently simulating the corresponding time evolution unless 
classical computers are as powerful as quantum computers (BPP=BQP). 

The first theoretical computational models based on a single time-independent Hamil- 
tonian go back to [m [39l |6l] . However, these Hamiltonian computers were not explicitly 
designed for realizing universal quantum computing. Margolus' model jMj has the attrac- 
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tive feature that it is laid out on a 2-dimensional lattice with translationally invariant, 
finite-range interactions. (In [18] it was argued that the part of the Hamiltonian respon- 
sible for the synchronization in a 1-dimensional variant is close to real interaction in sohd 
states.) However, this scheme does not satisfy the requirement (1) since its initial state has 
to be prepared in a superposition. Building upon Margolus' idea, a translationally invariant 
Hamiltonian universal for quantum computing even if the initial state is restricted to be 
a canonical basis state was given in 08]. This model requires 10- local, finite-range inter- 
actions among qubits on a 2-dimensional rectangular lattice wrapped around a cylinder. 
Subsequently, it was established in [U] that nearest-neighbor interactions among qutrits on 
a 2-dimensional lattice suffice. However, the Hamiltonian of [47] is translationally invariant 
only when translated over several lattice sites. A different approach was taken by Vollbrecht 
and Cirac in |107j , showing that one can implement universal quantum computation with a 
translationally invariant, nearest-neighbor Hamiltonian on a chain of 30-dimensional qudits. 

I present two different simplified HQCA constructions on one-dimensional qudit chains. 
In both models, I think of the qudit chain as composed of two registers, data and program. 
The work qubits I compute on are located at a static location in the data register. Driven 
by the autonomous Hamiltonian time evolution, the program sequence contained in the 
program register moves past the work qubits and the gates are applied to them. After I 
let the system evolve for a time not larger than a polynomial in the length of the program, 
I measure one or two qudits in the computational basis to read out the output of the 
computation with high probability. 

My first construction is for a chain of 10-dimensional qudits and is related to the ideas 
of |107j . The mechanism behind the progress of the program sequence in this particular 
model can be thought of as the diffusion of a system of free fermions on a line. My second 
construction uses qudits with dimension d = 20 and is inspired by [49] , utilizing a technique 
of [5] to transport the program. Here, the mechanism for the progress of the computation 
can be thought of as a quantum walk on a line. 



5.2 The HQCA in ID with cell size d = 10 

I present a simple universal HQCA on a chain of qudits with dimension d = 10. First, I 
encode the progression of a quantum circuit U on N qubits into a set of states {ipa) of a 
chain of qudits with length L = poly{N). Second, I give a translationally invariant nearest- 
neighbor Hamiltonian on this chain of qudits, which induces a quantum walk on the set 
of states Ivo-)- Finally, using a mapping to a system of free fermions in ID, I prove that 
when I initialize the qudit chain in an easily determined computational basis state and let 
the system evolve for a time r < rio = O(-Llog-L) chosen uniformly at random, I can read 
out the result of the quantum circuit U with probability piQ > | — O {j^^Z^ t)y measuring 
one of the qudits in the computational basis. I then show that this is enough to ensure 
universality of my HQCA for the class BQP. 

5.2.1 The Construction 

The gate set {Toffoli, Hadamard} is universal for quantum computation [88j. With only 
polynomial overhead, one can simulate a circuit consisting of these gates using only the gate 
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Figure 5-1: a) A quantum circuit consisting of two sequences of gates acting on nearest 
neighbors, b) The previous circuit with a third sequence of identity gates added. 
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(5.1) 



if it can be apphed to any pair of qubits. Let me consider implementing universal quantum 
computation on a qubit chain using only nearest neighbor gates. Let me also restrict the 
use of the W gate so that the control qubit has to be to the left of the target qubit. 
Using only polynomially many additional swap gates S, one can still do universal quantum 
computation on a qubit chain. Thus given a quantum circuit U' on N' qubits with poly{N') 
generic two-qubit gates, I can transform it into a circuit C/ on a chain of = poly{N') 
qubits with nearest neighbor gates W (with control on the left) and S without loss of 
universality. I then add identity to my gate set and further transform the circuit U to have 
the following form (see Figure [5^ . Rewrite the circuit as K sequences of nearest neighbor 
gates Uk,g G {W,S,I}, where gate Uk^g belongs to the k-th sequence and acts on the pair 
of qubits {g, g + 1): 



U = {Uk,n-i ■ . . Uk,i) ■ ■ ■ (C/i,7V-i . . . f/i,i)- 



(5.2) 



I wish to encode the progression of the circuit U into the states of a chain of qudits with 
dimension 10, with length L = poly{N). The basis states of each qudit \q) = \p) ® \d) are 
constructed as a tensor product of a 5-dimensional program register p € { . , 5, /} 

and a 2-dimensional data register d € {0, 1}. I start by writing the initial product state 

L 

W) = ®{\P,)®\d3)),, (5.3) 
i=i 

with pj and dj as follows (here I give an example for the circuit in Figure 15-lb ) : 
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The qudit chain has length 

L = 2M = 2KN, (5.5) 
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where K is the number of gate sequences in (j5.2p . The left half of the top (program) 
register contains K pointer symbols ► at positions kN for k = 1 . . . K and empty symbols 
. everywhere else. The right half holds the program in the form 

/ . . . Ui,N-i I U2,i . . . U2,N-i I...IUk,i... Uk,n-i , (5.6) 

V V V 

gate sequence 1 gate sequence 2 last sequence 

with C/fc g G {VF, 5, /} and each sequence preceded by an identity gate. The bottom (data) 
register contains A'^ work qubits (labeled Wn in the table) at positions M + n for n = \ . . . N 
and qubits in the state |0) everywhere else. I designate as the readout qubit. 

I now give two simple rules, many applications of which generate the set of states {| Vo-)} 
from the initial state |(/?). I label the states by cr, the description of the sequence of rules 
I choose to apply to the initial state. I invite the reader to work out what happens to \ip) 
(|5.4p as the rules are applied, noting that there are usually several possible rules that can 
be applied to a given state. The first rule says that the symbols A € {W, 5, 1} (from now 
on I call them gates and think of them as particles) in the program register can move one 
step to the left, if there is an empty spot there: 
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The second rule concerns what happens when a gate meets a pointer symbol: 



► 


A 


X 


y 



A 


► 


A{x,y) 



(5.7) 



(5.8) 



In this case, the gate moves to the left, the pointer ► gets pushed to the right, and the 
gate As {W,S,I} is applied to the qubits in the data register below. 

The initial state has K pointers, one for each sequence of gates in the circuit. I con- 
structed the initial state (|5.4p in such a way that as a gate Uk,g from the /c-th gate sequence 
moves to the left, it meets the k-ih. (counting from the right) pointer ► exactly above 
the work qubits to which the gate was intended to be applied {g,g + 1). However, one also 
needs to consider what happens when a gate meets a pointer for a different sequence. If this 
happens over a pair of extra qubits in the state |0) |0), the qubits stay unchanged because 
the gate is either the controlled gate W , a swap gate or the identity. The second possibility 
is that a gate meets a pointer above the boundary of the work qubits (i.e. |0)^ \ m+i 
or \wn) M+N \^) M+N+i)- I ensured that it is going to be the identity gate by inserting one 
in front of each gate sequence in ()5.6p . This implies that the leftmost (or the rightmost) 
work qubit again stays unchanged. Therefore, I am sure that the state |(/?o-) in which the 
gate particles have all moved to the left half of the chain contains the result of the quantum 
computation U in the state of the work qubits, while the additional qubits in the data 
register remain in the state |0). 

Let me now allow all the rules to be applied backwards as well, opening the possibility of 
returning back to the initial state \ip) (undoing the computation). Although the number of 
possible sequences of rule (and backward rule) applications then becomes infinite, the space 
of states {|</?o-)} is nevertheless finite-dimensional. In every state Ivo-), the M gates {W, S, 1} 
in the program register that started at positions {M + 1, . . . , 2M} occupy some combination 

C = {a\ , . . . , 0,)^ } of the L sites of the chain. Given the rules ()5.7p - (|5.8p . the order of 

(C) (C) 

the gates in the program register cannot change, i.e. a) < am ior k < m. Because the 
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positions of the pointers ► and the state of the data register are also uniquely determined 
by {a^f'\ . . . ,a^^^}, I can label the states I constructed by \(pc) with C = 1 . . . (^). 

I now give the universal translationally invariant Hamiltonian as a sum of translationally 
invariant terms 

«'»=-|:(^+^*)c„«,' p-^' 

where R corresponds to the rules (j5.7p - (j5.8p and acts on two neighboring qudits as 

^ = Yl [M ■ ) ( ■ ^ \p,,p, ® ^dud2 + I ^ ►) ( ► ^ \p„p, ® M,d2] > (5-10) 

Ae{W,S,I} 

where p stands for the program register and d for the data register of the respective qudit. 

Recall that my model of computation consists of initializing the qudit chain in the state 
\ip) and evolving the system for a time r < tiq which I will determine in Section 15.2.21 
Finally, I want to show that when I measure the output qubit wn in the data register, I 
will read out the result of the quantum computation U with high probability. 



5.2.2 Required Time Analysis 



The time evolved state l^pir)) (obtained from \(p) by evolving with Hiq for time r) is a 
superposition of the states {ipc)- I can write it as 



l^i^)) = X/ V'ciT) \Trc)program \dc)work \o!)extra, 
C=l 



(5.11) 



where \ttc) is the state of the program register of the chain in the state \fc)j the corre- 
sponding state of the work qubits is \9c) and \a) is the state (all zero) of the extra data 
qubits. The state \6c) of the work qubits holds the output state of the computation U, if 
all of the gate "particles" have moved to the left of the work qubits. Let me now choose 
some number / and pad the qubit chain with (/ — 1)M empty sites on the left and M sites 
containing / in the program register on the right: 



pad\ 







3(/-l)Af 



(8) \f) 



/ 





(5.12) 



The original chain had length 2M, so the length of this padded chain is L = (/ + 2)M. 

The program register of the initial state |(/?^'^'^) = \(pc=i) now has 2M gate "particles" 

{W, S, 1} at positions = fM + m with m = 1 . . . 2M. On this modified chain, every 

state l^pc) in which the first M gate "particles" are located in the first fM sites of the chain 
(C) 

(om < fM for m < M) contains the finished computation in the state of its work qubits 
\9c)- Note that for all these states, the state of the work qubits is the same and equal to 
\du) = U {wi . . . wn), as it does not change under the extra identity gates I added. I can 
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now rewrite the time evolved state (jS.lip as 



V'c) I ~l~ \'^)prog ® \6u)work ^ \o!)extra, (5.13) 
VM>f^ ) 

where the sum in the first term is over the set of positions of the 2M gate particles in which 
the M-th particle is still near the right end of the chain and the computation is thus not 
finished yet. Meanwhile, 

\T^)prog = ^ 'fci'T) \T^c)prog- (5-14) 

is a superposition of the program register states which correspond to an executed compu- 
tation. 

Recall the definition of the quantum complexity class BQP (Definition [6] in Section 
ll.2.2p . There is a verifier circuit U associated with each problem. Denote x a problem 
instance. As the problems in BQP are decision problems, x G {Lyes U Lno), where Lyes is 
the set of instances with the answer 'yes'. Now, when x £ Lyes, the probability of the circuit 
U outputting 'yes' is not smaller than |. On the other hand, when x G L„o, the probability 
of the circuit outputting 'yes' is not greater than ^. Let me assume the worst case for 
the circuit U, i.e. that the circuit outputs 'yes' on a good proof of a 'yes' instance with 
probability pu = ^- In the language of spins, the circuit U outputs 'yes' when I measure 
spin up on the output qubit. Therefore, the expected value of measuring az on output qubit 
of the circuit U is bounded from below by 

(^a)S^r* >lxpu + {-l)x{l-pu) = 2pu-l = l (5.15) 
when X G Lyes- Analogously, when x G Lno, it is bounded from above by 

{<rii)Z''''' < -2pu + I = -I- (5.16) 



To solve BQP problems with my automaton, I need to distinguish the 'yes' from the 
'no' cases, i.e. I need to show that the expectation value of measuring a^^^ on the output 
qubit of my automaton at a random time r < rio is greater than zero in the 'yes' case, 
and smaller than zero in the 'no' case. First, call pio the probability to find a state where 
the computation is finished. Let me consider a yes instance (x G Lyes)- Recalling (j5.13p 
and observing that \Tr)prog is orthogonal to the states of the program register in which the 
computation is not finished, I have 

{<TLi)yes = (V^(r)l 41 \^ir)) = p,o {9u\ 4^ \(^u)^ork +(1 - Pio) W\ 41 W) (5.17) 

^ V ' 

output of U 

where is the normalized first term in ()5.13p . The first term is the circuit output (|5.15p . 
therefore 

{Ou\ crLi = (41)™* > 2pu - 1. (5.18) 
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The second term can be bounded from below (adversarially, i.e. for every time the compu- 
tation is not finished, the output qubit gives the opposite of the correct answer) by 

(v^Vit k'> > -1- (5-19) 
Putting it together, the expectation value is bounded by 

> pio(2pc/-l)-(l-Pio) = 2pioPr/-l. (5.20) 

Analogously, for the x ^ L case, I obtain 

{<y^^l)no > -2pwpu + l. (5.21) 

I wiU now prove that when I choose the time r uniformly at random in (0,tio), with 
Tio = poly{M), the probability of finding a state with the computation executed (with 

aM < fM) is pio > I - O (4) with L = (/ + 2)M. 

Let me analyze the time evolution of \ipP°''^) (j5.12p under Hio (j5.9p . I can restrict the 
analysis to the program register of the chain, as the content of the data register in the 
time-evolved state \^{t)) is completely determined by the content of the program register. 
The data register does not hinder the time evolution of the program register in any way. 
In fact, there exist bases, in which Hiq is identical (as a matrix) to Hiq restricted to the 
program register. Moreover, let me consider a further mapping of the system restricted to 
the program register to a line of qubits with length L = (f + 2)M as follows. Map the 
states {►, . } to the state |0), and the states {W,S,I} to the state |1). The mapping of 
HiQ to this system is a sum of hopping terms 

L-1 

H, = - ^(|10) (Oil + |01) (10|),-,+i. (5.22) 
i=i 

Using the Wigner-Jordan transformation I can define the operators 

b] = a^..a|_l®|l)(0|^.®Vl,..„i, (5.23) 
bj = aJ...cT|_i®|O)(l|^.0%+i,...,L. (5.24) 

As 6j and bj have the required properties {bi,bj} = 6ijl and b'j = = 0, they can be 
viewed as the creation and annihilation operators for a fermion at site j. Rewriting (|5.22p 
in terms of dOSjl - ilOi]) . I obtain 

L-1 

Hf = -Y^ b]bj+i + h.c, (5.25) 

i=i 

a Hamiltonian for a system of free fermions in second quantization. Following the mapping, 
the initial state \ip'P°''^) of the qudit chain thus corresponds to the state of the fermionic 
system |^') = b^jrj^,j_^i ■ ■ ■ b]^f^2)M 1^) '^ith 2M fermions on the right end of the line (here |0) 
is the state with no fermions). I now use the following Lemma proved in Appendix [Pl 

Lemma 3. Consider the state \^) of 2M fermions on the right end of a line with L = 
{f + 2)M sites. Let the system evolve for a time chosen uniformly at random between and 
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Tio with the Hamiltonian given by (j5.25p and measure the number of fermions in the region 
1 <^ 2; < fM. The probability to measure a number greater than M is pio > jq^ — O 

Let me choose / = 22 and tiq = 0(L log L) = 0{M log M). Following the mapping I 
did from my qudit chain backwards, this implies that when I initialize the qudit chain of 
length L = 2AM in \ip'P'"^) as in (|5.12|) and let it evolve with Hiq (j5.9|) for a random time 
■7" ^ ''"10) the probability for the chain to be in a state where the gate particles have moved 
sufficiently to the left for the computation to be done {om < f^) is 

Therefore, equations (j5.20p and (j5.2ip now read 

1 „ / 1 



logM 



Therefore, one can recognize any language in BQP using the HQCA I described above. 

As an aside, note that there is a way to determine that I obtained a state in which 
the computation has been done with certainty (and thus getting rid of the second term in 
(|5.17p ). I could have chosen to measure all the program qudits to the right of the first work 
qubit and check whether all the S and W are gone. This happens with the above probability 
pio, and the postselected state of the work qubits now surely contains the output of the 
circuit U. Note also that I can think of the state of all the work qubits as the circuit output, 
as compared to only the last work qubit. Nevertheless, thinking only about the last work 
qubit is enough to ensure universality of my HQCA for the class BQP. 



5.3 The HQCA in ID with cell size d = 20 

I now present the second construction, a HQCA for a chain of 20-dimensional qudits. As 
in Section 15.2. H I describe an encoding of the progression of a quantum circuit U into a 
set of states of a qudit chain. However, the geometry of this set of states j-i/^j) will be now 
much simpler, as I can label them by a "time" label t = 1 . . . L = poly{N), thinking of the 
set of states as a "line" . The Hamiltonian i?20 I construct induces a quantum walk on this 
"line" of states. I conclude by proving that when I let the initial state I'i/'o) evolve with 
for a time r chosen uniformly at random between and T20 = O(ilogL), I can read out 
the result of the quantum computation U with probability p20 ^ measuring 
two of the qudits in the computational basis. 



5.3.1 The Construction 

I encode the progression of a quantum circuit U in the form ()5.2p (see also Figure 15-ip into 
a set of states \tpt) of a qudit chain with length L = {2K — 1)(A^ + 1) + 2. As in Section 
15.2. 1^ each qudit consists of a program register and a data register. The data register is 
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again two-dimensional, but the program register can now be in the following 10 states: 



W, S, I 



► 



the program sequence, 

marked characters in the program sequence, used to propagate 

the active spot to the front (left) of the program sequence, 

apply gate symbol, 

shift program forward, 

a turn-around symbol, 

empty spot (before/after the program). 



Similarly to ()5.3p and (15. 4p . the initial product state iV'o) = 
by (the following is an example for the circuit in Figure [5^Tb ) 



(8)t=i {\Pj) ®\d3))j is given 
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(5.28) 



In general, the data register contains N work qubits (labeled Wn in my example) at positions 
(i^— l)(A''-|-l)+2-|-n for n = 1 : (counting from the left). Qubit wjy is the designed output 
qubit for the computation, i.e. once the computation is done, wn contains the output of U. 
Next, the data register contains qubits in the state |1) at positions (A; — 1)(A^ -|- 1) -|- 2 for 
k = 1 . . . 2K and qubits in the state |0) everywhere else. The I's serve as sequence boundary 
markers. The program register has empty symbols . on the left, and then it contains the 
program in the form 



/ C/i,i . . . UiM-iII U2A . . . U2N-1II ■■■II Uka ■ . . Ukm-u 



(5.29) 



1^' gate sequence 



gate sequence 



last gate sequence 



with the program written from left to right. In the example given in (j5.28p . the first gate 
sequence (see Figure 15-lb ) is WS and the second gate sequence is SW. Finally, the last 
qudit in the program register is in the state O , marking an active spot in the computation. 

I now give the rules to obtain the sequence of states {ipt) from iV'o)- These rules are 
constructed so that there is always only one of them that can be applied to a given state 
\^pt), thus giving me a unique state iV't+i)- (Also, using the rules backwards, one obtains a 
unique \ipt-i) from {ipt))- The first three are 



1 : 

2 : 

3 : 
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(5.30) 



O A 



where A, B stands for either W, S or I. These rules ensure the passing of the active spot 
from the back end (right side) of the program to the front (left side), without modifying 
the data register or the order of the gates in the program sequence. Next, 
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(5.31) 
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After the active spot has moved to the front of the program, there are two possibihties. The 
turn symbol O can change to the apply gate symbol ► (rule 4a), or to the shift program 
symbol > (rule 4b), depending on whether the data qubit below contains the sequence 
boundary marker state 1. Afterwards, for the states containing the apply gate symbol ►, 
I have: 



5a : 
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y 
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A{x,y) 



6a : 



► 
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1 





(5.32) 



When applying rule 5a, the apply gate symbol ► moves to the right, while a gate from the 
program sequence is applied to the qubits in the data register below. Applying the rule 
repeatedly, the ► symbol moves to the right end of the program sequence. As an example, 
I now write out the state \1p12) that I obtained from the state {ipo) applying rules 1, 2 (6 
times), 3, 4a and 5a (3 times) from the state |^o)- 
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where | . . . ^ . . . ) stands for the state of the three work qubits after the gates 1^12 and then 
S23 were applied to them. Let me now take a closer look at the marker qubits (all qubits in 
the data register except for the work qubits Wn) and the application of rule 5a. The marker 
qubits stay unchanged for all l-^i). The gate applied to pairs |0) |1) and |1) |0) of marker 
qubits or the pairs of qubits |1) \qi) and {qn) |1) (the left and right ends of the work qubit 
sequence) is always I, because of the identity gates I inserted between sequences of gates 
in the program ()5.29p . Finally, the qubit pairs |0) |0) between the 1 markers do not change 
under the swap operation or the W gate (a controlled gate). 

After the apply gate ► symbol gets to the end of the sequence, it changes into the turn 
symbol O via rule 6a. Note that the boundary markers in the data register are spaced in 
such a way, that the ► symbol will arrive at the right end of the sequence when the qubit 
below is in the state 1. Note that at the very end of the line, rule 6a needs to be modified to 
involve only the two particles directly above each other. Using rule 6a, ► will then change 
into the turn symbol O . After applying rules 1, 2 (6 times) and 3, the active spot again 
moves to the left of the program. Because the O symbol is now above a marker qubit, 
rule 4b can be used, and I get a state with the shift program symbol O. Finally, here are 
the last two rules: 



56 : 



> 


A 




A 


> 



66 : 



> 





















(5.34) 



where again A stands for either W, S or I . Rule 5b makes the program shift to the left 
while the > symbol moves to the right. Finally, rule 6b deals with what happens when the 
> symbol arrives at the end of the program sequence. Because of the way I constructed 
the data register in iV'o)) the data qubit below the > symbol will then be in the state 0, so 
that the [> symbol changes to the turn symbol O . The reason why I need to look at the 
qubit in the data register below the O symbol in rules 6a and 6b is that when I apply the 
rules backwards (making \ipt-i) from \ipt)), again only one of them applies for each \ipt)- 

After applying rule 1, 2a (6 times), 3 and 4b, the \> symbol appears again and starts 
shifting the program further to the left. After several rounds of this, when the program 
shifts to the left by + 1, rule 4a can be used again (as the O symbol will be above a 
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1 marker qubit), and subsequently, the ► symbol facilitates the application of the second 
sequence of gates to the work qubits. 



After many applications of the above rules, I arrive at the state for which none of 
the (forward) rules apply. 

, , ^ I O I W S I I S W 1 o.^ 

1"^^^ " [ 1 1 \...e'...) 1 1 J- 

This is the state in which the program has moved to the left of the qudit chain, and all 
sequences of gates have been applied to the qubits in the data register. The state \ . . .9' . . .) 
is thus the output state of the circuit U and the last of the work qubits (wn) holds the 
output of the quantum computation. 



Starting from ()5.28|) . I have constructed the set of states \tpt) for t = . . . L with L = 
0{K'^N'^) = poly{N). As t grows, these states encode the progress of a quantum circuit U. 
What is the geometry of this set of states? They are labeled by a discrete label t, with the 
state I'i/'t) obtainable only from the states \ipt-i) and IV't+i) using the above rules and their 
backward applications. Therefore, the states \^pt) can be thought of as position basis states 
on a line of length L + 1 

li^t) ^ \t)une^ (5-36) 

where t = . . . L. 



Let me choose a Hamiltonian H20 for this system as a sum of translationally invariant 
terms: 

L-l 66 

'^- = -EE(n. + ^i),._.^„ (5.37) 

i=l k=l 

where the terms correspond to the rules l-6b (|5.30p . ()5.3ip . (|5.32p and (|5.34p and act on 
two neighboring qudits as 

^'i = El® •)^^ 0\j,^^p^^Id,42, (5.38) 
Ae{w,s,i} 

P2 = Yl (5.39) 

A,Be{W,S,I} 

P3 = E 10 A)( .®\p„p,^hud„ (5.40) 

A£{W,S,I} 
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and 



P4a = I . ►)( . 0|p^,p,®lrf, (5.41) 

P46 = I . >> ( . O ® Id, ® |0) (0|,^ , (5.42) 

Psa = Yl M ►)(► ^lpi,P2®^rfi,'^2, (5-43) 
A£{W,S,I} 

P5b = Yl \^ ^\pi,P2'^^d,,d2, (5-44) 

A£{W,S,I} 

Pea = 10 .)(► ■\p,,p,^\^){Md,^^d2, (5.45) 

Peb = 10 .)(> ■U,p,O|0)(0|,^®I,,. (5.46) 

When thinking of the set of states \ipt) as the set of positions of a particle on a hne (j5.36p . 
H20 becomes 

L-l 

Hiine = - X] (I*) + 1| + I* + 1) (*l ) • (5-47) 

This is the Hamiltonian of a (continuous-time) quantum walk on a line of length L + 1. 
Therefore, H20 induces a quantum walk on the "line" of states \ipt)- 



5.3.2 Required Evolution Time Analysis 

After initializing the qudit chain in the state |^o) fmd evolving with H20 for time r, the 
final step is to read out the output of the computation. As in Section [5.2.21 I need to ensure 
that the probability of finding the chain of qudits in a state where the computation was 
performed completely is high. To raise this probability, I choose to pad the program (K 
sequences of gates) with another 5K sequences of identity gates and redo the construction in 
the previous section. The length of the qudit chain thus becomes L = (2(6i^) — l)(A^+l)-|-2. 
The states \ipt>L/6) (with L modified) now all contain the result of the quantum circuit U 
in the readout qubit wat, as the relevant gates have been applied to the work qubits in those 
states. Note that as the extra identity gates pass by, the state of the work qubits does not 
change. 

The readout procedure consists of two steps. First, measure the qudit Pl-k{n+i) 
the program register (the qudit with distance from the right end of the chain equal to the 
length of the original program) . Let me call p2o the probability to measure . (which would 
mean the program has moved to the left of the qudit I just measured). When this happens, 
I am assured that the work register is in a state in which the computation is done. Second, 
I measure wn, the last of the work qubits, and read out the result of the computation U. 
I now prove that when I choose to measure Pl-k{n+i) at a random time < r < T20 with 
T20 = poly{N), the probability p2o of obtaining the state . is close to |. 

To simplify the notation, let me label the states as \t). In this basis, the Hamiltonian 
(j5.37p is the negative of the adjacency matrix of a line graph with L + 1 nodes. I now use 
the following lemma about the quantum walk on a line proved in Appendix O 

Lemma 4. Consider a continuous time quantum walk on a line of length L + 1, where the 
Hamiltonian is the negative of the adjacency matrix for the line. Let the system evolve for 
a time r chosen uniformly at random between and T20, starting in a position basis state 
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|c). The probability to measure a state \t) with t > L/6 is then p2o > | — O ("^^i^j ' 

This implies that when I initiahze the qudit chain in the state iV'o) (corresponding to 
the leftmost state on the line |c) = |1)) and let it evolve with H for a random time r < T20 
with T20 = 0(L log L), the probability to find a state with t > L/6 is close to |. Therefore, 
when I measure the program qudit Pl-K{n+i)^ I will obtain . with probability close to 
|. Finally, when I subsequently measure the work qubit wn, 1 will obtain the result of the 
quantum circuit U. 

Note that I can also avoid this postselection procedure and simply measure the output 
qubit. The analysis of the outcome would then follow what I did above in Section 15.2.2^ 
resulting in (j5.27p again, with M replaced by L. 
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Chapter 6 

Conclusions 



In this thesis I presented three ways to look at local Hamiltonians in quantum computation. 
Some of these Hamiltonians are useful for building a quantum computer, some of them can 
be simulated classically, while the properties of others are hard to determine even on a 
quantum computer. 

First, in Chapter [2] I showed that local structure is a necessity for Adiabatic Quan- 
tum Computing to work, and also that a simple 3-local Hamiltonian can be used in the 
AQC model to simulate any quantum circuit much more effectively than had been known. 
Second, in Chapter [3] I presented a numerical method to find the ground state of a trans- 
lationally invariant Hamiltonian with an infinite tree geometry of interactions, investigated 
two particular AQC Hamiltonian models, and found encouraging results about the phase 
transitions in these systems. Third, in Chapter [J] I obtained several new results about the 
complexity of Quantum /c-SAT and Local Hamiltonian, implying that finding the ground 
state properties of rather simple local Hamiltonians is unlikely to be possible even on a 
quantum computer. Finally, returning to the idea of using local Hamiltonians for build- 
ing a quantum computer, in Chapter [5] I presented two Hamiltonian Quantum Cellular 
Automata, showing a way to perform quantum computation in continuous time with a 
translationally invariant, time-independent Hamiltonian on a line, reminiscent of universal 
Turing Machines. 

Adiabatic Quantum Computing 

— how to use it and how to make it fail 

My results about how not to design quantum adiabatic algorithms are a response to re- 
searchers who made a bad choice of the initial Hamiltonian and then showed that the 
algorithm doesn't work. I showed that throwing away local structure in the Hamiltonian 
maps the AQC algorithm into running the unstructured search problem backwards. The 
running time of any such attempt must then take time of order at least \/2", dashing the 
hopes of obtaining an exponential speedup. On the other hand, it has not been ruled out 
that a properly chosen initial Hamiltonian could lead to algorithmic success. A recent nu- 
merical study of the NP-complete Exact Cover problem up to 128 bits shows promising, 
inverse-polynomial scaling of the eigenvalue gap governing the required runtime of the AQC 
algorithm |113j . Besides using AQC for solving locally constrained optimization problems, 
it remains an interesting alternative to the conventional quantum circuit model, because it 
is possible to simulate quantum circuits using AQC Hamiltonians. Moreover, I showed that 
the previously estimated required running time r oc L^^ (where L is the number of gates in 
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the circuit and the norm of the Hamiltonian is also proportional to L) for such simulations 
was grossly overestimated. When using an AQC Hamiltonian made from 3-local projectors 
I constructed in Section 14. 5^ the running time needs to scale only as r • \\H\\ oc L^log^L. 
In fact, I believe the log factors in the required running time can be omitted, as I required 
them to prove strong convergence of the underlying quantum walk (see Appendix [C]), while 
showing weak convergence would be enoughEl- Moreover, such computation is inherently 
protected against decoherence by an energy gap of the order L~^. On the other hand, we 
do not yet know how to make the AQC model fault tolerant [51]. It will be interesting to 
see whether the first useful quantum computer will be based on the quantum circuit model, 
use measurement-based computation [24], utilize topological quantum computing [55j, or 
whether it will be based on "analog" adiabatic evolution [36j . 

Matrix Product State ansatz based numerics 
— investigating spins on an infinite tree 

It is hard classically to find the spectral properties of (and to simulate the time evolution 
with) AQC Hamiltonians for many qubits. This is why I chose to investigate two simple 
translationally invariant systems. The method I used is based on a local approximate de- 
scription of a quantum state, the Matrix Product State ansatz. I improved the imaginary 
time evolution method for finding the ground state energy of a Hamiltonian within the 
MPS ansatz, and showed how to apply it for a system with an infinite tree geometry of 
interactions. I then investigated the phase transitions in the transverse field Ising model 
and the Not 00 model on the Bethe lattice. The phase transitions I found are not first-order. 
This may have positive implications for the minimum gap scaling of the corresponding AQC 
Hamiltonians on finite trees. While the two models I considered do not encode computa- 
tionally interesting optimization problems, they are an important stepping stone for further 
research. Showing that the systems with infinite tree geometry have similar properties 
as large finite trees, adding randomness and loops to the underlying graph structure, and 
finding the properties of more interesting Hamiltonians are now open for research. 

Matrix Product State methods were thought to work especially well in ID. However, 
after Aharonov et. al. [5] showed that Quantum 2-SAT on a line for particles with di- 
mension d = 12 is QMAi complete, it became clear that MPS methods are not a panacea 
for ID systems, especially for spins with higher dimensions. Moreover, recently Schuch et. 
al. showed that the ground state of a particular frustration free Hamiltonian in ID can 
encode solutions to NP-complete problems, even though the Hamiltonian has a polynomial 
eigenvalue gap and the ground state is exactly described by a low-dimension MPS [84J. 
Classical MPS-based methods for finding this ground state are thus bound to fail. How 
is this possible? The imaginary time evolution method I proposed is not complete. Its 
weakness must lie in the method for bringing the state back to the MPS ansatz, implying 
that the procedure can get stuck in local minima, although the required dimensions of the 
MPS description of the state are small. However, this didn't happen for the two Hamiltoni- 
ans I investigated. DMRG (and MPS) based methods remain a very useful tool for spin-| 
systems in ID and on trees. Nevertheless, new methods such as PEPS are necessary in 
higher dimensional geometry. At this front, exciting new results have been recently found 

^Here weak convergence means that the time averaged probability distribution Pt{x\0), averaged between 
r = and r = to, converges to the limiting distribution 7r(a;|0) as jX^igSH (P'^oi^W ~ '""(^lO))! — where 
SR is a region covering two thirds of the line. On the other hand, strong convergence means that Pt{x\0) 
converges to 7r(a::|0) at every point x, i.e. J2xeA \P-ro{^W ~ '"'(^10)1 < £■ 
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by Verstraete at al. [T021I68] . 



Quantum Satisfiability 

— an interesting analogue of a classical problem 

The complexity of the Local Hamiltonian problem was sorted out by the paper by Oliveira 
and Terhal [74J, who showed that 2-local Hamiltonian on a 2D grid is QMA complete. As 
fc-local Hamiltonian is the quantum analogue of MAX-fc-SAT, it is interesting to look at 
Quantum fc-SAT problem, constructed from a sum of fc-local projectors, which is the ana- 
logue of classical Satisfiability. The solutions to Q-fe-SAT are states which are annihilated 
by all of the terms in the Hamiltonian. This problem seems easier than fc-local Hamiltonian, 
but already for A; = 3 it contains the NP-complete problem 3-SAT. Is Quantum 3-SAT then 
QMAi complete? This big question remains unanswered here. As showed by Cook and 
Levin, the answer to a classical 3-SAT instance can encode the evaluation of a classical 
circuit, because classical information can be freely copied. However, quantum information 
cannot be cloned, forcing us to store the information about a computation in orthogonal 
states such as \tpt) \t)^ by adding a clock register to the workspace. So far, no one has been 
able to show how to encode the clock register in such a way that its transitions would be 
checkable by 3-local projectors while being coupled to the application of a quantum circuit. 
So far, my attempt using projectors on one qutrit and two qubits, the QMAi-complete prob- 
lem (3,2,2)-SAT comes closest to this goal. Wholly new ideas, such as Eldar and Regev's 
triangle construction, which I successfully built on in Section 14.51 are required to answer 
the question of complexity of Q-3-SAT. There I have shown that a Q-3-SAT Hamiltonian 
can be used as a Hamiltonian computer to perform any quantum computation. However, 
the question of the complexity of the original Q-3-SAT problem (as a question about a 
Hamiltonian: is there a zero energy ground state of H) remains open. It might even be 
possible that a classical verifier circuit exists for it. 

In the course of investigating Q-3-SAT, I obtained another result greatly strengthening 
previous results about the Local Hamiltonian problem. The 3-local Hamiltonian was proved 
to be QMA complete in [5l]. When the Hamiltonian in |54j is rescaled to have norm 0(1), 
the separation in the question about its ground state energy (whether Eq is lower than 
a or greater than b) scales like b — a cc L~'^^ . On the other hand, when using my new 
3-local Hamiltonian construction of Section 14.21 I showed that even for b — a oc (and 
||-ff|| = 0(1)), the 3-local Hamiltonian problem remains QMA-complete. If one could show 
this for b — a (X 1, it would result in the quantum analogue of the PCP theorenH, the basis 
of many results about hardness of approximation and a great achievement of theoretical 
computer science [111 [11]. It is unlikely that a method based on perturbation theory gadgets, 
such as one employed in when transforming 3-local Hamiltonian to 2-local Hamiltonian [53], 
can lead to a quantum PCP theorem. After rescaling the Hamiltonian back to norm 0(1), 
each round of perturbation theory necessarily shrinks the gap in the spectrum. Although 
the separation 6 — a is not the gap A = E2 — Ei the Hamiltonian, it is related to it as 
follows. The low energy states in the 'yes' instance of a Local Hamiltonian problem encode 
a correct computation with a positive result. The excited states encode either a wrong 
computation or a negative final result. The eigenvalue separation between these is A. On 
the other hand, in the 'no' instance of Local Hamiltonian, the ground state must again have 

^ PCP (Probabilistically Checkable Proofs): Any 3-SAT instance can be reformulated so that any pur- 
ported proof of the new instance is probabilistically checkable by looking at only a constant number of 
randomly chosen clauses. 
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the computation history wrong or the final result negative. It is very likely that the ground 
state energy is going to be smaller than the gap A in the 'yes' instance. It surely was the 
case for the Hamiltonians we have encountered so far, where A oc L~^^ while the largest 
separation was b — a (x (after rescaling the Hamiltonian to \\H\\ = 0(1)). Is it then 
possible to ever obtain a better scaling of the necessary separation b — a? If it is tied to the 
A as I described, the outlook is bleak. If the ground state of a Hamiltonian with gap scaling 
even like A tx could encode any quantum computation, it would mean we could use the 
Hamiltonian as a computer and obtain the result in time shorter than what is physically 
possible. Thus, when searching for a PCP-like quantum theorem, it seems necessary to find 
new Hamiltonian encodings which do not tie A and 6 — a, or whose underlying idea is not 
that of Feynman's Hamiltonian computer. 

I also showed that Quantum 2-SAT in ID for particles with dimension d = 11 is QMAi 
complete. It would be interesting to look at the problem from the other end and classify the 
hardness of Quantum 2-SAT on a line for lower-dimensional particles, perhaps relating it 
to the matrix consistency problem in ID [62|. My preliminary results show that the ground 
states of random instances of Quantum 2-SAT in ID for particles with dimension d = 4, 
with 4 clauses per spin could be hard to approximate classically. 

Hamiltonian Quantum Cellular Automata 

— universal computing with a simple Hamiltonian 

I proposed a model of quantum computation based on a time-independent, translationally 
invariant Hamiltonian in Chapter [5l The best realization of a HQCA that I found requires 
10 dimensional particles placed on a line with nearest neighbor interactions. I initialize 
the program and data by preparing a computational basis state of the chain, and use a 
diffusion-like mechanism to run the program through the data. What errors could plague 
this model? First, it is prone to an Anderson localization problem in the motion of the 
program, requiring sensitive tuning of the transitions. Other error preventing measures 
(such as encoding an already fault tolerant circuit into the HQCA) can be taken, but the 
model cannot be made fault tolerant, unless the number of 'gate particles' in the system 
is kept fixed. This could be done by using actual particles whose number is conserved. 
However, three kinds of 'gate particles' are needed. Unless this is remedied, this model will 
remain purely theoretical. On the other hand, it shows interesting results about the hardness 
of simulation of simple translationally invariant Hamiltonians. It would be interesting to 
find out how much can the dimension of the particles be lowered, keeping the quantum 
computing universality of the model. One possible new research direction I have embarked 
upon so far is a programmable scheme in 2D, based on 2-local nearest neighbor interactions 
of particles with dimensions 3 and 5. For this, I must increase the size of the unit cell (for 
translational invariance) to about 20 sites. 

There is a lot to be learned when looking at the underlying machinery of quantum 
mechanical systems - local Hamiltonians. They can lead to discoveries of new quantum 
algorithms such as the NAND tree algorithm [32], classical simulation methods such as 
iTEBD |106j . proofs of complexity results such as the QMA completeness of Quantum 
(12,12)-SAT on a line [5j, or a scalable Hamiltonian computer design based on a feasible 
physical Hamiltonian such as the superconducting AQC of [52]. I hope this thesis will serve 
as a pointer to these directions not just for myself, but for the reader as well. 
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Appendix A 

Kitaev's propagation Hamiltonian 



In this Appendix, I look at Kitaev's propagation Hamiltonian (jl.53p in more detail and 
prove a lower bound on its second lowest eigenvalue. 

Consider the system consisting of two registers, the first of which holds n work qubits 
and the second of which holds L + 1 clock qubits: 



— T~^work ^X" 'Hclock- 

Kitaev's propagation Hamiltonian (|1.53|) is 



prop 



Tjt 

prop 



rop^ 
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2 ^pro 
t=0 



where 
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t+l,t 
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is the projector onto the state \t) of the clock register, 
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t+l,t 
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\t + l){t\, 
\t) {t + l\, 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



(A.5) 
(A.6) 



increase and decrease the state of the clock register by one, and Ut acts on two of the work 
qubits. To diagonalize this Hamiltonian, it is convenient to make a basis transform 



\s) fg) \t) 



\^s,t) = {UtUt-i...U2Ui)\s)<E)\t) 



(A.7) 



where |,s) are computational basis states of the work qubit register. In this new basis, Hprop 
(|A.2p has a block diagonal form. 



A 



A 



prop 



A 
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with each (L + 1) x (L + 1) block 



A 



1 _ 1 

2 2 



.1 1 _1 
2 ^ 2 

_ 1 
2 



1 1 
'2 2 



(A.9) 



corresponding to a basis state |,s). Note that the Hamiltonian does not mix sectors with 
different \s). This matrix can be decomposed as 



A = 1-]^{B + Po + Pl)^ 



(A.IO) 



where B is the adjacency matrix for a line of length L + 1 and Pq a-nd are projectors on 
its ends. The eigenvectors of A are plain waves [56], 



I0o) 



1 ^ 

' + 1 



t=0 
L 



L + 



-f^cos(p,(t + ^))|t) 



with momenta 



TTJ 



j = 0,...,L, 



L + 1' 

corresponding to the eigenvalues 

Xj = 1 — cos Pj . 
The smallest eigenvalue is Aq = 0, and the first nonzero eigenvalue is 



Ai = 1 — cos 



vr 



For large L, this eigenvalue is bounded from below by 



(A.ll) 
(A.12) 

(A.13) 
(A.14) 
(A.15) 
(A.16) 



for some constant ci. The eigenvalue gap between the ground state and the first excited 
state of Hprop is thus bounded from below by 



-'--'-prop — J 2 



(A.17) 



The dynamics of this system is a quantum walk on a line with the two additional boundary 
terms. 
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Appendix B 



Dyson Series for an AQC 
Hamiltonian 



Consider a time-dependent AQC Hamiltonian (see Chapter [2]) 

H{s) = (1 - s)Hb + sHp = Hb +1^ {Hp - Hb), (B.l) 

1 ^ ' 

A B 

where s = t/T and Hb does not commute with Hp. Assume also that all the terms inside 
Hb commute with each other, and so do the terms within Hp. My goal in this section is 
to find a formula of the form 

approximating the time evolution operator UtQj^^t.to corresponding to (jB.ip up to order 
(At)^ and possibly higher orders as well. 

The time evolution operator for a time-dependent Hamiltonian H{t) that doesn't com- 
mute at different times can be expressed using the Dyson series 

C/to+At,to=I + y" IT / dh dt2--- dtnH{h)H{t2)...H{t„). (B.3) 

„=1 V / Jto Jto J to 

In the following, I take h = 1. 

The zeroth order term in ()B.3P for H (|B.ip is just Uq = I. The first order term is 

Ui = i-i) I H{h)dh (B.4) 
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and the second order term is 



rto+At / r-ti \ 

U2 = - H{ti)i H{t2)dt2]dh (B.5) 

Jto \Jto / 



Neglecting all terms of order higher than (Af) 7^ in Uq + Ui + U2, I obtain 

V2 = l-i(A*M-i(At)(| + ^)B (B.7) 

2 2T ^ ^ 2T2 

Rewriting A and B back in terms of Hb and //p, I obtain 

V2 = I 



^ ( 1 - ^ ) (AtyiHsHp + FpFb) 



2 



^(^) (At)2i/p//p. (B.8) 



I now want to find a formula of the form (|B.2p . approximating (jB.Sp . The simplest 
version is 

^ {l-iaiHB){l-ibiHp), 



which approximates (jB.SP to order (At) when 



ai = (|)At, (B.IO) 
61 = (l-|)At. (B.ll) 



This agrees with the usual first-order Trotter formula. 
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Next, let me look at 



1 - ia2HB - 1(1- ihHp - ^Hlj I 1 - ic2HB - ^^H^ 



When choosing 



this formula agrees with (jB.SP to order (At)^. It has extra corrections terms, not found in 
the usual second order Trotter formula derived for a time-independent Hamiltonian H 



to' 
T ^ 

hi = ('^)Ai. (B.17) 



4 = cl = fl-ll^, (B.16) 
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Appendix C 



Continuous-time Quantum Walks 
in ID 



First, in Section [C.ll I analyze the continuous-time quantum walk on a line and prove two 
useful lemmas about the mixing of this walk used in Section [5.3.21 and Appendix[Dj Second, 
in Section IC.2I I analyze the quantum walk on a cycle and prove another mixing lemma used 
in the proof of universality of the train switch Hamiltonian Computer in Section 14.51 



C.l Quantum Walk on a Line 

Consider a continuous time quantum walk on a line of length L, where the Hamiltonian is 
the negative of the adjacency matrix for the line 

L-l 

H, = -Y^{\j){j + l\ + \j + l){j\). (C.l) 

The eigenvalues of this Hamiltonian are 

A, = -2cos(-^), (C.2) 

for j = 1 . . . L, while the corresponding eigenvectors |</>'-''^) = Ylk=i 1^) have components 



Consider the time evolution of a particular basis state | c) . The probability of finding the 
system in the basis state |m) at some time r can be found by expanding |c) and \m) in the 
basis of the eigenvectors (IC.3P : 

L 

Pr{m\c) = \{m\ e-^- |c)|' = ^ e-^(^^-^'=)>(^),/.(^-)*0(^)*#. (C.4) 

j,k=i 

Because the time evolution (according to the Schrodinger equation) is unitary, this prob- 
ability Primlc) does not converge. On the other hand, let me define the time average of 
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Primlc) for time < r < T20 as 

Pt-2o('^|c) = — / pT{rn\c)dT. (C5) 
"^20 Jo 

As I will show below in Lemma U this average probability distribution does converge to 
a limiting distribution 7r(m|c), defined as the T20 00 limit of the average probability 
distribution ()C.5p . All the eigenvalues ()C.2p are different, so I can express the limiting 
distribution as 

L 

7r(m|c) = lim Pr2o("^|c) = V j^^^^j j^-^^l , (C.6) 

Ton — ^ryT] ' ^ ' ' ' ' 

which in this case is 



r20— +00 



-(-|c) = . (C.7) 

According to the following lemma, the average probability distribution (jC.SP converges 
to the limiting distribution 7r(m|c). 



Lemma 5. Consider a continuous time quantum walk on a line of length L, where the 
Hamiltonian is the negative of the adjacency matrix for the line. Let the system evolve for 
time T < T20 chosen uniformly at random, starting in a position basis state \ c) . The average 
probability distribution Pt-2o{-\c) converges to the limiting probability distribution it{-\c) as 

\Pr20im\c) - 7r(m|c)| < O (—) . (C.8) 

Proof. First, recall Lemma 4.3 of [1] for the total variation distance of the probability 
distribution ^^-20 from the limiting distribution, saying 

V|Pr2o("^|c) -vr(m|c)| < — V ' [ . (C.9) 

Using ()C.2p and ()C.3p . I can bound the expression on the right of ()C.9p . When j is close 
to k, i.e. \j — k\ < Ci, I can write 



/,(i)l2 



On the other hand, for \j — k\ > Ci 1 can bound 

U(i)|2 



< 2. (C.IO) 



Ao — At. L -\- \ 



< (C.ll) 
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with Ci and C2 constants independent of L. Inserting into ()C.9p . I have 



m=l 



m\c) — 7r(m|c)| < 1 = O — 

T20 ^"20 



(C.12) 



□ 



Using Lemma [5l I now prove a useful result utilized in the time analysis of the d = 20 
HQCA in Section dO]) . 



Lemma 4. Consider a continuous time quantum walk on a line of length L, where the 
Hamiltonian is the negative of the adjacency matrix for the line. Let the system evolve 
for a time r < T20 chosen uniformly at random, starting in a position basis state \c) . 
The probability to measure a state \t) with t > L/6 is then bounded from below as p2o > 

l-of— 

6 \T20 

Proof. The probability to measure a state \t) with t > L/6 at time r < T20 chosen uniformly 
at random is 



P20 



m\c). 



(C.13) 



Starting with (IC.8|1 . I have 
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(C.14) 
(C.15) 

(C.16) 
(C.17) 



Therefore, the probability of finding the chain in state \4't>L/Q) s-t ^ random time r < T20 
is thus bounded from below by 



5 ^ L 
P20 > 7T - O — 

6 \T2Q 



(C.18) 
□ 



Also, Lemma H] can be easily generalized for any desired probability q. Measuring a 
state which started in a position basis state |c) at a random time t < Tq chosen uniformly 
at random, the probability to measure a state \t) with t > {1 — q)L is then bounded from 
below hy p>q-0 . 
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C.2 Quantum Walk on a Circle 



If the geometry of the system is a closed loop of length L instead of a line, the Hamiltonian 
(|C.ip gets an additional wrap-around term. 

L-l 

Hioop = - m (1| + |1) (L|) - (li) (i + 1| + \j + 1) 01) • (C-19) 

The eigenvalues of this Hamiltonian are 

Aj = -2cos(pj) , (C.20) 
corresponding to plain waves with momenta 

= ^, (C.21) 
for j = ... L — 1. The corresponding eigenvectors \4>^^'^) = Ylk=i ^k'^ 1^) have components 

rf^ = (C.22) 
= i = l,...,L-l, (C.23) 

These can be combined to make real eigenvectors. For my analysis, it will be enough to 
consider a line with even length, and only the cosine plain waves: 



y^cos(ipj/c), j = l,...,^. (C.25) 



The limiting distribution on the cycle when starting from site c is 

7r(m|c) = ^-^- (C.26) 

for all points m except for m = c (return back) and m = L + 1 — c (the point across the 
cycle), where I have 

Am\c) = \-^- (C.27) 

As in Section IC.ll I will again utilize Lemma 4.3 of [3] to prove the convergence of the 
time-averaged probability distribution to this limiting distribution. On the right side of 
dCH), I now have 
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Figure C-1: A contour plot of |Aj — A^l ^ with Xj from (|C.20p . The sum in ()C.3ip is over 
the marked region. 



The sum over the non-equal eigenvalues 



(C.29) 



is now more complicated, because of the degeneracy of the spectrum. I plot |Aj — Xk\~ in 
Figure IC-li Because of the symmetrical way Xj ()C.20p arise, assuming L is divisible by 4, 
I can express ()C.29P as 



L/4 L/2-j 
j=0 k=j+l 

L/2 L/AL/2-j 



< 



j=l k=j+l 



Xj — Xk\ 



(C.30) 
(C.31) 
(C.32) 



5k 



The term Aq comes from j = 0. I can bound the sum by an integral, taking x = ^ 
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obtaining 



L/2 



A 



fc=i 



For X G [O, i] , I can bound 



L 1-2 
< — 



cos 



dx 



27rf ) ~ 2 7^ 1 - cos(27rx) ' 



1 — cos 27r.T > 27rx , 



resulting in 



Ao<^ r% = ^{L'-2L)=0{L'). 
Att Ji Att ^ ' 



(C.33) 



(C.34) 



(C.35) 



I bound the other term, Ai, in ()C.32p by an integral as well: 

L/4 L/2-j 

^ j=i k=j+i cos (27r|) - cos (^27ri 

< ll. f^'^ dy — 

~ '2 J ^ Jj^+y COS [lirx) — cos (27ry) 

Again, I can lower bound the difference in eigenvalues for y G [O, ^] and x G [y, ^ 

cos (27rx) — cos (27ry) > 27r(x^ — y^), 

allowing me to write 

L2 



■ 4 I, ^ /■ 2 f dx 



r, 11 

An Ji ^ 



2y 



i dy 



— / [log (1 - Ay) + log (1 + 2yL)] 
An Ji_ 2y v , ' 



(C.36) 

(C.37) 

y] by 
(C.38) 

(C.39) 
(C.40) 
(C.41) 



As y in ()C.4ip is at most | — j;, I can bound R by 

<logL. 

Finally, this results in 

L^logL f^^T: dy L^logL 



^1 < 



2n 



2y 



2n 



log 



4 L 



< 0(L2 log^L). 



(C.42) 



(C.43) 



<1okL 
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Putting (10351) and (lOiSll into ([032]) . I obtain 



So<OiL^log^L). 
Lemma 4.3 of [4] (see (IC.9P ) then reads 



E 



2 1 

\pr2o{m\c) - 7r(m|c)| < — -O(L^log^L) = O 

To 



L log^ L 



(C.44) 



(C.45) 



Thus, for Tcirc = eO(Llog^L), the time-averaged distribution becomes e-close to the hm- 
iting distribution. Using the bound on the total variation distance I just proved, it is 
straightforward to obtain the following lemma which I use in Section 14. 5t 

Lemma 6. Consider a continuous time quantum walk on a cycle of length L (divisible by 4), 
where the Hamiltonian is the negative of the adjacency matrix for the cycle. Let the system 
evolve for a time r < Tq chosen uniformly at random, starting in a position basis state |0). 
The probability to measure a position state \t) farther than L/6 from the starting point (the 



farther two thirds of the cycle ) is then bounded from below as po > 



3L 



O 



L log^ L 



Proof. The proof is analogous to the proof of Lemma [5] in Section IC.ll Let me call the 
farther two thirds of the cycle (see Figure 14-111 in Section 14. Sp the success region (SR). 
When I choose the time r < Tq uniformly at random, the probability to measure a state \t) 
with t S SR is 



(C.46) 



m&SR 



Using the bound on the total variation distance (|C.45P I just proved and the formulae for 
the hmiting distribution (10261) . ([0271) . I have 



O 



L log^ L 



> 



> 



> 



J2 \Pro{'m\c) - 7r(m|c)| 
\pTo{m\c) 



m=l 

E 



vr m c 



E 

Po 



Pra{m\c) 



TT m C 



(C.47) 
(C.48) 

(C.49) 
(C.50) 



Therefore, the probability of finding the chain in state {ipteSR) at a random time r < To is 
thus bounded from below by 



2 1 



L log^ L 



(C.51) 



□ 



It is thus enough to wait a random time not larger than 0(Llog^ L) to find the state of 
the system in the success region with probability close to two thirds. 
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Appendix D 

Diffusion of Fermions on a Line 



Here I prove Lemma [3l a result about the mixing of a discrete free fermion gas, used in 
Section 15.21 

Lemma 3. Consider the state 

l^o) = 4m+i^/m+2 • • • 4m+2M |0) ■ (D-1) 

of 2M fermions on the right end of a line with L = [f + 2)M sites. Let the system evolve 
for a time chosen uniformly at random between and tiq with the Hamiltonian 

L-l 

Hf = -J2 b]bj+i + h.c. (D.2) 
i=i 

and measure the number of fermions in the region 1 < a; < fM. The probability to measure 
a number greater than M is piQ > — O • 

Proof. Let us start with the outhne of the proof. We look at the fermionic system in both 
first and second quantization to obtain an expression for the time evolution of the creation 
and annihilation operators in the Heisenberg picture, mapping it to a quantum walk on a 
line. We then consider the observable X, the number of particles sufficiently far from the 
right end of the line. We will show that when we choose the time to measure X uniformly 
at random between and tiq, the expected value we will obtain is approaching a number 
close to 2M. To show this, we will express the expected value of X in the time-averaged 
state of the system using the results from a quantum walk on a line. Finally, because the 
number of particles in the system is 2M, we will deduce that the probability to measure a 
number less than M is then small. 

Observe that Hf is the Hamiltonian of a free fermion gas on a line in second quantization 
(a special case of the XY model). The time evolution of the state |^o) can be obtained by 
looking at the problem back in the first quantization, where we write l^'o) ^ 



\(pfM+l) ® \<PfM+2) ® • • • (8) \4>fM+2M) , (D-3) 



with \(f)j) = \j) in the position basis and [•] the standard antisymmetrization operator. 
We first solve for the time evolution of the corresponding one-particle wavefunction \<pj{T)) 
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with the Hamiltonian 



L-l 



Hi = -Y.{\j) (i + ii + ij + i) 01), 



(D.4) 



and then obtain the solutions for the many-particle problem by antisymmetrization as 



(D.5) 



The eigenfunctions of Hi (quantum walk on a line) are plain waves (as in ()C.2p and (jC.SP ). 
and the time evolved states {(pjir)) thus readily available. Let us define the unitary matrix 
u{t) by 



/W)=£Mr)|A;) 



(D.6) 



k=l 



Returning to the second quantized system, the time evolution of the creation and annihila- 
tion operators in the Heisenberg picture is then 



(D.7) 



k=l 



k=l 



Consider now the observable X, the number of particles in the first fM sites of the line 
with length L = {f + 2)M 



Its expectation value at time r is 



fM 

m=l 



fM 



(Di 



EM) = Y.{^>{T)\nm\^{T)). 



(D.c 



m=l 



The number operator for site m is hm = h'mhm- We can go to the Heisenberg picture and 
use ()D.7p to write 



{■^{T)\hm\^{T)) = {■^^\bl{T)bm{T)\^^) 



L L 



Yl ^"^c{T)u*^dir) (*o| blbd l^o) 

c=l d=l 
L 



c=l 



\Umc{r)[ 
•' PT(m\c) 



(D.IO) 
(D.ll) 

(D.12) 

(D.13) 
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where each term \umc{T)\'^ = Prin^lc} can be thought of as the probabihty of finding a 
particle at site m at time r when it started from the site c and performed a quantum walk 
on a hne, according to (|D.4p . Inserting this into ()D.9p . the expected number of particles 
not in the rightmost part of the chain at time r is 

L / fM \ 
EriX)= [Y.Pr{m\c)\. (D.14) 

c=/M+l \m=l J 



Let us now choose the time r uniformly at random between and tiq. The average 
value of X (the expectation value in the time-average state) is 

Er,,{X) = — I Er{X)dT. (D.15) 



no Jo 

For a quantum walk on a line, the time-averaged probability ()C.5P of finding a particle 
that started at position c at final position m converges to the limiting distribution (|C.7|) 
according to Lemma 1 (jC.SP proven in Appendix [Cj Using this fact, we can show that the 
expectation value E'tio {X) in the time-averaged state converges to the limiting expectation 
value 

m<fMc>fM 

as 

\E,^^{X)-E{X)\<o{^y (D.17) 

Recalling the limiting probability distribution for a quantum walk on a line of length L 
()C.7p . we have 

m</A/c>/A/ 

2 1 
= /Mx2Mx— --h2Mx— (D.19) 

= (D.20) 

Putting this into ([DT7|) . the average value of X when the time r < rio is chosen uniformly 
at random is bounded from below as 



f + 2) Vno 



Er,M) > 2M ( ^ ) - O ( ^ ) . (D.21) 



We want to find the probability of measuring X > M. First, the maximum possible value 
we could measure at any time is 2M, the number of particles in the system. Second, the 
average value Et-^o^X) at time r chosen randomly is close to 2M. Therefore, the fraction 
A of times at which we measure a number significantly lower than 2M must be small. Let 
us bound A in the worst case scenario. This is when each unsuccessful measurement yields 
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X = M, and each successful measurement gives us 2M. We then have 

AM+(1-A)2M > ^rio(^), (D-22) 
A < — . (D.23) 

Hence we arrive at the desired bound on the probability to measure X > M: 

□ 
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Appendix E 

Transitions in a two level system 



In this Appendix I give the analysis of transitions in a simple two-level system referenced 
in Section 12.3.11 Consider a two level system with Hamiltonian 



His) = Eois) \Ms)) {Ms)\ + E,is) \Ms)) {Ms)\ 



(E.l) 



which varies smoothly with s = t/T. Here \(f)o{s)) and |0i(s)) are orthonormal for all s. 
The Schrodinger equation reads 



ds 



TH{s) 



(E.2) 



The two energy levels in the system are separated by a gap 

g{s)=Ei{s)-Eo{s), (E.3) 
which I assume is always larger than 0. Let me introduce 6 (with the dimension of energy) 



as 



and let 



e{s)= / g{s')ds', 



(E.4) 



I pick the phases of |0i(s)) and \(j)o{s)) such that ((/'i(s)| ^ {(piis)) = (</'o(s)| ^ I0o(s)) = 0. 



Plugging IE.5I into the Schrodinger equation gives 



dcQ 

ds 
dci 

d7 



^ JT0 



dl 
d_ 

d^ 



or equivalently, 



dco 
d^ 
dci 

d^ 



(E.6) 
(E.7) 

(E.8) 
(E.9) 
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where 



fiO) 



d 



9 



dH 



de 



dH 



ds 



(E.IO) 



Now let 0{1) = 0.1 started with ci(0) = and I want to find the transition amplitude at 
s = 1 which is 



coe'^'fde 



-Co/ 



1 
T 



iT 

icofe 



J 

re 



de 



'dd 



iTe^ df 



Now I Co I < 1 and |ci| < 1. As long as the gap docs not vanish |/(^)| and 

therefore |ci(^)| = O (^). The probability of transition to the excited state 
system with a nonzero gap is thus 



(E.ll) 
(E.12) 
(E.13) 

are bounded, 
:"or a two-level 



\c,{ef = o{T- 



(E.14) 
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